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ABSTRACT 


A multi-frequency quadrature phase shift keyed (MFQPSK) signal has 
been developed at NPS to be used in computer-to-computer communications. 

This report discusses the simulation of a MFQPSK signal transmitted 
from a moving transmitter platform through a near vertical acoustic channel 
as seen by a moored receiver. The simulated received signal is tested against 
MF QPSK signal theory. The simulation was developed to be an experimental 
tool for testing various Doppler, synchronization, and coding algorithms/ tech- 
niques for a MF QPSK communication signal. The degradation of output 
signal-to-noise ratio due to Doppler shifts caused by the moving transmitter 
is analyzed. An algorithm for estimating Doppler compression/expansion in 


the received signal is evaluated. 
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I. INTRODUCTION 


Data links used in many modern military communication systems re- 
quire computer-to-computer communications which transmit bits from one 
terminal to another with a low error rate. The communication channel may 
be a lowpass channel, if the transmission medium is wire or optical fiber, or a 
bandpass channel, if the transmission medium is radio frequency (RF) link or 
acoustic channel. The signal used to carry the data bits must be a properly 
modulated signal with sufficient energy. It must also be positioned in the fre- 
quency spectrum to propagate effectively through the intended medium. 
Multi-Frequency Modulation (MFM) is a signal modulation format that is 
readily adaptable to a variety of data link scenarios. MFM does not require 
special purpose MODEMS to translate between the digital and analog do- 
mains and can emulate most existing signal modulation formats as well as 
generate entirely new formats. MF'M’s descriptive language is that of Digital 
Signal Processing (DSP) [Ref. 1]. Simulation and analysis of a Naval Post- 
graduate School (NPS)-developed MFM signal received through a bandpass 


channel is the subject of this thesis. 


A. BACKGROUND 

Classical modulation methods for bandpass channels use amplitude and/ 
or phase to carry signal information on a carrier wave in the channel. When 
the information source is a finite alphabet, as with data or quantized analog 


sources such as digitized speech or video, then only a finite number of signal 


states are needed to represent or code the source. Phase shift keying (PSK) is 
an efficient method for coding these signal states or message [Ref 2]. 

MFM is a modulation technique that is ideally suited for computer-to- 
computer data transmission and reception because its basic structure is one 
of time and frequency slots. In MFM, the signals are directly encoded, modu- 
lated, decoded, and demodulated using DSP techniques within the host com- 
puter. In MFM, signal “packets” located in a given frequency spectrum and 
time are created as shown in Figure 1. Figure 1 shows a single signal packet. 
Each packet consists of L bauds of K tones. The length of a baud is AT seconds 
during which K discrete tones are transmitted. The phase of each tone is the 
coded information. These LK subsignals form an orthogonal signal set. Each 
subsignal may be independently modulated with phase information. 

Dr. P. H. Moose at NPS, in conjunction with the Naval Ocean Systems 
Center (NOSC) in San Diego, has developed Multi-Frequency Quadrature 
Phase Shift Keying (MFQPSK) signals to be used in high data-rate acoustic 
burst communications from moving platforms using medium frequency, rela- 
tively wideband (25 percent bandwidths) links [Ref. 1, 3]. A simulation of the 
MF QPSK signal from a moving platform as seen through a bandpass channel 
was needed to analyze and test various Doppler, synchronization, and coding 
techniques and/or algorithms before final implementation of the modulation 
scheme is actually realized. Thus, the subject of this thesis is the simulation 
and analysis of a MFQPSK signal from a moving platform as seen through a 


bandpass channel. 
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Figure 1. MFM Signal Packet (after Ref. 1: p. 3.) 


B. OVERVIEW 

The NPS-designed MFQPSK signal will be transmitted from a moving 
platform. In Chapter II a simple model of the moving transmitter platform’s 
dynamics is developed; its performance is compared to actual track data. This 
model is implemented in the simulation. 

A mathematical description of the transmitted signal and the design 
parameters of the MF QPSK signal are introduced in Chapter III. The channel 
parameters which delay and compress/expand this transmitted signal and 


which characterize the received signal, the additive Gaussian noise which is 


added to the received signal after sampling, and the final simulated sampled 
received signal are derived and presented in Chapter ITI. 

Once the simulation was complete, testing and analysis was performed. 
First, the simulation was tested against MFQPSK signal theory to insure 
that the simulated received MFQPSK signal is consistent with theory. The 
simulated received signal should be consistent with theory so that the simu- 
lation may be used as an experimental tool for testing various Doppler, syn- 
chronization, and coding algorithms. Chapter IV consists of this analysis and 
simulation results. 

Second, the amount of residual Doppler mismatch which the signal can 
tolerate due to the moving transmitter was analyzed. Chapter V provides the 
analysis of the output signal-to-noise ratio degradation due to the Doppler 
mismatch. Then, in Chapter VI, an estimation of the Doppler compression 
factor within a baud is derived using the Discrete Fourier Transform of a re- 
ceived signal’s baud. The simulated outcome of this estimated Doppler versus 
the Doppler due to the moving transmitter within a baud is also illustrated in 
Chapter VI. 

In the simulation, the actual Doppler compression/expansion factor due 
to the dynamics of the moving transmitter is computed; therefore, it is 
known. The start time of the signal reception is also known; thus, the signal 


is always perfectly synchronized in the simulation. 


II. MOVING TRANSMITTER PLATFORM DYNAMICS AND 
BEAM PATTERN 


A. BACKGROUND 

Since it is assumed that the MF QPSK signal will be transmitted from a 
moving platform, provision must be made at the receiver to compensate for 
Doppler shifts. The Doppler shift factor is computed based on the geometry in 
the simulation using the radial velocity of the moving transmitter relative to 
the stationary receiver. The complete set of equations for the Doppler shift 
factor is presented in the next chapter. The radial velocity is the derivative 
with respect to time of the slant range to the receiver. Realistic slant range 
values are necessary to compute realistic Doppler shift factors. Slant range is 
computed using the time-varying x-position, y-position, and z-position of the 


transmitter relative to the receiver. 


B. ACTUAL TRACK DATA VS. SIMULATED TRACK DATA 

To realistically model the dynamics of the moving transmitter, actual 
track data was provided by NOSC. The parameters that were given in the 
actual track data were x-position, y-position, z-position, and slant range rela- 
tive to the receiver. A simple random walk model was used to model the track 
data. This model assumes the transmitter is moving along a straight line 
with random fluctuations in the x, y, and z velocity components. Figures 2 
through 7 illustrate how well the random walk model corresponds to the 
actual track data on six different runs. The z-position seems to consistently 


have the largest error. The largest error in the z-position is 25 feet, (see 
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Figure 2a. Run 1: Actual and Simulated Track Data 
X—Position and Y—Position vs. Time 
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Figure 2b. Run 1: Actual and Simulated Track Data 
Z—Position and Slant Range vs. Time 
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Figure 38a. Run 2: Actual and Simulated Track Data 
A—Position and Y—Position vs. Time 
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Figure 3b. Run 2: Actual and Simulated Track Data 
Z—Position and Slant Range vs. Time 
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Figure 4a. Run 3: Actual and Simulated Track Data 
X—Position and Y—Position vs. Time 
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Figure 5a. Run 4: Actual and Simulated Track Data 
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Figure 5b. Run 4: Actual and Simulated Track Data 
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Run 5: Actual and Simulated Track Data 
X—Position and Y—Position vs. Time 
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Figure 6b). Note that the z-position graphs have a smaller vertical scale than 
the x- and y-position graphs. Ideally, since the transmitter is attempting to 
maintain a constant depth, the location should not fluctuate much in the 
z direction as was the case in Run 5 shown in Figure 6b. However, overall 
there is little error for these six runs in the model output for the slant range; 
therefore, the random walk model was considered to be an adequate model of 
the transmitter’s dynamics for the purpose of simulating Doppler shifts, time 


delays, and propagation losses of the MFM acoustic communication signal. 


C. TRANSMITTER PLATFORM DYNAMICS MODEL 

The following three equations are the random walk model used to simu- 
late the transmitter’s dynamics. The names of the variables used in the simu- 
lation are the same as those used in the equations below. When applicable, 
the parameters from Figure 1, which are printed in boldface, are equated to 


their simulation variables, printed in capital letters. The dynamics model is 


described by: 


X(LL) = X(LL-1) + (VX(LL) * DELT(LL)) (1) 
Y(LL) = YCLL-1) + (VY(LL) * DELT(LL)) (2) 
Z(LL) = Z(LL-1) + (VZ(LL) * DELT(LL)) . (3) 


The slant range to the receiver 1s computed using the following equation; 


R(LL) = ( X(LL)2 + Y(LL)? + Z(LL)4 )9-5 (4) 
where, 


X(LL): The position of the transmitter relative to the receiver in the x 
direction during the LLth baud 


Y(LL): The position of the transmitter relative to the receiver in the y 
direction during the LLth baud 
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Z(LL): The position of the transmitter relative to the receiver in the z 
direction during the LLth baud 


R(LL): Slant range of the transmitter relative to the receiver during 
the LLth baud 


LL=1: The baud number of the transmitted signal 


DELT(LL)=AT: The LLth baud length in seconds 

VX(LL): The velocity in the x direction during time DELT(LL). VX(LL) is 
Gaussian distributed with mean = VXAVG and variance = 
VXVAR. VXAVG and VXVAR are input by the user. 

VY(LL): The velocity in the y direction during time DELT(LL). VY(LL) is 
Gaussian distributed with mean = VYAVG and variance = 
VYVAR. VYAVG and VYVAR are input by the user. 


VZ(LL): The velocity in the z direction during time DELT(LL). VZ(LL) is 
Gaussian distributed with mean = VZAVG and variance = 
VZVAR. VXAVG and VXVAR are input by the user. 


Note: LL and 1 are used interchangably throughout. 


The rectangular coordinate system of the transmitter’s motion is relative 
to the receiver. The origin (0, 0, 0) is perpendicular to the receiver in the 


plane of the transmitter (the x-y plane) as shown in Figure 8. 


D. TRANSMITTER BEAM PATTERN 

If the receiver is not within the transmitter’s transmission beam, then 
the message will not be received. The simulation does not use the actual 
beam pattern of the transmitter; therefore, it does not know if the receiver is 
within the beam pattern. The simulation assumes the receiver is within the 
transmitter’s beam the entire time of transmission (i.e., from the beginning of 


the first signal packet to the end of the last signal packet). 
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Figure 8. Coordinate System of the Transmitter and Receiver 


The angle between the transmitter’s slant range and z-position relative 
to the receiver is computed and output graphically by the simulation. This 


angle, THETAD(LL) in the simulation, is 
THETAD(LL) = arccos( abs(Z(LL)) / R(CLL) ) * (180.0 / PI) (5) 


where PI = 3.141592654. Figure 9 illustrates the geometry of the line of sight 
angle, THETAD(LL) = 0, and THETAO = @, the half beam width of the 
transmission beam. 

The simulation user can graphically observe the value of THETAD(LL) 
for each baud and determine for a given THETAO if the receiver would 


actually be able to receive the transmitted signal. The simulation does not 
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stop if THETAD(LL) is greater than the input half beam width of the trans- 
mitter, THETAO, but a warning message is output to the screen. The initial 


position (X0, YO, ZO) and THETAO are inputs to the simulation. 
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Figure 9. Geometry of the Transmission Beam 


Al 


IlY. DERIVATION OF THE SIMULATED RECEIVED 
MFQPSK SIGNAL 

A. DESCRIPTION OF THE TRANSMITTED MFQPSK SIGNAL 

Referring to Figure 1, the following definitions are used in MFQPSK 
[Ref. 1]: 

T: Packet length in seconds 

AT : Baud length in seconds 

k, : Baud length in number of samples (not in Figure 1) 

L: Number of bauds per signal packet 

At: Time between samples in seconds 

f= 1/At: Sampling frequency in Hz 

Af=1/AT: Minimum frequency spacing between MFM tones in Hz 

K: Number of MFM tones in a baud 

1: Baud number 


k: Harmonic number of the MFM tone 
®), : Symbol or phase on the kth tone of the 1th baud 


Since At = AT /k,, the sampling frequency is fy = ky * Af. Consequently 
there are a maximum of k,/2 tones spaced Af Hz apart in frequency covering 
the range from 0 Hz to f,/2 Hz. Here f,/2 is the Nyquist frequency [Ref. 2]. 
The K tones carry the phase information during each baud. The relationship 
between fj, the ith frequency, and k;, the ith harmonic, is k; = f; / Af. Some of 
the tones may not be used (1.e., their amplitudes are set equal to zero) during 


any or all bauds of the packet. To generate the transmitted bandpass signals 
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between frequencies f) and fg, only tones between harmonics kj, = f; / Af, the 
minimum harmonic number or tone, and kg = fo / Af, the maximum tone, will 
be given non-zero amplitudes. The harmonic numbers less than kj and 
greater than kg are given zero amplitudes. The number of tones in a baud is 
K = ko - ki} + 1. These K contiguous tones are transmitted with non-zero am- 
plitudes. The signal bandwidth is W = K * Af. Thus, the time bandwidth 
product of the entire signal packet is TW = L * AT * Af * K = LK, which is the 
total number of symbols that can be sent in one signal packet. 

A mathematical description of the transmitted MF QPSK signal is neces- 
sary to understand some of the parameters used to simulate the received 


MF QPSK signal. The 1th baud of the transmitted signal is described by: 


X) (u)= > Xn a a0 = AT (6) 


where x};,(u) = Aj, cos( 2m k Afu + ®), ) [Ref. 1]. Here, u is time referenced to 


the beginning of the baud. Actual real time is t = t, + (I*AT) + u where t, is 


the time at the beginning of the 0th baud (i.e., the beginning time of the first 
signal packet). 


The discrete time signal corresponding to the lth baud is generated by 
sampling (6) at the sampling intervals At = 1/f,. Thus, the discrete time signal 


is given by: 


x) Gay) = Sy Xn, (Hie oO =n < (k -1) Cr) 
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where x]k(n) = Aj; cos( (2a k n)/ ky + ®), ). Here, n is discrete time refer- 
enced to the beginning of the baud. 

Note that a baud interval of time AT seconds contains exactly k cycles of 
tone k. Therefore, adjacent tones differ by one in the number of cycles they 
generate during a baud. 

The phase, ®),, may be given one of the four values 7/4, 37/4, -3n/4, or 
-t/4, which are in quadrants 1 through 4 respectively. These four values make 
up the symbol set used to code the information or message on the signal. The 
simulation variable for ®), is PHI(LL,K). 

The design parameters used in the MFQPSK signal designed and devel- 
oped by Dr. P. H. Moose at NPS are listed in Table I. These parameters are 
for a signal packet in a 16 to 20 KHz bandpass channel. These parameters 
are also used in the simulation to uniquely characterize a chosen baud type 
Gi.e., baud types 1 through 5). The simulation variables in capital letters are 
equated to the parameters they represent. Recall that LL = 1, the baud 
number. 

B. TIME DELAY AND COMPRESSION/EXPANSION OF THE 

TRANSMITTED SIGNAL 

When the MFQPSK signal is transmitted, there are various factors 
which can affect the signal such as a moving transmitter, the channel or the 
medium, and the receiver. In the simulation, these parameters are computed 
and applied to the transmitted signal, thus producing a model of the received 
signal baud. By delaying and compressing/expanding the signal in time, the 


frequencies of the transmittedsignal are shifted in the received signal. Some 
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TABLE I 


DESIGN PARAMETERS FOR SIGNAL PACKET 
IN A 16-20KHZ BANDPASS CHANNEL 


BAUD TYPE: 1 2 3 a. D 


AT = DELT(LL) 1/240 1/120 1/60 1/30 1/15 
(sec) 

Af = DELF(LL) 240 120 60 30 15 
(Hz) 

ky = KMIN(LL) 68 135 Zoo oof 1073 

f, (Hz) 16320 16200 16140 16110 16095 

ko = KMAX(LL) 83 166 332 664 1328 

fo (Hz) 19920 19920 19920 19920 TIZO 

ky = KX(LL) 206 212 1024 2048 4096 

fy (Hz) 61440 61440 61440 61440 61440 


of the parameters or variables, which delay and compress or expand the 


transmitted signal in time, are illustrated in Figure 10 for the 1th baud. 
Recall x;(u) is the 1th baud of the transmitted signal, which is described 


by (6). In Figure 10, x;(u) and y)(u) are represented with a rectangle for illus- 
tration purposes only. Recall x;(u) is the superposition of x},(u), for kj<k<ko, 


over a total of K tones. The frequency of xj,(u) is 


O, = M, * (k/k,) (8) 
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Figure 10. The 1th Baud of the Transmitted Signal Delayed and 
Dilated in Time 


where , = 2x f, and f, is the sampling frequency of the transmitted signal. 


Table I lists the values of these parameters and the applicable variables 
which are used in the simulation. 

When x};,(u) propagates through the transmitter system, it may be 
delayed by the transmitter system electronics. In the simulation, this delay is 


denoted by 7, = TAUK(K), and is a constant for all k. (In the simulation, 


TAUK(K) may be set to a constant or to any known function of k.) Acoustic 
transmission of x;(u) begins at u = 7, and ends at u = AT + TJ. 

When the signal, x)(u - T,), leaves the moving transmitter, it is delayed 
and compressed or expanded in time due to the movement of the transmitter 


and the distance it travels through the channel to the receiver. (The time 
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delay introduced by the electronics of the receiver has not been included in 
the simulation.) This time delay and compression/expansion factor will be re- 
ferred to as £(u). Derivation of £(u) is presented in Appendix A. The time 


delay and compression/expansion factor due to the medium is described by 


£(u) = (T + au) /(1 + gy) (9) 
where 


T= The time for the 1th baud of the transmitted signal to travel through 
the channel and arrive at the output of the receiver. 


a, = The Doppler compression factor due to the moving transmitter 

platform. 

The parameters, T) and a), depend on the slant range to the receiver, 
R(LL), and the speed of sound of the acoustic channel, C(LL), for the 1t® baud. 
Recall R(LL) is computed using equation (4). The speed C(LL) is Gaussian 
distributed with mean = CO and variance = CVAR. The parameters CO and 
CVAR are input by the simulation user. The parameters, T,; and q), are 


related to R(LL) and C(LL) in the following manner: 


Ces lene) — RCI) GUL) (10) 


and ; 
or= ALPHAG.L) =RGL)/ CLL) (11) 


where 
{(X(LL) * VXAVG) + (Y(LL) * VYAVG) + (Z(LL) * VZAVG)] 


R(LL) = 
R(LL) (12) 


Recall X(LL), VXAVG, Y(LL), VYAVG, Z(LL), and VZAVG were described in 
the previous chapter with equations (1) through (3). TAUL(LL) and 
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ALPHA(LL) are the simulation variables for T; and a), respectively. When the 
transmitted signal, x)(u - 7), arrives at the receiver output, it is compressed/ 
expanded in time by £ (uw). 


The received signal is 
k 


2 
y) (a) D3 (u) ; Hesse 


t (13) 


where 


y}, (a) = Fy, Ay, cos(2n k Af (u - T, - £(u)) + Dy), (14) 


and F);, is the attenuation of tone k due to the propagation through the 
channel. The signal, y)(u), will begin at time u, and end at time ug. The 
signal, x)(0), transmitted at time u = 0 is the same signal, y)(u}), that arrives 
at the receiver at time uj; and likewise, the signal, x,)(AT), transmitted at time 
u = AT is the same signal, y)(ug), which arrives at the receiver at time ug (see 
Figure 10). Since it is the same signal, the times that this signal arrives at 
the output of the receiver may be equated to the corresponding times that 


this signal was sent to yield 


uy = (| + Q)) hy + 7 (15) 


and 


ug = (1+ a) AT+#+(L+a)t+T . (16) 


It is obvious that y)(u) is also a superposition of signals, y),(u). However, 


y},(u) is at frequency oy,', where 


ay,’ = @, / (1 + a) ea 
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due to the moving transmitter. At the receiver, y;(u) will be sampled at 


My = My / (1 + Op) (18) 


where @,, is the Doppler factor associated with the mth Doppler channel. If 
Qt) = O,,, then the frequencies of the received signal are exactly in the center of 
Doppler channel m. However, in general, a will fall between Doppler 
channels and there will be some residual Doppler mismatch a - a,. If the 
channels are spaced at intervals Aa, then the maximum mismatch will be + 
Aa/2. In the following section, a derivation for the spacing of the Doppler 


channels Aa is presented. 


C. DOPPLER ESTIMATION AND DOPPLER CHANNELS 

Several processing approaches were considered for estimating the chan- 
nel model parameters. However, the ultimate goal is to send a signal through 
a bandpass channel using MFQPSK modulation on a succession of bauds 
which constitute a signal "packet". In the presence of white noise, the opti- 
mum receivers for these signals are filters matched to each of the tone/phase 
combinations. Recall the four phases are {n/4, 37/4, -3n/4, -1/4} for a MFQPSK 
signal. Therefore, a filter matched to the transmitted signal with 0), = 1/4 
will have a positive output when 7/4 is sent, a negative output when -37/4 is 
sent, and an output of zero when 37/4 or -z/4 is sent, provided that the 
receiver is synchronized with the signal. Similarly, the output of a filter 
matched to the transmitted signal with phase equal to 3n/4 is positive when 


37/4 is sent, negative when -7/4 is sent, and zero otherwise. Also filter pairs 
matched to the ith frequency of the transmitted signal, f; = k;Af, produce zero 


outputs for all phases at the frequency f; = kjAf when k; # k;. 
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To summarize, a matched filter system for MFQPSK signals with K 
tones and baud length AT = 1/Af consists of K pairs of filters, one filter of each 
pair matched to phase z/4 and the other filter of the pair matched to phase 
37/4. Each filter pair output demodulates the phase information encoded on a 
particular tone. 

Since the maximum Doppler shift will be present on the highest fre- 
quency or tone in the baud, the response of the filter pair matched to this tone 
will be analyzed to show how the Doppler channel spacing, Aa, is derived and 
how a, is computed in the simulation. 


Let hj(u) and h,(u) be the filters matched to the signals: 


XIko,i(U) = Aj_,cos(2n ko Afu + 7/4);0SusAT (19a) 
and 


Xko,q(U) = Aj,,cos(2m kg Afu + 3m/4);0su SAT. (19b) 


where the indices 1 and q denote in-phase and quadrature phase, respectively. 


Thus, hj(u) and h,(u) are described as follows: 


h;(u) = 2 cos(2m kg Af (AT - u) + 7/4); OS us AT (20a) 


and 


hg(u) = 2 cos(2n kg Af (AT - u) + 32/4) ; 0 su SAT. (20b) 
Now the received signal from equation (14), given 7/4 is the phase of the tone 
ko, is 


Yiko(¥) = FikoAlkycos(2n ko Af (u - To ° £(u)) + W/4);uysSusus . (218 


The received signal, given 3n/4 is the phase of the tone ko, is similar to (21) 
with ®), = 37/4 instead of 7/4 as above. 
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As illustrated in Figure 11, the output of the matched filter pair, in gen- 
eral, is the convolution of h(t) and y};,(t) given by 


T 
max 


Zak | = | h. (t) Yi, (u-t) dt 
0 


(22a) 
and 


+ max 
ty (w= [h, (yp, (u-t) dr, 
° (22b) 


where Tax = Ug - Uy. Evaluating (22a) and (22b), the following expressions 


are obtained for the output of the matched filter pairs due to tone ko: 


ay 
Zi (u) = Bik, cos [2k Aftu -£(u-t)- a AT)] dt 
: (23a) 
and 
max 
Zk. (u) = Bik, sin [2nk, Aftu - £(u - Tt) - Ne AT)] dt 
: (23b) 


3] 


yu) 
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Figure 11. Matched Filter Pair for Tone ko 


Inserting equation (9) for £(u) in the above two equations yields 


max 

a 
Aes = MIke lies | cos [(2nk, Af) (u-t, + OT - (1+0,) a AT)) (1+0,) eke 
(24a) 


and 


max 
“eg me Fi Mk, | sin ((2nk, Af) (u-t, + opt - (1+0,) + AT) (1+0,) “dz 
: (24b) 


Recall that the received signal arrives when u, S$ u S$ ug with 


u, =((1 + a) TF] + %] and ug = [(1+ a) AT + (1+ Q)) T, + TH) 
Also recall from the previous section that a, is the Doppler compression/ 
expansion factor due to the moving transmitter. Note that for a) < 0, Tmax = 
Ug - u; S AT and for a) > 0, Tray = Ug - Uy > AT. Using T,,,, = Al / (1 oe 
the upper limit in the integration in the convolution integrals (24a) and (24b), 
will not affect the integrals significantly because a is of the order 10-3; 


therefore, dividing by (1+ a)) is very close to dividing by 1. 
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Sampling Zko iu) and Zko,qiU) at u = ug and evaluating their integrals 


in (24a) and (24b) from 0 to Tyay yields 


| sin [2rk,o] 

=F, A, AT ————  =Z, . 
Alea? usu, Ik, Ik, Eh) dak a Zi i % 
(25a) 


and 
| (1-cos [2xk,a,]) 
“Iko,q MY! uss = Tk ky 8* O*%) “r= ky a (ay) 
(25b) 
The total power, P},, in the kth tone of the 1th baud is one-half the ampli- 


tude, F), Aj, of the received signal squared. That is 
FA = (2 Py,)9-5. (26) 


Also assume 1+) = 1; therefore, AT (1+a)) 1s approximately equal to AT. With 


these substitutions, equations (25a) and (25b) become 


(sin [(2zk.a ]) 
7 0.5 Zl 
(27a) 
and 
(1-cos [2xk.q, ]) 
: 0.5 2% 
Aik. (a) SET, AT a 
(27b) 


A graph of Zj,, (0) and Z),, (04) from (27a) and (27b), respectively, is shown 


in Figure 12. 
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Figure 12. Matched Filter Outputs Versus Doppler Shift 


Note that at oy = 1/4 kg), | Zo 3! = | Ziq5,q!; this is illustrated in Figure 


12. This means that the response of the in-phase and quadrature channels 


are equal even though an in-phase 1/4 symbol was sent; therefore, a signifi- 
cant decoding error would be made if a = 1/(4 kg). It is obvious that, at q) = 0, 


Z)k5,i(0) is maximum and Z),,,4(0) is zero; in which case, no decoding errors 
would be made. The ratio of 2), ; to Zj,5.q below, evaluated for various a's 


less than 1/(4 kg), produces the maximum width of the Doppler channels, Aa. 


34 


The ratio is described by: 


mils J sin(2nk,0,) 


Aik.4 1- cos(2nk, a) 





(28) 


For q = 1/(16 kg), the ratio above is 2 5.03 which is 14 dB. This is called inter- 
symbol separation. 
The value of 14 dB was considered the minimum amount of inter-symbol 


separation to be tolerated. Recall that the maximum amount of Doppler 


mismatch is @)-Q@,,=tAa/2. Therefore, Aa, the Doppler channel spacing, 
is equal to 1/8 kg) and the magnitude of the maximum amount of 


residual Doppler mismatch is Aa/2 = 1/ (16 kg), where Gm 1s the Doppler 


factor associated with the mth Doppler channel. In the simulation, 


Aa = DELALF(LL) = 1(8*KMAX(LL)). 
Solving a - O,, = + Ao/2 for ap, yields 
CO aeeNOy 2 (29) 


where @,,, the Doppler factor due to the mth Doppler channel, is equal to 
m Aa. Substituting a,,=m Aa into (29) yields the following expressions for the 


lower bound, mj, and upper bound, mg, of m, the Doppler channel number as 


my = (a) / Aq) - 0.5 (30a) 


and 


My = (a), /Aa)+0.5. (30b) 
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Since m is an integer between m, and mg, it is computed as 
m= int[(m, + M9)/2]. (31) 


In the simulation, m = M(LL) and is computed using (31); and 


O,, = ALPHAM(LL) = M(LL) * DELALF(LL). 


D. THE ATTENUATION FACTOR DUE TO AN ACOUSTIC 
CHANNEL 


When the transmitted signal, x},(u - 7;,), travels through the medium or 
channel, its amplitude will be attenuated by a factor F), due to the channel. 
In the simulation, the user has the input option to have this attenuation 
applied to the received signal or to let all the received signal tones have equal 
amplitudes. The quantity AA(LL,K) is the amplitude of the 1th baud and kth 
tone of the simulation's received signal. If the user chooses to have equal 
amplitudes, then AA(LL,K) is set equal to 1 for all LL and all K. If the user 
desires the attenuation to be put on the received signal, then AA(LL,K) = Fy, 
which is computed below. 

In general, the transmission loss depends on frequency, depth, pressure, 
and temperature. However, when specific propagation conditions are of no 
interest and only a rough approximation of the transmission loss is adequate; 
the spherical-spreading law plus an added loss due to absorption is equal to 


the transmission loss [Ref. 4]. Therefore, the transmission loss in dB is 


TL = (20 log;) R)+ BR (32) 
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where 
R is the depth of the receiver from the surface in feet and 


B is the attenuation coefficient due to absorption in dB per foot. 
Thus, the following expression for F}, is used in the simulation: 
Fy, = 10 -(TL/20), (33) 


The attenuation coefficient, B, is described by [Ref. 4] 


o1f | 40f 


z + (0.000275) f + 0.003} (3000) 2 
1+f 4.1004 


(34) 


where 
f = (k Af) / 1000 in kilohertz. 


This expression applies for a temperature of 39°F (4°C), which is an average 
seawater temperature, and a depth of about 3000 feet. The constant 0.003 is 
added to take care of the attenuation at very low frequencies [Ref. 4]. In the 
simulation, ABSORP = B. 

The absorption of seawater decreases by about 2 percent for every in- 
crease of 1000 feet in depth [Ref. 4]. Note in the simulation, the transmitter is 
at a depth of 1000 feet below the ocean's surface. This depth is set with the 
simulation variable TXDEP. Thus, the receiver is R = TXDEP + Z(LL) below 
the ocean's surface. TXDEP is not an input of the simulation; therefore, to 
change TXDEP from 1000, one must change TXDEP in the simulation code. 
The simulation computes f$ using equation (34), then increases f by 2 percent 


for every 1000 feet that the receiver's depth exceeds 3000 feet from the 
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ocean's surface or decreases B by 2 percent for every 1000 feet that the 


receiver's depth is less than 3000 feet from the ocean's surface. 


E. THE SAMPLED RECEIVED SIGNAL WITHOUT NOISE 


The received signal, yj(u), is sampled at an interval 
Au = [AT Glee ce) kee (35) 


(Recall wy = 2a k, /[(1 + a,,) AT], (18), is the sampling frequency.) The 
parameters AT and ky are defined in section A. The Doppler factor, a,,=m Aa, 


associated with the mth Doppler channel was derived in the previous section. 


The first sample will be taken at G,, which is the best estimate of uj, the be- 


ginning or synchronization point of the 1th baud at the receiver. This is illus- 


trated in Figure 13. The synchronization error of the 1th baud is Au,=uj-0). 
The simulation variable names for the above parameters are u,=U1, 
0,=UHATI, and Au,=DELUI1. 
The sampled received signal for the lth baud is (13) with u = nAu + 4. 
When 0, = u), the sampled received signal for the 1th baud is (13), 
k, 


y(n) = yj (nAu+u,) = 2Xy;,(nduen,) 021s (ky - 1) 
(36) 


where 


y},(nAu + uz) = Fy, Aj, cos[2mk Af (n Au - (£(n Au + u,) - £(u,))) + ®y, J. (87) 
Substituting (9) for £(u) and computing £(n Au + uj) - £(uj) yields 


£(n Au + uj) - £(u,) = (aq, n Au)/ (1 + Q). (38) 
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Figure 13. The Sampled Receiver Output 


Now substituting Au = [AT (1 + a,,)] /k, into (38) produces 
£(n Au + uj) - £(u,) = [a, n AT (1 + a) /k, (1 + q). (39) 


Collecting the results from (39) and nAu = [n AT (1 + a,,)] / k, together and 


substituting them into (37) yields the following sampled received signal for 


the 1th baud and kth tone: 
y}K(n) = Fy, Ay, cosl(2ak/k,) ((1+0,,) / (14+a))) n + D;,]. (40) 


If a factor v), is added to n to account for random timing jitters of the 1t) baud 


and kth tone due to synchronization error, then (40) becomes 


Y]K(n) = Fy Ai cos|(2nk/k, ) ((1+ay) / (1+a))) (n - Vik) + 0)1] (41) 
fOr Urseniess (i= ly), 
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where 


Au, k 
—————e 
Ik J AiiG@ecms 
2 (42) 
Finally summing (41) over all the tones in the 1th baud gives 
ko 
y(n) = = Yn) UES ts (col) 
k=k 
(43) 


Summarizing, (41) is the final computational form of the received signal 
which is denoted by YY in the simulation. The received signal, YY, is then 
summed over all the k's for each baud from KMIN(LL) to KMAX(LL) to yield 
YRX in the simulation, which is y)(n) above. In the simulation, all the packets 
are contiguous. The continuous sampled received signal is YYRX(II) where II 
ranges from 0, which is the first sample of the first baud (LL=1) in the first 
packet, to NPTS, which is the last sample of the last baud (LL = BDTOTL = 
the total number of bauds) in the last packet. YYRX(II) is equal to YRX plus 
NOISE with all the packets contiguous and the bauds within each packet 
contiguous, where NOISE is white Gaussian noise added to the received 
signal. NOISE, which depends on the input signal-to-noise ratio (SNR), is 


described in the next section. 


F. SIGNAL PLUS ADDITIVE NOISE 
Inevitably noise will be added to the signal, either from the environment, 


the electronics, or both. This noise is assumed to be additive white Gaussian 


noise. Since it is assumed that the received signal, y)(u), has been ideally 
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bandlimited to one-half the sampling frequency, there is no power in 


frequencies greater than or equal to f,/2. Recall that the signal is sampled at 
the transmitter at f,. If w(u) is white noise and has the white power spectral 


density equal to N,/2, then 


E[ wu) ]=E[w(n)]=0 (44) 


and 
o2 = VAR[ w(u) ] = VAR[ w(n) ] = N, f,/2, (45) 


where w(n) is the white noise sequence [Ref. 1). 
Let the receiver input SNR be defined as the signal power in bandwidth 
W divided by the noise power in bandwidth W. The average tone signal power 


in the 1th baud is defined as 
k 


2 
PAVG, = (V/K) »P,, 
k=k 
(46) 
where P), is the total power in the kt tone of the lt baud. (Recall that K is 
the total number of MFM tones in the 1th baud.) Therefore, the wideband 
input SNR for the 1th baud [Ref. 1] is 


SNRWB, = (K PAVG))/(W N,) = PAVG)/(Af N,) = PAVG) k, / (2 6]”), (47) 


because Af = f, /k, and N, = (2 0]*)/ fy. The narrowband input SNR for the 


Ith baud is 
SNRNB, = P};, / (Af N,) = Py, ky / (2 61”). (48) 


Note that if all the tones have equal amplitudes (i.e., AA(LL,K) = 1), then 
PAVG, is equal to P},; in which case, SNRWB, = SNRNB). 
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Solving (47) for 0/2 yields 


ky 
o?=((VK) ), PJk,/(2SNRWB,) . 
(49) 


Recall P}, = (AA},2)/ 2, where AA}, is the amplitude of the kth tone and 1th 
baud of the received signal. In the simulation, NOISE is added to the received 
signal, y)(n) = YRX, which is white Gaussian noise with zero mean and vari- 
ance = 6)* computed in (49). A single input wideband SNR is input by the 
simulation user for all bauds, |. The following simulation variables are used 
for the above parameters: 

6)* = NOSVAR(LL), 

K= Kets dei): 

AA), = AA(LL,K), 

k, = KX(LL), and 

SNRWB, = SNRIN. 
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IV. TESTING THE SIMULATION AGAINST MFQPSK THEORY 


A. BACKGROUND 

The purpose of developing this simulation was to create an experimental 
tool to test and analyze various Doppler, synchronization, and coding tech- 
niques and/or algorithms to support the implementation and testing of an 
MFQPSK acoustic link. Before the simulation could be used in this capacity, 
it had to be tested against MFQPSK signal theory. There was no actual test 
data to use for comparison; therefore, tests were performed to verify that the 
simulation output, the received signal, agrees with the theory of MFQPSK 
signals. It is anticipated that the simulation will be tested further when 


actual test data is available. 


It has been shown that, in theory, the output SNR (SNRoy-) equals the 
input narrowband SNR (SNRNByw) for MFQPSK signals in additive white 
Gaussian noise memoryless channels that are demodulated coherently with a 
Discrete Fourier Transform (DFT) [Ref. 1]. The results of testing the simula- 


tion's ability to reproduce this result are described and presented below. 


B. THE TESTING METHODOLOGY 
The following approach was taken to verify that the simulation performs 


at an acceptable level (1.e., SNRoyp = SNRNByn). A set of 5 signal packets 


were generated with the simulation with a given input wideband SNR 


(SNRWByn). Each packet contained a single unique baud type; therefore, 


each simulation run consisted of the 5 different baud types. Referring to 


Table I, the 5 different baud types are 
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Baud Type 1: 256 sample points, 
Baud Type 2: 512 sample points, 
Baud Type 3: 1024 sample points, 
Baud Type 4: 2048 sample points, and 


Baud Type 5: 4096 sample points. 


Recall from Chapter III that if all the tones have the same power within a 
baud (i.e., the amplitudes of tones within a baud are all equal), then the input 
wideband SNR is equal to the input narrowband SNR for that baud. All the 
runs used for this analysis were produced with normalized amplitudes within 
each baud. The phases that were transmitted were generated randomly by 
the computer code. 

It is known that the minimum input SNR required to accurately decode 
MFQPSK is approximately 15 dB. This analysis considered input narrowband 
SNRs from 0 to 20 dB. The initial conditions for the transmitter's position 
and velocity were chosen from the NOSC track data. These initial conditions 
are 

X0 = -5.0 feet, VXAVG = 1.0 ft/sec, VXVAR = 0.0025 (ft/sec)?, 

YO = 5.0 feet, VYAVG = 5.4 ft/sec, VYVAR = 0.0025 (ft/sec)?, 

ZO = -5400.0 feet, VZAVG = 0.0 ft/sec, and VZVAR = 0.0025 (ft/sec)?. 


The initial conditions for the speed of sound was C0=4900 ft/sec and CVAR=0. 
Note the transmitter was placed almost directly above the receiver to reduce 
the amount of Doppler shift in the received signal. A description and 


examples of the simulation inputs and outputs are presented in Appendix B. 
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An average output SNR was estimated for each baud using the real and 
imaginary parts of the received signal's DFT. A k,-point DFT was performed 
on each baud of the received signal to decode the phase of the received signal. 


(The real and imaginary parts of the received signal's DFT contain the phase 


information sent on each of the tones between harmonic numbers k, and ko 
in each baud.) If ©), = 2/4 was transmitted (which is in the first quadrant), 
then the received signal DFT's real and imaginary parts should be located in 
the first quadrant (i.e., the real and imaginary parts should be positive). 
Likewise, for a transmitted phase angle of 37/4, -3n/4, or -x/4 the received 
signal's DFT phase angle (i.e., real and imaginary parts) should lie in the 
second, third, or fourth quadrant, respectively. The conditional distributions 
within each baud were produced by categorizing the DFT's real and imagin- 
ary parts on a tone-by-tone basis given the phase that was transmitted on 
each tone between harmonic numbers k, and ko. There are two distributions 
of the received signal DFT in each of the four quadrants, a distribution of the 
real part and one of the imaginary part. The real and imaginary parts are 
statistically independent [Ref. 1]. 

For example, in baud type 5 there are 4096 samples of the 256 tones en- 


coded with phase information between harmonic numbers k,=1073 and 
ko=1328. There should be approximately 64 (i.e., one fourth of 256) tones 
transmitted and received in each quadrant because the four transmitted 
phases were generated randomly from a uniform distribution. The distribu- 
tion of each DFT component should be close to a Gaussian distribution. 

Once the DFT of K tones in each baud were conditionally partitioned by 


knowing the transmitted phase (i.e., quadrant) into the four quadrants, then 
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the sample mean and variance were computed on each of 4 real part distribu- 


tions and each of the 4 imaginary part distributions. The sample means are, 


R; and Ie, where the subscript 1 denotes the quadrant number. The sample 


means for each quadrant in a baud are 


N. 
] 
R, = (UN) » R, fori=1,2,3,or4 
Jel (50) 
and 
N; 
5 = (1/N ) Zz I. fori=1, 2,3, or4 
Jel (51) 
where 


N; = the number of times the phase corresponding to quadrant 1 was 
sent in that baud 


Rij = the jth real part of the received signal's DFT corresponding to a 
harmonic number that was transmitted with phase in quadrant 1 


lij= the jth imaginary part of the received signal's DFT correspond- 


ing to a harmonic number that was transmitted with phase in 
quadrant i. 


The sample variances of a given quadrant in a baud are 
_ 
Z =< 
S* = [OND] DR, - Rp 
jas (52) 
and 
Ni 
2 gi 
Si, =(vN,-D] D0, - 1)’. 
i (53) 
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The output SNRs for the ith quadrant within a baud were estimated by the 
sample mean squared divided by the sample variance as follows [Ref.1] 
“~~ = 
SNROUT, , = (R)*/S%. 
(54) 


and 
a (55) 


The mean squared value of the real and imaginary parts of the DFT is a mea- 
sure of the signal power received in that quadrant. The variance of the real 
and imaginary parts is a measure of the noise received in that quadrant. 
These values form an estimate of the output SNR. Appendix C contains the 
statistics for the real and imaginary distributions of each of the four quad- 
rants and the overall statistics for each baud. 

To derive an estimate of the output SNR for each baud, the real and 
imaginary distributions of each quadrant were averaged over all four quad- 
rants. There were 2K points in the overall mean and in the overall variance 
due to the K real parts and K imaginary parts. SNROUT is the overall mean 
squared divided by the overall variance. Figure 14 is a plot of SNROUT ver- 
sus the input narrowband SNR, where SNROUT is described by 














oN 
SNROUT = 


(56) 
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Figure 14. SNRo, 7 vs. SNR, for Five Different Bauds 


where 


N.=K 


4 
=] 


_— 


is the total number of MFM tones in each baud (see Chapter III, Section A), 
Note that the absolute value of the sample means must be used in the 
estimate (56), because, if they are not used, the overall mean will always be 
close to zero. This is due to the fact that the eight sample means are 
approximately equal in magnitude, but four have opposite signs due to their 


quadrant location. 
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C. RESULTS 

If the simulation results agree with the MFQPSK theory, then the input 
narrowband SNR should equal the output SNR. As can be seen in Figure 14, 
baud types 2, 3, and 4 agree within one dB with the MFQPSK theory. The 
longest baud, baud type 5, and the shortest baud, baud type 1, produce esti- 
mated output SNRs that are within two dB of the input SNR. As the input 
SNR increases, it can be observed from Figure 14 that the output SNRs for 
baud types 1 and 5 converge towards the value of the input SNR. 

In conclusion, the output SNR is approximately equal to the input nar- 
rowband SNR; therefore, the simulation reproduces what the MFQPSK the- 


ory predicts. 
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V. OUTPUT SNR DEGRADATION DUE TO DOPPLER 


A. BACKGROUND 

When the MFQPSK signal is transmitted from a moving platform, 
Doppler will shift the received signal's frequencies. The received signal (at 
frequency @,', (17),) is sampled at Wy, (18), baud by baud. In general, q), the 
time Doppler compression/expansion factor, does not equal @,,, the Doppler 
factor set for the mth Doppler channel; hence, there will be some Doppler mis- 
match. This mismatch will cause degradation in the output SNR. 

The output SNR degradation versus the Doppler mismatch was simu- 
lated and compared to theory. If inter-symbol interference (ISI), which is 
energy leaking from one tone to the next tone, is negligible, then the output 
SNR should degrade due to the Doppler mismatch as sinc(ze/4) from (27a) 


where € = @)/Aa. Derivation of the theoretical degradation and the comparison 


to the simulation results follows. 


B. ANALYSIS APPROACH 

The simulation was run several times. Each run generated one signal 
packet with an input narrowband SNR = 15 dB. (The values of 15 dB is con- 
sidered to be the minimum input SNR required to decode the MFQPSK 
signal.) A packet consisted of five bauds of baud type 3 (i.e., k, = 1024 points 
and K = 64 tones). Baud type 3 was chosen based on the results of the 
simulation versus MF QPSK theory in the previous chapter. Recall that baud 
type 3 was in close agreement with the theory of MFQPSK signals (i.e., 
SNRout = SNRNBy). The Doppler channel, m, was set to zero for all the 


o0 


packets. If m = 0, then q - a, = a) - m Aa = q). Thus, by increasing q), the 
amount of Doppler mismatch increases. 

Each packet was generated with a different fixed value of a) =e Aa, 
where € was fixed at 0, 0.25, 0.5, 0.75, 1.0, ..., 2.5, 2.75 and Aa = 1/(8k9). There 
were 12 packets generated, one for each value of e. The phases in the pass- 
band of each baud were selected randomly in the simulation. The initial con- 
ditions on the transmitter position and velocity and the speed of sound were 
the same as used in the earlier simulation analysis of the previous chapter. 

For each baud, the conditional statistics based on knowing the transmit- 
ted phase (i.e., quadrant) and an estimate of the output SNR (56) were com- 
puted with the method described in the previous chapter. A packet's 
estimated output SNR, SNROUT(), is simply the average of the estimated 
output SNR's from each of the five bauds within that packet. SNROUT() 1s 


A S in 
SNROUT(e) = (1/5) », SNROUTC). 


i=] (57) 


Theoretically, if ISI is neglected, the output SNR, SNRoyr should de- 
grade due to Doppler, a, from (26a) as (SNRyy sinc#(2n ko a))). Substituting 
a=e Aa =e/(8 ko) yields 


SNRoyrt = SNRyy + 20 logig [sin(me/4)/(ne/4)] (58) 


where SNRoy Tis in GB. For this analysis, the input SNR, SNRyw, was equal 
a 
to 15 dB. SNROUT(e) and SNRoyrt were plotted versus € to compare the 


simulation results with the theoretical degradation for an input SNR of 15 dB 
(see Figure 15). 
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Figure 15. SNRout vs. ALPHA for SNR yw = 15 dB 


C. RESULTS AND CONCLUSIONS 

There 1s a significant difference between the simulation results and 
theory as shown in Figure 15. It was speculated that this difference is due to 
ISI. As a; increases, the tones "smear" into adjacent tones (i.e., ISI). To prove 
whether or not the difference was due to ISI, the analysis was repeated with 
an input SNR of 5 db. By lowering the input SNR to 5 dB, the signal level is 
not much higher than the noise level within the band; therefore, symbol inter- 


ference from adjacent tones is at about the same level as the noise level of 


that tone. The noise level of each tone, at SNR wn = 5 dB, did not increase as 
much relative to increases in a] as when SNRyw = 15 dB; thus, the output 
SNR at a given tone should not degrade as rapidly as qj increases. Figure 16 
is a plot of the theoretical degradation and the simulation results for both 
input SNR’'s (5 and 15 dB). By examining Figure 16 one can see that, indeed, 
the output SNR at 5 dB degrades at a much slower rate than the output SNR 
at 15 dB. The statistics of each packet for both 5 and 15 dB input SNR ona 
baud by baud basis are presented in Appendix D. 

The analysis of the simulation results for 15 and 5 dB seemed to sub- 
stantiate that ISI can not be neglected. An estimate and correction of Doppler 
shift may be used to reduce the ISI and lessen the output SNR degradation at 
higher SNRyy’'s. An estimate of the Doppler shift can be fed back into the re- 


ceiver sampling system to compensate for the Doppler shift. 
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Figure 16. SNRour vs. ALPHA for SNR, = 15 dB and 5 dB 
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VI. A DOPPLER ESTIMATE FOR FINE CORRECTIONS 


A. BACKGROUND 

The Doppler shifts due to the moving transmitter should be estimated so 
that they can be removed from the received signal in the receiver system to 
reduce the errors in decoding the MFQPSK signal. One technique for estimat- 
ing the Doppler compression/expansion factor was examined. An adaptation 
of a Doppler estimation method recently proposed for ordinary QPSK [Ref. 5] 
was developed for use with MF QPSK and is the subject of this chapter. 


B. IMPLEMENTATION OF THE DOPPLER ESTIMATION 

This Doppler estimation method approximates the Doppler time 
compression/expansion factor within a single baud. To implement this esti- 
mator, several simulation runs were produced. A simulation run consisted of 
a single packet. In each packet, five bauds of baud type 3 (i.e., k, = 1024 
points and K = 64 tones between harmonics 269 and 332) were generated 
with randomly selected phases. The initial conditions of the transmitter's 
position and velocity are the same as used in the simulation versus MFQPSK 


theory in Chapter IV. Each packet was generated with m, the Doppler 


channel number, equal to zero (i.e., a, = 0, which means the residual Doppler 
mismatch, © | - Q, is equal to a) and a fixed a = a, = € Aa where € is between 
0 and 1.5 and Aa = 1/(8kg). For this analysis (using baud type 3), ko is 332; 


therefore, Aa = 1/2656. Recall « = a) is the Doppler compression/expansion 


factor. A set of runs or packets were produced for four different input 


OO 


narrowband SNR’'s (10, 15, 20, and 40 dB). The amplitudes for all the packets 
were normalized (i.e., AA(LL,K) = 1 for all LL and K). 

Comparison of the phase for a given tone in the first half of a baud to the 
phase of the same tone in the second half of the baud yields an estimate of the 
Doppler compression/expansion factor within that baud. If there is no Doppler 
shift and no noise on the signal, the corresponding phases from the first and 
last half of the baud should be equal. The following method was used to 
compute Oy which is an estimate of the Doppler compression/expansion factor 
within the 1th baud. To obtain the phases of the first and second half of the 
received signal's baud, a 512 (i.e., k,/2) point DFT was taken on the first 512 
points of the baud and another 512 point DFT was taken on the last 512 
points of the baud. Note that, when the bauds were generated with the simu- 
lation, the odd tones within the passband (harmonics between 269 and 332) 
were zeroed out by setting their amplitudes to zero (all the harmonics outside 
the passband are always zero). If the odd harmonics in the passband are not 
zero, then the output of the DFT’s of the first and last halves of the baud will 
have interference between harmonic numbers because 1/2 of an odd harmonic 
does not produce an integer harmonic, causing an extra half cycle of signal 
between the harmonics, k’, of the 512 point DFT's. With an ideal received 
signal (i.e., no noise and no Doppler), if the odd harmonics in the passband 
are nonzero, the corresponding phases of the harmonics of the first and last 
half of the baud will not be equal due to the interference of the signals 
between each harmonic in the 512 point DFT's. 


A Doppler estimate, Qy, was then computed for each baud by averaging 


the differences between the received signal's phases from the first half and 
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the received signal's phases from the second half of the baud. Thus, the 


Doppler compression/expansion factor estimate for a baud 1s 


@, = (1MK/2)] >, [PHASE of (S2*(k') $1k))] / (2nk’) 
- 


(59) 
where, 


S2*(k') = Conjugate of the DFT of the second half of the baud 
S1(k') = The DFT of the first half of the baud 


k'= A harmonic number of the first and second half DFT's 
(Note k' = k/2, where k are the harmonics of the original 
ky point DFT of the baud.) 


K = The number of MFM tones in the original baud. 


An average Doppler estimate, 0, was calculated for each packet by averaging 


the five Doppler estimates, Oy, within each packet and was plotted versus o 


for each of the four input SNRs (see Figure 17). Appendix E contains the 


estimates and statistics for each packet at each input SNR. 


C. RESULTS AND CONCLUSIONS 

From the analysis for a = a between 0 and 1.5 Aa, & does not estimate o 
accurately for input SNRs less than 15 dB. Recall that an input SNR = 15 db 
is the minimum for successfully decoding this MFQPSK signal. At an input 
SNR of 15 dB, & estimates o well for « between 0 and 0.8 Aa; and for input 
SNRs of 20 dB or greater, a estimates a for a between 0 and Aa. These 
results are shown clearly in Figure 17. 

The estimate @ is not acceptable for large Doppler shifts (i.e., Doppler 
shifts > Aa); therefore, this Doppler estimation/correction method may only be 
used for fine corrections of Doppler within a received signal baud and Doppler 
channel, where the maximum Doppler mismatch is + Aa/2. Note from 


Figure 17 that the estimation is perfect for 0.5 Aa. 
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Figure 17. ALPHA-HAT vs. ALPHA for SNR = 10, 15, 20, and 40 dB 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


A numerical code to simulate a MF QPSK signal received from a moving 
transmitter through a bandpass channel has been developed. It resides 
presently on the NPS IBM mainframe computer. The simulation is written in 
FORTRAN 77 and the code is included in Appendix F. 

The model of the transmitter platform dynamics in the simulation is lim- 
ited to straight line motion with random fluctuations. Actual transmitter 
platforms have motion in the pitch, yaw, and roll planes too and thus a six- 
degree-of-freedom model may be more realistic for future studies of platform 
motion effects on MF QPSK signalling. 

The simulation should be tested further with actual test data to ensure 
that it produces a realistic received signal. Validation with actual test data 
may show that other channel characteristics or parameters which have not 
been considered in the simulation are important and must be included. 

An important discovery of the simulation was that the output SNR 
degradation due to Doppler shift was primarily due to ISI. ISI cannot be 
neglected; therefore a theoretical analysis of ISI should be developed to real- 
istically analyze MF QPSK in the presence of Doppler. 

The Doppler estimation/correction method analyzed and presented in 
chapter VI only appears useful to remove small Doppler shifts (.e., for fine 
adjustments only). Analysis of methods to remove larger Doppler shifts (i.e., 
coarse adjustments) across several bauds or a packet should be performed for 


MF QPSK signals. 


og 


Synchronization error estimation and correction has not yet been 
addressed with this simulation, but should be in the future. 

This NPS developed MF QPSK signal has been implemented in hardware 
by T. Gantenbein and Dr. P. H. Moose using differential phase coding [Ref. 6]. 
It is conjectured that with differential coding the Doppler shifts, ISI, and syn- 
chronization errors will not affect the received signal as significantly as they 
affect the individually phase coded signal which is the coding technique 
presently implemented in the simulation. Therefore,the differential phase 
coding should be coded in the simulation. Doppler, ISI, and synchronization 
estimation/correction methods should be studied that are compatible with dif- 


ferential phase coding. 
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APPENDIX A 
DERIVATION OF £(u) 


Given the model of the channel's time delays below in Figure 18, assume 


x(u) is a superposition of N impulses; 


N 
xu) = >» A. ) (u-u,) 
i=l 
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Figure 18. Time Delay Model of Channel 


An impulse sent at time u = 0 arrives at y at time 7 + T + QT. There- 
fore, the output signal from the channel y(t; + T; + aT) is the impulse x(0); 
that is, y(t}, + T) + @)T,) = x(0). An impulse sent at time AT arrives at the 


output of the channel at time 


AT + T] + T + a(AT + Th). (61) 


Therefore, the output signal from the channel 


WAT + tT] +7) + OAT + TK)) (62) 


6] 


which is the impulse x(AT); is 
yVAT + |] +) + AAT + T)) =x(AT) . 
Now letting 
u=AT+ IT] +7 + AAT + T%}) 
which can be rewritten as 
u=AT(1+ a)+%(1+a)+7, 


then solving for AT yields 


AT =((u-7)/(1+q))-TJ=u-T- (t+ ayu)/ (1 + @)). 


Thus, 
y(u) = xCAT ) = xCu - TR - (T+ aya) / (1 + @))) ) 
or 
y(u) = x(u - Tr - £(u)) 
where 


£(u) = ((Tj + oyu) / (1 + @))). 


62 


(63) 


(64) 


(65) 


(66) 


(67) 


(68) 
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APPENDIX B. AN EXAMPLE OF A SIMULATION RUN 


A. INPUT FROM THE SCREEN 


Beene sehen Hee INT TEAL POSTTION x0,Y0,Z0 (CFT) 
Tepe Reece vEeR RELATIVE STO JHE TRANSMITTER ... 
nu —- 10.0000 (= 20.0000 ee =O 000 000 


EMTER THE TRANSMITTER"S VELOCITY AVERAGE AND VARIANCE 
Cie eeaelheemhOhomorl/SECy) AND CFITZSEC)XX2 
CI.E. VXAVG, VXVAR, VYAVG, VYVAR, VZAVG, VZVAR) 


VXAVG = 2.00000000 VXVAR = 0O.500000007E-0O1 
VYAVG = 4.00000000 VinvoRe Ue GO 9999S 9GE=02 
VZAVG = 0.000000000E+00 VZVAR = 0.29999999G4E-02 


itch lhe eAveERAGE STEED OF SOUND IN FI/SEC 
BAUS tHe VAR EANCE IN CFIZSEC)XX2 .... 
CO = 5000.00000 CVAR = 0.000000000E+00 


Pen ime PRANSMITTER"S DOWH LIWK TRANSMIT WED 
foo Den NDR eANGLE (DEG)  ... 
THETAO = 4.00000000 


Deer OU WAT WORMALIZEDVAMPLITUDE FOR THE RECEIVED SIGNAL? 
PeSMoOEMEN EN ese Yeo OR OsHO 5. 
IAMP = 1 


BRiex THE #eOr PACKETS IN THE TRANSMITTED SIGHAL 
@ 


HPAKS = 
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THIS PROGRAM ENCODES A QPSK MULTIFREQUEN Cy Sarena 
THE PHASES ARE SHOWN BELOW FOR ONE FREQUENCY 


SELECT ONE OF THE FOLLOWING METHODS EGR ENContiit 
THE PHASES FOR EACH FREQUENCY FOR EVERY SBAUD 
HITHIN- EACH PACKET 253 


Etiver 


1: THE PROGRAM RANDOMLY SELECTS ALL THE PHASES 

2 PACKETS 

2: YOU INDIVIDUALLY SEUEC) Sie aE nAsies 
Ih 


METHOD = 


ENTER THE BAUD TYPE # CORRESPONDING 


FOR ALL 


TO THEGEOLEOMNiiG 
| eRe 


BAUD LENGTH FOR PACKET NUMBER: 
] BAUD LENGTH CDELT) = 17240 SECONDS 
Z BAUD LENGTH CDELT) = I7i20SSEConr. 
5 BAUD LENGTH «DELI. =" Ileus SEConDs 
G BAUD LENGTH. CDELT) = a7 Ss0R eeu 
IS) BAUD LENGTH (DELT) = 715" Sc ceonte 
BDIYE EG)? ] 
ENTER THE NUMBER OF BAUDS IN PACKET NUMBER: ] 
NBAUDS(1) ] 


ENTER THE BAUD TYPE # CORRESPONDMNG 


TO0.THE FORLOHIEE 
Ci tee. 


BAUD LENGTH FOR PACKET NUMBER: 

1 BAUD LENGTH (DELT) = 17290 SECONDS 
2 BAUD LENGTH (DELI) = 17120 (SEeenes 
3 BAUD LENGTH CDELT) = 1760 SEGEIiS 
4 BAUD LENGTH (DELT) = 730 SEGUMis 
S BAUD LENGTH CDELT) = i715 SECON: 

DDL TP eGze 7a 

ENTER THE NUMBER OF BAUDS IN PACKET NUMBER: 2 


HBAUDS(C2) 


] 
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WOULD YOU LIKE THE DISCETE FOURIER TRANSFORM 
OF (He sOUTEUiIRS CHAD oR SCERTER JT:YES OR O:NO) 
IDFT = 1 


WOULD YOU LIKE THE DFT OQUTPUT WINDOWED ? 
BNTER = “VES sORmUE NO. . 
IWNDOW = Al 


PLEASE ENTER THE DESIRED INPUT WIDE BAND 
SIGNAE=TOsNCLSE RATIO IN GDB... 
SNRDB = 15.0000 
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B. LIST OUTPUT TO DISK FILE 30 


48 
LL K PHICL ES IPHI 

68 =0.7 65 
69 Oa 7 oD 
-70 OR es toys, 
ca aes 
72 =22956 
cS =239 960 
74 =2 mone 
i =07 3 
76 2856 
ee =O 65 
78 fey) ee) 
9 Ory 35 
80 =U 55 
81 Ua ees, 
82 0 62 
83 SCuooe 
eS5 =O eo 
136 estate, 
13 Ba Sells, 
135 256 
137 =o Sle 
140 FEN, 
141 =02705 
142 =O0h7 55 
14 O27 65 
papas) 

Zo as 

= 02-665 


1 | i | 
NMPNOONGCOOCN 
~ 
Oo 
6 a) 


NMP NNM CONN PP 
CN 
ui 
On 
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510 0O.397997E+01 O0-000000E700 

511 0.242Z562E+01 O.000000ETU0 

DFT OUTPUT DATA FOR BAUD & 2 
K REAL PART IMAG PART MAGNITUDE 
135 0.21869614E+03 -.13437286E+03 0.25667896E+03 
136 -.20494591E+03 -.20220082E+03 0.28790259E+03 
137 -.14575168E+03 -.15294678E+03 0.21127296E+03 
138 -.21255940E+03 0.19105444E+03 0.28580273E+03 
139 -.15577037E+03 -.13739233E+03 0.20770425E+03 
140 -.13250822E403 0.14138666E+03 0.19377464E+03 
141 0.17392702E+03 -.22279741E+#03 0.28264697E+03 
142 0.13396306E+03 -.13225531E+03 0.19183899E+03 
143 0.16901073E403 0.14665094E+03 0.22376578E+03 
144 0.19780035E+03 -.15813335E+03 0.25327560E+03 
145 -.20230388E+03 0.27242383E+03 0.33932495E+03 
146 0.20827090E+03 -.18137456E+103 0.27617651E+03 
147 ~.19098985E+03 0.22660436E+03 0.29620264E+03 
148 O.15571500E+03 0.14765985E+03 0.21459401E+03 
149 0.22490648E+03 0.18650G17E+03 0.29217578E+03 
150 0.18717580E+03 0.16107140E+#03 0.24693380E+03 
151 -.13020734E+03 -.21927161£+03 0.25501762E+03 
152 0.21125522E+03 -.22968059E+03 0.31206079E+03 
153 0.16326271E+#03 -.17587959E+#03 0.23997571E+03 
154 -.18001019E+03 0.15380461E+03 0.23676892E103 
155 -.19258081E+03 -.13529289E+03 0.23535405E+03 
156 -.1902/963E+03 0.16097437E+03 0.24923732E+03 
157 -.21412985E+03 0.18872530E+03 0.28542749E+03 
158 -.20541913E+#03 -.22383812E+03 0.30381006E+03 
159 -.19955086E+03 0.20583072E+03 0.28668237E+03 
160 ~-.21075299E+03 0.22362857£+03 0.30728906E+03 
161 0.17198572E+03 0.17881973E+03 0.24810397E+03 
162 -.17657633£+03 0.15956850EF+03 0.23799471E+03 
163 -.18814537E+03 0.14442966E+03 0.23718896E+03 
164 -.17734966E+03 -.17915202E+03 0.25208797E+03 
165 -.18872591E+03 -.18354022E+03 0.26325732E+03 
166 -.18190575E+03 0.12815350£+03 0.22251521£+03 
The DFT is windowed (CIWNDOW = 1) 
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C. GRAPHICAL OUTPUT USING DISSPLA 
1. DFT Windowed Between ki and kg (IWNDOW=1) 
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2. DFT Not Windowed (IWNDOW=0) 


OF QUTPLT Cle Urle aCe ST istWAL 
FOR BAUD NUMBER 


NAGHI TUDE 
75: < 





0 780 into 12 so) 2D0 750 00 27750 
Saab ae 


eee ee Cee ce |} ae ST SshlAaL 
FOR BAUD NUNBER 


>. ee F- 
_— — — — 
=— —— i U2 EE 8 OOo ee 
— 7: =e Q 
== —— — | — | 
—— 
eal 
eee Gar ee 
ee emer a Toe 
= aa ell 
-—— — 
[ea eee one 3 
Lo a. - 
S a nc 
=—jey= _ - 


PHASE (DEG) 





85 


NACHT TUDE 
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APPENDIX C. STATISTICS OF SNRourt VS. SNRIN ANALYSIS 


A. 


INPUT SNR = 0 dB 
INPUT SNRNB = 1.000 = 0.000 DB 
BAUD TYPE 1: KX = “256 SAMPLE POINTS; K = 16 TONES 


®* GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT x 


NUMBER OF TONES = 4 
MEAN OF DFT REAL PART = Smet O99 5S 
VARIANCE OF DFT REAL PART = Salo 3 / 
QUADRANT SNROUT OF DFT REAL PART = 1.1347 = 0.5486 DB 
MEAN OF DFT IMAG PART = 129216241 
VARIANCE OF DFT IMAG PART = 6640.81 
QUADRANT SNROUT OF DFT IMAG PART = Geolea = 4.0005 DB 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xx 


NUMBER OF TONES = 2 
MEAN OF DFT REAL PART = aie re 0 0 
VARIANCE OF DFT REAL PART = 429.54 
QUADRANT SNROUT OF DFT REAL PART = 16.9056 = 12.2803 DB 
MEAN OF DFT IMAG PART = 82.06499 
VARIANCE OF DFT IMAG PART = 27 or. 50 
QUADRANT SNROUT OF DFT IMAG PART = 2.4158 = seoou0. DB 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xXx 


NUMBER OF TONES = 6 
MEAN OF DFT REAL PART = womeo oo 0 
VARIANCE OF DFT REAL PART = VS a8 1 
QUADRANT SNROUT OF DFT REAL PART = tees gs - -p.9069 DB 
VeeieOr Dili wiiinG PART =" -15/7.14986 
VARIANCE OF DFT IMAG PART = Peon 59 
QUADRANT SNROUT OF DFT IMAG PART = 21.8667 = 13.3978 DB 


XX GIVEN THE TRANSMITTED PHASE IS IN THE STH QUADRANT xXx 


NUMBER OF TONES = G 
MEAN OF DFT REAL PART = MS Ic 4S 
VARIANCE OF DFT REAL PART = Foo «41 
QUADRANT SNROUT OF DFT REAL PART = 2.4520 = Sooo DB 
Hees OF DERPSIMAG PART = -139.99750 
VARIANCE OF DFT IMAG PART = Proc os 1 
QUADRANT SNROUT OF DFT IMAG PART = 5675" = io > 4? DB 


REX OVERALL CREAL + IMAG) STATISTICS FOR BAUD TYPE | Ss: 


TOTAL NUMBER OF PODMIS,. 2k = oye 
BAUD MEAN = Pla a650 56 
BAUD VARIANCE = Hoo .15 
BAUD SNROUT = 2.2/43 = 3.56846 DB 
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INPUT SNRNB = 


BAUD TYPE 2s) ekx 


1.000 


S512 SAMPEE FOTN: 


0.000 DB 


K 


3¢ TONES 


** GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DET eRe 


DFT 


AL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


2 


BG Wrasse | 
1499349.4] 


le. 2653" = 


1.0151 DB 


263:,06552 
54654.25 


ll cya 


1.04648 DB 


xX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


OF 


Diet 
Dean 
DFT 


DFT 
Died 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


“wou ud u 


8 


=12673)4920 
60110 .95 


0.2697 


=—P=5.7 7 Ges 


137 .GelS 
19688..26 


0.9620 


= -0.16846 DB 


**® GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xXx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SHROUT 


OF 
OF 
OF 


OF 
OF 
Om 


DFT 
DFT 
DFT 


DFT 
DFT 
Bie) 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


wou oi a 


10 


=191 87025 
LG ale o 


1.8793 


= 2./ 255.05 


=28O,7 9902 
SUS TZ.35> 


2. fous 


= 4.3956 DB 


** GIVEN THE TRANSMITTED PHASE IS IN THE GTH QUADRANT xx 
MUIMB ER Cie iiss 


MEAN 
VARIANCE 
QUADRANT SHROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


Dien 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


2) 
156 abo Sec 
4049449 .66 


0.592 


= = 2j42756 DE 


=21 52559500 
21217 286 


2.1356 


= 3.596. DE 


X¥X OVERALL CREAL + IMAG) STATISTICS FOR BAUD TYPE 2 XxXx 


TOTAL NUMBER OF POINTS, 
BAUD MEAN 
BAUD VARIANCE 

BAUD SNROUT 


88 


2K 


64 


191 ..450Zze 
Zac) ls 56 


ime (4 


1.1308 DB 


INPUT SNRNB = 


BAUDet yRE St KX 


te 


_ 
-— 


000 


LUZS saree POINTS; K = 


0.000 DB 


64 TONES 


*X GIVEN THE TRANSMITTED PHASE IS IN THE IST QUADRANT Xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


OF 
de 
OF 


OF 
OF 
OF 


DEI 
DET 
DFT 


Da 
DET 
Diet 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


13 


2991S) 
171068.06 


0.5108 


= eee? 6DB 


(fe tn ae kagtes 
(desl. 61 


U. 9679S 


= -0.0496 DB 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT Xx 
HUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DED 
DFT 
DET 


DFT 
DET 
DET 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


15 


BONICat1o0/ 
moro 7 9.00 


2.3405 


= So750 DS 


393.86084% 
CEN a ee ee, 


1.6876 


= Zret2o DB 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xXx 
NUM DER SOF TONES 


MEA 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DET 
Ding 
Dit: 


Dine 
DFT 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


1 oot 


20 


Goo) 7 
locuNo. 0 


0.8958 


= -0.72/74 DB 


meee 1 9G 
TU0S0Gs 31 


0.7447 


Ss eer co Ol: D5 


*X GIVEN THE TRANSMITTED PHASE IS IN THE 4TH QUADRANT xXx 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


NUMBER OF TONES 


Die 
Bint 
DET 


Dial 
Dil 
DET 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


mo ou 


16 


Sco, 0 924 / 
97164.81 


Les 


= 0.4565 DB 


wag. 95506 
125734.50 


eavey 9 


= S07 US ops 


<A Ov eERPee CREAR + IMAG) STATISTICS POR BAUD TYPE 3 xxx 
NOTAL KUMBERSOF POLHTS, 


89 


2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 


ot lat 


12s 


368 .994938 
ors can 9 


| i 2 es 


0.8540 DB 


INPUT SNRNB = 


BAUD TYPE 4: KX 


My 


000 


0.000 DB 
2048 SAMPLE POINTS; 


K 


12s TONGS 


xX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT Xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


F DFT 


DEI 
DFT 


DFT 
DET 
DET 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


38 


165. 55200 
G80832..57 


poe loy 


= 0.8596505 


G87 .95532 
636912. 256 


0.5450 


= -2.6363508 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


OF 
OF 
OF 


OF 
OF 
Ole 


DET 
Dial 
DFT 


DET 
Digit 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


34 


-919.91406 
GPeat dao? 


leno 7 


= 2.5181) DE 


HOS SF (8) Sys) 
GUS56 2. 00 


220999 


= 3.1 1D7 3508 


xX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT Xx 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


NUMBER OF TONES 


Diet 
Dial 
DET 


Dew 
Dri 
DET 


REAL 
REAL 
REAL 


IMAG 
INAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


BS 


-749.60986 
DiS De 6i2 


1.0355 = 


0 S06 “ie 


=6 4 dio6 2.06 
4939 185.0 


0.8446 


= SOS Guerre 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 4TH QUADRANT Xx 


MEAK 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


NUMBER OF TONES 


Bp 
Ble I) 
Diet 


DFT 
DFT 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


2i 


845.749316 
506747 .94 


LAG 5 


= 1.4968 DB 


=O)07e25 9955 
GPea0 Dili 


0.7808 


= -1.0796 DB 


¥XX OVERALL CREAL + IMAG) STATISTICS FOR BAUD TYPE @ xxx 
TOTAL NUMBER OF POINTS, 


90 


2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SHROUT 


256 


I S6Reorae 
460690.81 


Ito 3t 


0.730)g0s 


INPUT SNRHB = 


BAUD TYPE 5: KX 


Le 


000 


4096 SAMPLE POINTS; 


Og 00 D5 


K 


Zoo SRUMES 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 15T QUADRANT Xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DFT 
Dink 
DET 


DFT 
DFT 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


64 


Pos moe? Lo 
2880635.00 
Un7 69S 


1674 .69238 
€977856.00 
1.0880 


=n 0279 DB 


= 025661 DB 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xXx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DFT 
DET 
Die 


Diet 
DET 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


61 


TaGeorre I 226 
2406180.00 
So 


$545 .89106 
32744981 .00 
0.5552 


= 0.9671 DB 


=a -2oo1)> DB 


**% GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xXx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SHROUT 


DET 
DE 
he |p 


DED 
Det 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


Br) 


-14947 .54858 
1710154.00 
NE 4045, 


SL 20 Cr 0915 
26866264%.00 
0.5437 


= 0.8823 DB 


=e ooo DD 


X¥ GIVEN THE TRANSMITTED PHASE IS IN THE 4TH QUADRANT xXx 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


HUMBER OF TONES 


DFT 
DE 1 
DET 


Diet 
le 
Die 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


Ue 


Pola US 5 
9270.00 
0.8680 


ee Osen >) DB 


—lSoleso7 92 
2939624.00 
02 G0S7 


=secU-92c0 DB 


POV ERA SCREAL + IMAG) STATISTICS FOR BAUD TYPE 5 Xxx 
TOTAL NUMBER OF POINTS, 


ome 


2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 


Die 


Pao ure ll 
2524291 .00 


0.38956 


wun coor UE 


B. INPUT SNR = 5 dB 


INPUT SNRNB = 5. V6 2ae >.000 "DE 
BAUD TYPE 1: KX = 256 SAMPLE POINTS) K = eu 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 


NUMBER OF TONES = “i 
MEAN OF DFT REAL PART = 29 sua 
VARIANCE OF DFT REAL PART = 2661.94 
QUADRANT SNROUT OF DFT REAL PART = Sesh = 5. coo) Die 
MEAN OF DFT IMAG PART = 112.22491 
VARIANCE OF DFT IMAG PART = c09emr? 
QUADRANT SNROUT OF DFT IMAG PART = 6.0008 = 7 Tae De 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xx 


NUMBER OF TONES = 2 
MEANMOF DET REALS PART — —O7 250000 
VARIANCE OF DFT REAL PART = lo 0260 
QUADRANT SNROUT OF DFT REAL PART = 56.4194 = 17.5143 DB 
MEAN OF DFT IMAG PART = Gi a0 010 
VARIANCE OF DFT IMAG PART = 883.68 
QUADRANT SNROUT OF DFT IMAG PART = §.3268 = 9.2048 DB 


xX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xx 


NUMBER OF TONES = 6 
MEAN OF DFT REAL PART = =1 Gaus 
VARIANCE OF DFT REAL PART = 4080.46 
QUADRANT SHROUT OF DFT REAL PART = | 35229 = 1. 320S5DE 
MEAN OF DFT IMAG PART = -l2yeee5sc 
VARIANCE OF DFT IMAG PART = 537.9) 


QUADRANT SNROUT OF DFT IMAG PART 45.8118 = 16.6098 DB 


XX GIVEN THE TRANSMITTED PHASE IS IN THE GTH QUADRANT xx 


NUMBER OR SOE oe— G 
MEAN OF DFT REAL PART = 15889991 
VARIANCE OF DFT REAL PART = COO lee 4 
QUADRANT SHROUT OF DFT REAL PART = See PalbN 7 GGnaeedB 
MEAH OF DFT IMAG PART = -1l>e592492 
VARIANCE OF DFT IMAG PART = JO 7-60 
QUADRANT SNROUT OF DFT IMAG PART = 3.6278 = 5. 5976S DE 


RX OVERALL CREAL + IMAG) STATISTICS FOR BAUDST Eee 


TOTAL NUMBER GF POINTS aa ae 
BAUD MEAN = LO) cely 
BAUD VARIANCE = oo 226 
BAUD SHROUT = 5.9656 = 7.7365 05 


Sia 


INPUT SNRNB = 


BAUR LIVPE 2: KX 


me 


162 


5.000 DB 
SLerooiioe FOLITS; 


K 


SZ ONES 


** GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DFT 
Dial 
DET 


DFT 
Dit 
DET 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


Wott ot a 


5 
156254995 
4717.56 


ye Home 


= 7.1498 DB 


227 .40594 
Wego. 6 9 


29 oe 


= ae DD 


** GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT x 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIAHICE 
QUADRANT SHROUT 


Be | 
DFT 
DFT 


DFT 
BAe | 
DFT 


REAL 
REAL 
REAL 


IMAG 
INAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


it i 


8 


-150.08290 
Poi | 


1.1848 


= Oeyes6 5 DB 


ste ne ees 
6250 401 


Se Oe 


= 5.94964 DB 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xx 
HUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DEL 
Di 
Dirt 


Dial 
Diet 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


Wot dl 


10 


el oom oe 5G 
GIES eRs 


Seco 


= 7eos6> DB 


= Gleeno4.9 5 
S aXs) (ral) 


6 Uce 


= 7.84927 DB 


X* GIVEN THE TRANSMITTED PHASE IS IN THE 4TH QUADRANT xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


i= 


De 
Dell 
Diet 


Dial 
Det 
ie I 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


2 


166 .15880 
Poe, 1 


2.1588 = 


Saodac DB 


SAWN 2 rests bl 
G09 D 


Sy eel 


= 7.7642 DB 


eee ERA CREAR OEMAGI STATISTICS FOR BAUD TYPE 2 *xx 


NOTA NUMBER OF POLIS, 
BAUD MEAN 
BAUD VARIANCE 

BAUD SNROUT 


93 


2K 


64 


SiG. 1 ODO 9S 
8934.94 


5 0468 


> 


S.9160 DB 


INPUT SNRNB = 


BAUDT Ti EG Oo Kx 


Se 


Gre 


». 000 Ds 
1024 SAMPLE Patiiis: 


K 


64 TONES 


*¥¥ GIVEN THE TRANSMITTED PHASE IS IN) TRE TS] eQUATE Ti 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SHROUT 


Die 
DEG 
DEG 


Dat 
Dies 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


2.702 


Si gwleo.G 
25/8547 


G.0967 


= 6,124e) DE 


X* GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DFT 
DET 
DFT 


DET 
DET 
Diu 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


1, 


-441.42603 
Sol Zio? 


5.7106 


= 7.5666 .DE 


380.12280 
2906791 


4.9709 


= 6. 969o 8b b 


*X GIVEN THE TRANSMITTED PHASE 1S IN THE 3RD QUADRANT Xx 
NUMBER WOE. TOES 


MEAN 
VARIAHCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SHROUT 


ie Al 
DED 
DET 


DET 
DFT 
DFT 


REAL 
REAL 
REAL 


IMAG 
TtHiAG 
ItiAG 


PART 
PART 
PART 


PART 
PART 
PART 


{I 


20 


369.127 20 
49952.497 


2.6 5a 


= G 2395 .0n 


=S20 moze 12 
S207 6 


2.9304 


= 4.6693 DB 


** GIVEN THE TRANSMITTED PHASE IS IN THE GTH QUADRANT Xx 
HUMBER OF) TOtles 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIARCE 
QUADRANT SNROUT 


OF 


ving 
DFT 
DET 


DF il 
Bile i 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
FART 


PART 
PART 
PART 


16 


S95 7 o S005 
507 35.597 


Re) ese, 


= 5.8429 DB 


-442.496045 
3973494.58 


4.92/70 = 


6. 925Gb eB 


¥¥% OVERALL CREAL + IMAG) STATISTICS FOR BAUD Si ees 
TOTAL, HUMBER OF PUGiiie. 


94 


2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SHROUT 


128 


Es VSS) | 


S501 G0 


Seog eg 


5 ../ 8 SGasub 


INPUT SNRNB = 


BAUD TYPE G: KX 


5 


Gi 
2048 SAMPLE POINTS; 


5) 5 NOON DNs, 


K = 128 TONES 


**% GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SHROUT 


DFT 
DFT 
DET 


DFT 
DEI 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


38 


7549.154930 
1513560 0ow 


Sf Deo 


= pe 49 DB 


59677 Ul/ 
ioocos.00 


2.5929 = 


Geilo?) DB 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xXx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


Die 
DFT 


DFT 
DET 
Diet 


fF DET inet 


REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


wot te 


34 


-842.05663 
150349.12 


4.7160 


8 
1 
aie 


D~ 


= oe 57 DB 


59./70776 
2 


0o.20 


= 70091 (DB 


¥X GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xx 
NUMBER QF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DFT 
DFT 
DFT 


DFT 
Dink 
Det 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


SS, 


~741.81006 
166615731 


2ecoey 


= Seo. DB 


=666.025)1./ 
Ips is Fake Bar aas, 


om sec 


= 4.8267 DB 


XX GIVEN THE TRANSMITTED PHASE IS IN THE STH QUADRANT xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


OF 
Giz 
OF 


OF 
Or 
Oke 


DFT 
DFT 
DET 


DFT 
DFT 
DET 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


21 


1 Soran 0:21 
1596 o0n 


G.0027 


= jrviazoo DB 


SOO Sal ei.) 
1969754206 


2.95464 


= 4.7047 DB 


OVER AeREAL + IMAGESTATTsnicS FOR BAUD TYPE 4 xxx 


HOTAt NUMEER OF POINTS, 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 


95 


2K 


256 


EH aS ys ie 


MF. AST OG 
MG a7 9.56 


5.74565 DB 


INPUT SNRNB = 


BAUDMI Reo KX 


SF 


lG? 


4096 SAMPLE POLNTs:; 


5.000 86 


K 


256 1 ONEe 


xx GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DET 
DFT 
Diet 


DEW 
DFT 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


64 


1508221269 
915649.00 
2.98493 


16053455105 
825299 225 
3.1304 = 


= 5. ISA0R DE 


4.9560 DB 


** GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xx 
NUMBER OF TONES 


MEAK 
VARIANCE 
QUADRANT SNROUT 


NEAN 
VARIANCE 
QUADRANT SHNROUT 


EoD ak 


Dial 
DFT 


DET 
DET 
DET 


EAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


61 
=1585256665 


702008. 57 

4) Ss 5. 1825905 
=. 1915377759 
= 1035020256 
= ] . 9SGioe = 2.87049 DB 


** GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT Xx 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SHROUT 


NUMBER OF TONES 


OF 


DEF 
DFT 
DFT 


DET 
Digi 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


Do 


=1497 3573535 
55e 94576 | 
G.0282 


= 6. US1250E 


=1549 259515517 
845222001 
2 elo 


= 3/2365" DE 


** GIVEN THE TRANSMITTED PHASE IS IN THE 4TH QUADRANT xXx 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


NUMBER OF TONES 


Digi 
DE 
DEY 


DE 
Be | 
DEi 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


(2 


1401.73390 
652/765). 94 
Bp Ue 2 


= 4.9209 DB 


=1 5265.5 50510 
9279 io? 
2. OMe 


= 3° 99S SDE 


X**X OVERALL CREAL + IMAG) STATISTICS FOR BAUD TYPE 5 XXX 


TOTAL, NUMBER OF SR Onis: 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 


96 


2K 


Spe 


1483.02075 > 


TI a e69 


rapier og he 


4.3956 DB 


C. INPUT SNR = 10 dB 


INPUT SNRNB = 10.000 = 10.000 DB 
Pauoetthewtle WX =s coo SAMPLE POLTHTS; K = 16 TONES 


*X GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 


NUMBER OF TONES = G 
MEAN OF DFT REAL PART = ras Troe ak 
VARIANCE OF DFT REAL PART = 840.84 
QUADRANT SNROUT OF DFT REAL PART = 10.2386 = 10.1024 DB 
MEAN OF DFT IMAG PART] = 10227 1296 
VARIANCE OF DFT IMAG PART = 665°. 56 


QUADRANT SNROUT OF DFT IMAG PART poe oeeoe—  lem0ls/ DB 


X* GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xXx 


NUMBER OF TONES = C 
MEAN OF DFT REAL PART = —SO7e 9999 
VARIANCE OF DFT REAL PART = ao? 
QUADRANT SNROUT OF DFT REAL PART = 183.6700 = 22.6404 DB 
MEAN OF DFT IMAG PART = 87.849000 
VARIANCE OF DFT IMAG PART = eo. 95 


QUADRANT SHROUT OF DFT IMAG PART 27.6602 = 14.9185 DB 


xX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xx 


HUMBER BOE TOMES: = 6 
MEAN OF DFT REAL PART = =e ood 90 
VARIANCE OF DFT REAL PART = 26 9.90 
QUADRANT SNROUT OF DFT REAL PART = 5.1394 = 7.10499 DB 
Hien Of DET TiAG PART Fino 66 


VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 


iS 1 9 
110.0076 = 20.4142 DB 


XX GIVEN THE TRANSMITTED PHASE IS IN THE GTH QUADRANT xXx 


HUMBER OF TONES = 4 
MEAN OF DFT REAL PART = 104.76740 
VARIANCE OF DFT REAL PART = too. <2] 
QUADRANT SNROUT OF DFT REAL PART = 14.5137 = 11.6178 DB 
MEAN OF DFT IMAG PART = -104.5949] 
VARIANCE OF DFT IMAG PART = 1164.08 
QUADRANT SNROUT OF DFT IMAG PART = ee kod la Ore 0S: DS 


XXX OVERALL CREAL + IMAG) STATISTICS FOR BAUD TYPE 1 Xxx 


TOTAL NUMBER OF POINTS, 2K = ae : 
BAUD MEAN = 9762696 
BAUD VARIANCE = 568.06 
BAUD SNROUT = 16.8477 = 12.2654 DB 


aie 


INPUT SNRNB = 1020000 = SiGe UC Uw rs 


BAUD TYPE 2: KX = 512 SAMPLE POINTS? f= eee oie 


*¥*X GIVEN THE TRANSMITTED PHASE IS IW THE 1ST QUADRANT xx 


NUMBER OF TONES 

MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 
QUADRANT SHNROUT OF DFT REAL PART 


MEAN OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 


5 
167 05995 
1490.65 
18.7182 = 12.7226808 


207 200189 
57a 
( Reepeliy 8.9386 DB 


tn) el | 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xXx 


NUMBER OF TONES 

MEAH OF DET READ Rei 

VARIANCE OF DFT REAL PART 
QUADRANT SNROUT OF DFT REAL PART 


MEAN OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 


8 
=1Go 25 249.5 


4.4490 = 6.4826 DB 


ot mt ot 
On 
So 
— 
fo) 
IN 
NO 


Wek s 
ll.5 2s 


*X GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT XxX 


HUMBER OPS uOies 

MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 
QUADRANT SNROUT OF DFT REAL PART 


MEAN QF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SHROUT OF DFT IMAG PART 


10 
=164225990 
LG Ses 
17.2714 5. 12.37 30m 


=219.98996 
5 0545356 
15.2519 = 11.827 4505 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 4TH QUADRANT xXx 


NOMBER OF TOMES 

MEAN OF DFT REAL Paw 

VARIANCE OF DFT REAL PART 
QUADRANT SNROUT OF DFT REAL PART 


MEAN OF. DET IMAGRE ART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 


2 
L256 655 
4046.22 
7.261 5e- 8.6696 DB 


=) 9sG9 959 
2 Vee 6 
l7.32Noe-) l2é.3é5 ewe 


*¥X OVERALL CREAL + IMAG) STATISTICS FOR BAUD TYRE0Z 


TOTAL NUMBER OF POINTS, 2K 
BAUD MEAR 

BAUD VARIANCE 

BAUD SNROUT 


98 


64 
1647-16551 
28 26) 925 
12,0005 = ~10.7920 DE 


Par Cieotkis = 10. 000F= 10-000 DB 


DAVE YPE SS: KX ="T024 SAMPLE PUINTS; K = 64 TONES 


** GIVEN THE TRANSMITTED PHASE IS IN THE IST QUADRANT xx 


NHUMBERGOR TONES 

MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 
QUADRANT SNROUT OF DFT REAL PART 


MEAN OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 


lS 
341 .20361 
WS. OD 
Groltgio = ose sity 10) 8} 


334.22290 
(Asis era a 
14.8627 = 11.7210 DB 


xX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT Xx 


NUMBER OF TONES 

MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 
QUADRANT SNROUT OF DFT REAL PART 


MEAN OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 


15 
-406 .96606 
TOT G5 «29 
iD. s90oe— = 11.66c9 DB 


Sie ESIC BR, 
2007.02 
Styne > “lie7sc8. DB 


** GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xXx 


NUMBER OF TONES 

MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 
QUADRANT SNROUT OF DFT REAL PART 


MEAN OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SHROUT OF DFT IMAG PART 


20 
eo0 e229 5 
15794 .66 
825575. = geeceu> DE 


Soom oo 21 
PVCs 66 
PU sc2oy— — 10S tooo DB 


**¥ GIVEN THE TRANSMITTED PHASE IS IN THE 4TH QUADRANT xx 


NUMBER OF TONES 

MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 
QUADRANT SHROUT OF DFT REAL PART 


MEAN OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 


iG 
55199591 
Od Ce 
i2Ze7oco = IT.04699" DB 


=407 50591 
A264 
Po2eteg.= Jiezril DB 


XXX OVERALL CREAL + IMAG) STATISTICS FOR BAUD TYPE 3 XXX 


TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 

BAUD VARIANCE 

BAUD SHROUT 


oo 


26 
364.57397 
Pelle? Gr. 1 
Wiewoocd,- 10.7525 DB 


INPUT SHRNB = 10.000 =" 20, 000M is 


BAUD TYPE 4: KX = 2048 SAMPLE POINTS; K = 128 TONES 


xX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 


NUMBER OF WOES 

MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 
QUADRANT SNROUT OF DFT REAL PART 


MEAN OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 


38 
747 .69604 
47469.11 
ll.7771 = JO. 720GeDE 


660.96704 
43777 .00 
9. 9796 = 9 9911 


** GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xXx 


NUMBER OF TONES 

MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 
QUADRANT SNROUT OF DFT REAL PART 


MEAN OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 


34 
=f DS Sauia 
G7 932 287 
13.3044 = l162599805 


792i oo 
OOS Dene 
15.7280 = Jiggecr ne 


*¥ GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT Xx 


HUMBER OF Ones 

MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 
QUADRANT SNROUT OF DFT REAL PART 


MEAN OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 


Se, 
=7GUea6 I 9 2 
SSO 29 
10.3216 = l02is7S5are 


=/ ieee 55> 
49063.26 
10.3259 = 10.2369 e05 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 4TH QUADRANT xXx 


MUNMBER OF TOES 

MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 
QUADRANT SNROUT OF DFT REAL PART 


fEAH OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 


21 
i Lowod 500 
DUG. 4.5 
11.9296 = 10,760 saS 


=69 7/27 01 
G75 Roa27 
10.2414 = 10.1036 DB 


XXX OVERALL CREAL + IMAG) STATISTICS FOR BAUDSI YEE sGaee 


TOTAL HUMBER OF “POINTS aan 
BAUD MEAN 

BAUD VARIANCE 

BAUD SHNROUT 


100 


256 
7397607 3. 
45891.41 


1179199 = 210. 7627 ee 


INPUT SNRNB = 


10.000 
BAUDPIYFE 5: KX 


10.000 DB 
4096 SAMPLE POINTS; K = 256 TONES 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xXx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DFT 
Di 
DET 


DFT 
DET 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


64 

50s .51 250 

2o2452..51 

eee Oy := e-9095 D5 


P0@ 18559 
265130.50 
9.2518 = yeood> DB 


**% GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xXx 
NUMBER SOF WOES 


MEAN 
VARIANCE 
QUADRANT SWNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DEL 
Diet 
DET 


Dit 
Diet 
DAT 


REAL 
REAL 
REAL 


IMAG 
ItHAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


61 
ebooe.2<o514 
242012.25 
POeC0vce = 1050033 DB 


TSG555 0 W865 
oc? ol 7 266 
6.4689 = erw es DB 


XX GIVEN THE TRANSMITTED PHASE IS IH THE 3RD QUADRANT xx 


MEAN 
VARIANCE 
QUADRANT SNROUT 


iEAN 
VARIANCE 
QUADRANT SNROUT 


ee Oe 1 OES 


DFT 
DFT 
DET 


DET 
Digi 
eas 


REAL 
REAL 
REAL 


IMAG 
ItAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


So 
plas 6791.5 


169614.31 

oro = mee 1567 DB 
= Gly eo co | 
= 265248.00 
= 577 = Sma 950 DB 


XX GIVEN THE TRANSMITTED PHASE IS IN THE GTH QUADRANT xx 
HUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


Digek 
DET 
ame 


DEW 
DFT 
Die 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


ve 

1450.37061 

A0Z104 .12 

10.4084 = 10.1738 DB 


miele SOG 
Ca 2e Olde | 
7.8682 = &.958/ DB 


SAP VERCEUSCREAL + TMAG) STATISTICS FOR BAUD TYPE 5 xxx 


TOTAL NUMBER OF POINTS, 
BAUD MEA 
BAUD VARIANCE 

BAUD SHROUT 


101 


2K 


os ol2 DB 


D. INPUT SNR = 15 dB 


INPUL SHRNB = 31,5623)" 315. 000 saa 
BAUD TYPE 1: KX = 256 SAMPLE POINTS; ees 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xXx 


NUMBER OF TONES = 4 
MEAN OF DFT REAL PART = 9176996 
VARIANCE OF DFT REAL PART = 265i 
QUADRANT SHROUT OF DFT REAL PART = 31.7097 = 15.0119 DB 


MEAN OF DFT IMAG, PART = 97.37494 
VARIANCE OF DFT IMAG PART = 209296 
QUADRANT SNROUT OF DFT IMAG PART = 45.1599 = 16.5475 DB 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xXx 


NUMBER OF TONES = 2 
MEAN OF DPT REALS Paki Og. 0999 
VARIANCE OF DFT REAL PART = IS 762 
QUADRANT SNROUT OF DFT REAL PART = 588.8657 = 27.7002 DB 


MEAN OF DFT IMAG PART = 89.00999 
VARIANCE OF DFT IMAG PART = 88.18 

QUADRANT SHROUT OF DFT IMAG PART = 89.8486 = 19.5351 DB 
xX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xx 


NUMBER OF TONES = 6 
MEAH OF DFT REAL PART = SOD. 20526 


VARIAHCE OF DFT REAL PART = G07 69 
QUADRANT SHNROUT OF DFT REAL PART = 17.8737 = 12.5221 DB 
MEAN OF DFT IMAG PART = -102.36160 


VARIANCE OF DFT IMAG PART = S560 
QUADRAHT SNROUT OF DFT IMAG PART = 294.1055 = 24.6850 DB 
XX GIVEN THE TRANSMITTED PHASE IS It THE GTH QUADRANT xx 


NUMBER OF TONES = + 
MEAN OF DFT REAL PART = 98.51746 


VARIANCE OF DFT REAL PART = Z2oG798 
QUADRANT SHROUT OF DFT REAL PART = 40.6241 = 16.0878 DB 
MEAN OF DFT IMAG PART = -98.39490 
VARIANCE OF DFT IMAG PART = S6On01 
QUADRANT SHROUT OF DFT IMAG PART = 26.3941 = 14.2151 DB 


XXX OVERALL CREAL + IMAG) STATISTICS FOR BAUD TYPE 1 XXX 


TOTAL HUMBER OF POINTS, 2K = 52 ; 
BAUD MEAN = 94.6193] 
BAUD VARIANCE = 179.44 
BAUD SHROUT = 49.8941 = 16.9805 DB 


INPUT SHRRHB = 


BAUD TYPE 2: KX 


Joc 
Bla SAMPREEPOLNIS:; 


Sy UU Ne 


K = 


32 TONES 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DFT 
DFT 
DFT 


Dnt 
DFT 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


» 
ise 599 1 
469.50 
6527905 = 


Oe 9.9.) 
Ii ge pepe aE 
Ze N27 0 = 


18.0475 DB 


1329295 DB 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


*X GIVEN THE TRANSMITTED PHASE IS IN 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DFT 
Deal 
Diteal 


DFT 
DFT 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


NUMBER OF TONES 


OF 
OF 
OF 


OF 
Or 
OF 


DFT 
DFT 
DFT 


DFT 
DFT 
Diet 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


8 
= 09.995 
S99 12 
15.4151 = 


L988 
G25. 56 
48 .0524% = 


10 


mou om 
ON 
ae) 
NO 
iN 
c 


=c 0 0FIG1 991 
961.04 
41.6297 = 


Pee 795 DB 


l6émmel7i DB 


THE 3RD QUADRANT xXx 


Piece oo DB 


16.1940 DB 


**% GIVEN THE TRANSMITTED PHASE IS IN THE 4TH QUADRANT xXx 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 


NUMBER OF TONES 


DFT 
Deo 
DFT 


DET 


REAL 
REAL 
REAL 


IMAG 


PART 
PART 
PART 


PART 


9 
ioral 7 G:/ 
P27 oe 
24.2659 = 


Sy OGu 6 6.61.9 


13.8496 DB 


VARIANCE 
QUADRANT SNROUT 


DFT 
DET 


IMAG 
IMAG 


PART 
PART 


G7 VielZ 


Bao = ively 02 DB 


XXX OVERALL CREAL + IMAG) STATISTICS FOR BAUD TYPE 2 Xxx 


POTAL HUMBER @DE POLNTS, .2Kh 
BAUD MEAN 
BAUD VARIANCE 


BAUD SNROUT 
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Sie 20° = 


64 
noe. 6/762/ 


894.04 
Po. 720 be DS 


INPUT SNRNB= = 3 


BAUD WR ees: ) Kx 


1.625542 


i> 000 TDs 
= 1026 SAMPLE PORN TS aa 


64 TONES 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xXx 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


NUM E GR Gees 
DFT "REAL EAR? 
DFT REAL PART 
DFT REAL PART 


PART 
PART 
PART 


DFT 
DFT 
DEG 


IMAG 
IMAG 
IMAG 


ie) 
550 (43823 
I392es 
22.7798 = WS 25745008 
S960.55 52 
2375760 


590.5885 =" 17 0405805 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xx 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


NUMBER OF TONES 
DFT REAL PART 
DFT REAL PART 
DFT REAL PART 


DE 
DET 
DEF 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


5 
=387 oso 
S9CS785 
44.1332 = 16.4476 DB 
56 8720725 |: 
2 Sages 


G6.5359 = 169677 Jee 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xXx 


MEAR 
VARIANCE 
QUADRANT SHROUT 


MEAN 
VARIANCE 
QUADRANT SHROUT 


NUMBER OF TONES 
DET REAL PAR 
DFT REAL PART 
DFT REAL PART 


IMAG PART 
IMAG PART 
IMAG PART 


DET 
DET 
DFT 


20 
-362.88428 
G90 Sac! 
26.3725 = 2G. 21 See 
= 59:9 e829 
S49 Sea2 


34.9457 = 15.4339 DB 


XX GIVEN THE TRANSMITTED PHASE IS IN THE GTH QUADRANT Xx 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


HUMBER SOR STONES 
DET REAL eae] 
DET “REA SRAR SG 
DET READ SEAR | 


DFT IMAG PART 
DFT INAG PART 
DFT IMAG PART 


16 
SD Ow SOS 2 
eu Oon tee 3) 
Gl. oveor= 


= 587 25 2. 
I96er 6 
of 9 0Goe 


16.160 beUe 


15.7565 2DB 


xX OVERALL CREAL + IMAG) STATISTICS FOR BAUD TY PESs =. 


TOTAL NUMBER 


OF RPOLHiS ea 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 


104 


SY 74 = 


es. 
JO oe Wie 


S551564 
15.7346 05 


Titleotnine = Ssl.6255-7015.000 DB 


BAUDP TYPE G4: KX = 2048 SAMPLE POINTS; K = 128 TONES 


**X GIVEN THE TRANSMITTED PHASE IS IN THE 15ST QUADRANT Xx 


NUMBER OF TONES 

MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 
QUADRANT SNROUT OF DFT REAL PART 


MEAN OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 


| | 


wot tl 


38 
744 .14038 
148493.46 
oy ome = 15.7177 DB 


696 .014940 
TS9Mer 01 
S4260116 =" 15.4072 DB 


**X GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xXx 


NUMBER OF TONES 

MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 
QUADRANT SNROUT OF DFT REAL PART 


MEAN OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 


out out oot 


34 
-774.32231 
15403.46 
38.9046 = 15.9000 DB 


769.15747 
Tass ol 
47.2322 = 16.7424 DB 


** GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xXx 


NUMBER OF TONES 

MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 
QUADRANT SNROUT OF DFT REAL PART 


MEAN OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 


i) 
tr Omer eo. 
L6G5)) 50 
Sone 5S. 1616 DB 


=1:e5 2 UG hoc 
Ps4a7 .06 
O94. 0555e= “25-5218 DB 


xX GIVEN THE TRANSMITTED PHASE IS IN THE 4TH QUADRANT Xx 


ROUMNBGERy UF TONES 

MEAN OF DFT REAL PART 

VARIANCE OF DFT REAL PART 
QUADRANT SNROUT OF DFT REAL PART 


MEAN OF DFT IMAG PART 
VARIANCE OF DFT IMAG PART 
QUADRANT SNROUT OF DFT IMAG PART 


Ht 


wou ot 


an 
onl eae ONS, 
TD0o7 so! 
S6.7005 = V5 76467 DB 


Sl oree lo deZ 
15lS5 a0 
Sosoee0 = 15.2960 DB 


XXX OVERALL CREAL + IMAG) STATISTICS FOR BAUD TYPE G Xxx 


HOT AL SNUMBER (OP SPOLNTS> 2K 
BAUD MEAN 

BAUD VARIANCE 

BAUD SNROUT 


tow ou ou 
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256 , 
S99 055 1 
14483.24 
Seer oo = 15.7748 DB 


INPUT SNRNB = 


SL e645 


BAUD TYPEM=e= KX 4096 SAMPLE POINTS; 


iS 00 ODE 


256 TONES 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST 3@UA0DKAii tex 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DET 
DET 
DFT 


DFT 
Die 
Det 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


64 


15 O8e59 551 
94514.62 
24.0731 


= 13.6155 08 


1544.17187 
S107 Sele 
2? ~5696 


= 14.3751 DB 


xX GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT xXx 
NUMBER OF 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DFT 
DFT 
Det 


DFT 
DFT 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


LOWES 


PART 
PART 
PART 


PART 
PART 
PART 


ona 


mo to 


=1539 55403 
77501 .494 
30-59 Ge 


1477 .20898 
103902256 
21.0018 


= 14.8567 DB 


= 13. 222605 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT Xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEA 
VARIANCE 
QUADRANT SNROUT 


DET 
Del 
DPi 


Dea 
DET 
DEG 


REAL 
REAL 
REAL 


IMAG 
INAG 
INAG 


PART 
PART 
PART 


PART 
PART 
PART 


Do 


Wom ul 


=14996 5051569 
53483 .48 
41.8469 


= 16.2166 DB 


=14659 Feo! | 
S51 S62 06 
250-7 


= 14.0888 DB 


xX GIVEN THE TRANSMITTED PHASE IS IN THE 4TH QUADRANT xXx 
HUMBER OF TONES 


MEAT 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


OE 
CF 
OF 


OF 
OF 
OE 


DFT 
DET 
Det 


DFT 
Den 
DEW 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


hic 


L677 so95.25 
65498 .62 
33.3441 


+1512 osc) 
D9 COU, 
29.7169 


=" lope o0 are 


SS Fa 


XXX OVERALL CREAL + IMAG) STATISTICS FOR BAUD TYPE 5 Xxx 
TOTAL NUMBER SC RSE Oise 


2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SHROUT 


thot Wl 
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Dl 


L502 Scot 
S161 .50 


ai ols ygl 


* 


14.4412 DB 


E. INPUT SNR = 20 dB 


INPUT SNRNB = 100.000 = 20.000 DB 
BAUD MI Peewee KX = 256 SAMPLE POJHTS; K = 16 TONES 


*X GIVEN THE TRANSMITTED PHASE IS IN THE 15ST QUADRANT xXx 


NUMBER OF TONES = 4 
MEAN OF DFT REAL PART = 9a 25 7°45 
VARIANCE OF DFT REAL PART = 84.15 
QUADRANT SNROUT OF DFT REAL PART = 98.9206 = 19.9529 DB 
MEAN OF DFT IMAG PART = 94.37244 
VARIANCE OF DFT IMAG PART = 66.41 
QUADRANT SNROUT OF DFT IMAG PART = 134.1027 = 21.2744 DB 


**X GIVEN THE TRANSMITTED PHASE IS It THE 2HD QUADRANT xXx 


NUMBER OF TONES = 2 
MEAN OF DFT REAL PART = sole 7OU00 
VARIANCE OF DFT REAL PART = God 
QUADRANT SNROUT OF DFT REAL PART = ®XXXx®xX = 32.7263 DB 
MEAN OF DFT IMAG PART = 89 .66499 
VARIANCE OF DFT IMAG PART = Wd 0 


QUADRANT SHROUT OF DFT IMAG PART 288.1604 = 24.5963 DB 


®*X GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xx 


NUMBER OF TOHES = 6 
iEAN OF DFT REAL PART = Sou mo slo 7 
VARIANCE OF DFT REAL PART = 297, 06 
QUADRANT SNROUT OF DFT REAL PART = 59.4730 = 17.7432 DB 
MEAN OF DFT IMAG PART Bore LOo 


VARIANCE OF DFT IMAG PART 
QUADRANT SHROUT OF DFT IMAG PART 


|i ae. 
840.5669 = 29.2457 DB 


xX GIVEN THE TRANSMITTED PHASE IS IN THE 4TH QUADRANT xx 


NUMBER OF TONES = 4 
MEAN OF DFT REAL PART = 95.00745 
VARIANCE OF DFT REAL PART = ge eae 
QUADRANT SHROUT OF DFT REAL PART = 119.8836 = 20.7876 DB 
MEAN OF DFT IMAG PART = —-94.95740 
VARIANCE OF DFT IMAG PART = FG 2S 8 
QUADRANT SNROUT OF DFT IMAG PART = 77.4774 = 18.8917 DB 


XXX OVERALL CREAL + IMAG) STATISTICS FOR BAUD TYPE 1 Xxx 


fOTAL NUMBERM UE SPOLIT Steck = 52 
BAUD MEAN = O2nGel4/ 
BAUD VARIANCE = Diora? 9 
BAUD SUROUT = 152.7199 = 21.8104 DB 
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DNPUT SSR B = 1002000 


BAUD TV REw a x 


512 SAMPLE PORTS: 


Z0 0002p 


K = 


52 WOES 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


Dia 
DFT 
Dew 


DFT 
DFT 
Din 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


i ee 


5 
176.45985 


147.68 


210.8442 
189 2035990 


= 23.2396 DB 


549.10 


6521219 


= 18 137207 


** GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT XxX 
NUMBER OF 91 OHIES 


MEAR 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DFT 
Dit 
DEF 


DFT 
DFT 
Dri 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


3) 


=1752549995 


De ao 0 


51.2051 
176 .44992 


= 17.09806DB 


ie. co 


157.6350 


= 2159820 BE 


** GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DFT 
DEe 
DFT 


DFT 
Dial 
Dew 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


10 


= Odeo 7 2-0 


197-205 


167.6235 
+19 lvoe 


= 22.2485 DB 


394.44 


120235977 


= 20.8134 DB 


** GIVEN THE TRANSMITTED PHASE IS IN THE GTH QUADRANT xx 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


NUMBER OF TONES 


OF 


DFT 
DFT 
DET 


DE 
DFT 
DEm 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


y 


17 Seo 


4055.12 


18 35936 


= 18.9428 DB 


=184 21101 
PAA 
160.4063 


32) 
= 22.0522 905 


X*¥*X OVERALL CREAL + IMAG) STATISTICS FOR BAUD Ti eee. 


TOTAL NUMBER OF SPORTS, 
BAUD MEAN 
BAUD VARIANCE 
BAUD SNROUT 
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2K 


Wout ot ot 


64 


181.85138 
WA TE 


Ilo w> 71 


» 


20.6803 DB 


INPUT SNRNB = 100.000 


BAUD TYPE 3: KX 


1029 SANPLE BOINTS; 


20.000 


DB 


K 


64 TONES 


**X GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 


QUADRANT SNROUT’ 


DET 
DFT 
Digal 


DRT 
Biel 
DET 


REAL PART 
REAL PART 
REAL PART 


IMAG PART 
IMAG PART 
IMAG PART 


Le 


ee 


Ue 


Se 6.2256 
Peal, 72 
74.3160 = 18.7108 DB 


SO Smo 5e) 
749 
166.9495 = 22.2257 DB 


~ £4 


*X GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT Xx 
NUMBER SGI I ORES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DET. 
DET 
Diet 


DET 
DFT 
DET 


REAL PART 
REAL PART 
REAL PART 


IMAG PART 
IMAG PART 
IMAG PART 


15 


-~376.68604 
1074.68 
Poe Jone 


pees. CUGG. DB 


SO ao wewoce 


2 eee 


144.9895 


= 2) .0459. DB 


XX GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xXx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


TE 
Biel 
DFT 


DED 
DFT 
DFT 


REAL PART 
REAL PART 
RE Atom eet 


IMAG PART 
IMAG PART 
IMAG PART 


ZU 


= 562007 5 
Lio 
83.3354 


=e oe Z0os eB 


= Do ee Ao 
HM) eae 
114.5647 


= 20.5905 DB 


**¥ GIVEN THE TRANSMITTED PHASE IS IN THE STH QUADRANT xXx 
NUMBER GE TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


ie 
ams 
DET 


DN 
DET 
DET 


REAL PART 
REAL PART 
REAL PART 


IMAG PART 
IMAG PART 
IMAG PART 


16 


Spo 9G | 


od eo | 


Poe loa 


=5/6= 4509 | 
1248.95 
LVS. 2006 


== Z2i.220)] DB 


= 20.5500 DB 


¥*¥X OVERALL CREAL + IMAG) STATISTICS FOR BAUD TYPE Si eX 


TOTALANUME ER WOR ePOTHTS, 2K 
BAUD MEAI 
BAUD VARIANCE 


BAUD 


SHROUT 


109 


Ze 


Sein Mes es 
MS, 0 


TVS 2599 


20 27 Gap 


INPUT SNRNB = 
BAUD TYPE 4%: 


= 
oo 


100 G00 


KX 2096 SAMPLE BOUTS, 


20000505 


eS 


lZomEOniEs 


*¥% GIVEN THE TRANSMITTED PHASE IS IN THE LST SQUADRANTexx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DFT 
DET 
DFT 


DFT 
Dew 
Dint 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


38 


Wot on ont 


79221963 
4625.64 
119 e627 


(Ee ley Slt 
4457.19 
114. 9¢55—— 


20.7570 0e 


20.6041 DB 


X*® GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT Xx 
MUMBER OF ieites 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DFT 
DFT 
DP 


DET 
DFT 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


34 


=/ OU Go ke 
SOS erly 
114.6950" = 


7 55S Oe > 
Ty 


20. .60290RE 


20°06 


145.7620 = 


21.6366 DB 


*¥ GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xXx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


XX GIVEN THE TRANSMITTED PHASE IS Tii@he 


MEA 
VARIANCE 
QUADRANT SNROUT 


MEAT 
VARIANCE 
QUADRANT SNROUT 


DET 
DFT 
Dil 


Bel 
DFT 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


HUMBER OF FORIES 


DFT 
DE 
DFT 


DE) 
ie Tf 
OLS af 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


wi 


35 


fal 


~738.78418 
Scour oC 
104.3598 = 


20.18538p5 


=1 SC moe ou 
G8549.96 
110.5407 = 


20.43.52 505 


GTH QUADRANT Xx 


150.57 007 
4882.50 
115 US ove 


=/26 306567 
4864.54 
108.3701 = 


2026216505 


20). 34990 


*¥XX OVERALL CREAL + IMAG) STATISTICS FOR BAUD STP ee 
TOTAL HUMBER OF POINTS, 


2K 
BAUD MEAN 
BAUD VARIANCE 
BAUD SHROUT 


110 


256 
7902 06009 = 
4535.54 
119.4913 = 207s Ds 


INPUT SNRNB = 100.000 


BDAWEIYPE 5S: KX 


20, 000708 


G096 SAMPLE POINTS; K = 256 JONES 


**X GIVEN THE TRANSMITTED PHASE IS IN THE 1ST QUADRANT xx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DFT 
Dil 
DFT 


DFT 
DFT 
DFT 


REAL 
REAL 
REAL 


IMAG 
TMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


64 
P50S. 59062 
SEO 717 
72.9934 = 18.6000 DB 


Poole Lco 
Bod 1 lS 
Vice 905° = 718.9697 DB 


*¥ GIVEN THE TRANSMITTED PHASE IS IN THE 2ND QUADRANT Xx 
HUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIAHCE 
QUADRANT SNROUT 


DET 
DFT 
DET 


DFT 
DET 
DET 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


on 
ESO Ul oO? 25 
2296/9 .01 
Fo uae = 196585 DB 


1489.814495 
65599.58 
Gowt6oS = 16.2077 DB 


** GIVEN THE TRANSMITTED PHASE IS IN THE 3RD QUADRANT xXx 
NUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SHROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


Dial 
DFT 
Bile 


DFT 
DFT 
DFT 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
eee fi 
PAR 


PART 
PART 
PART 


3B) 
milm00.62695 
eo. 9S 
Pode 95 = 2il.iS89 DB 


=LG65 7.2676 
26.27 0.100 
ConouUe = 19 -cSs2> DB 


** GIVEN THE TRANSMITTED PHASE IS IN THE GTH QUADRANT xx 
HUMBER OF TONES 


MEAN 
VARIANCE 
QUADRANT SNROUT 


MEAN 
VARIANCE 
QUADRANT SNROUT 


DET 
DE 
Dea 


DET 
DFT 
Bet 


REAL 
REAL 
REAL 


IMAG 
IMAG 
IMAG 


PART 
PART 
PART 


PART 
PART 
PART 


fie 
1493.18042 
220849.69 
POOR 956¢c = 20.0413 DB 


-1510.18042 
2700 1.50 
formesol = 18.8586 DB 


*** OVERALL CREAL + IMAG) STATISTICS FOR BAUD TYPE 5 X®xx 
TOTAL NUMBER OF POINTS, 2 


BAUD MEAN 
BAUD VARIANCE 
BAUD SHROUT 


a 


K 


Siz 
TP 50Gaa159 1.6 
26 500.95 
Semoe2od —- 19.3107 DB 


APPENDIX D. STATISTICS OF SNRout VS. DOPPLER ANALYSIS 


A. INPUT SNR = 15 dB 


SNROUT VS. DOPPLER TEST => EPSPaah = ier 
INPUT SNRHB = S1l36c5. = V5 . HOGER DB 
BAUD TYPE 3: KX = 10246 SAMPLE POINTS>= K == seen 


XX OVERALL CREAL + IMAG) STATISTICS FOR THE Vs sAure. 


TOTAL NUMBER UF POLIS, che LZé 
BAUD MEAH = 958.1967 
BAUD VARIAHCE = 3614.89 
BAUD SNROUT = 35.4955 = 15 50S 


** OVERALL CREAL + IMAG) STATISTICS FOR THE 2ZehD bau 


TOTAL HUN BERDin eS OUT Seach 128 
BAUD MEAN = 359) SS0n 
BAUD VARIANCE = G296e25 
BAUD SHROUT = 20.065 l= 14.7803 DB 


xX OVERALL CREAL * IMAG) STATISTICS FORBES SR eo Ure 


TOTAL HUBBER. OF SPOUTS. 22h = 128 
BAUD MEAN = 353-9637 
BAUD VARIANCE = S652. 
BAUD SNROUT = 534.9958 = a of7 © (DD 


xX OVERALL CREAL + IMAG) STATISTICS FOR THE Sih Eaune 


TOTAL HUMBER OF POLRIS, 2hr= 1235 
BAUD MEAN = 360.04858 
BAUD VARIANCE = 5868.04 
BAUD SHNROUT = 33.5194 = Ho 2 52550 


XX OVERALL CREAL + IMAG) STATISTICS FOR THE 5TH BAUD Xx 


TOTAL NUMBER OF POINTS, 2K = 128 
BAUD NEAN = 344 .496604 
BAUD VARIANCE = 5204.43 
BAUD SHROUT = 22.1965 = 13.57/67 908 
XXX TOTAL OVER ALL 5 BAUDS WITH EPSILON = 0.00 Xxx 
TOTAL NUMBER OF BAUDS = =) 
MEAN SHROUT = To696 08 
STANDARD DEVIATION OF SHROUT = 0.7865 DB 


eZ 


Ove o OPP LERwiEST s= EPSILON = 0.c5 
INPUT SNRNB = Sis. 6ao-5= 15.000 DB 
PaUUeTYruso: 9 KX = L024 SAMPLE POINTS; K = 


HH 


HX 


HK 


HK 


OVERALL AGREAL + IMAG) STATISTICS FOR THE 


TOA EANUMBER OF POINTS, cK = 128 
BAUD MEAN = Ss 6mo0 005 
BAUD VARIANCE = 4178.20 
BAUD SNROUT = 30.4355 = 


OVERALL.~CREAL + IMAG) STATISTICS FOR THE 


TOTAL NUMBER OF POINTS, 2K = 128 

BAUD MEAH = 357 .837490 

BAUD VARIANCE = G7 9..k 9 

BAUD SHROUT = 27,9630) = 

OVEWINISCREAL + IMAG) STATISTICS FOR THE 
TOTAL NUMBER OF POIHTS, 2K = 12s 

BAUD MEAN = 352.09834 

BAUD VARIANCE = SG 0.9 

BAUD SHROUT = So ouoe. = 

OVERALL CREAL + IMAG) STATISTICS FOR THE 
TOTAL NUMBER OF POINTS, 2K = 128 

BAUD MEAN = ooo 7 o406 

BAUD VARIANCE = a oie 

BAUD SHROUT = 29.6199 = 


64 TONES 


IST BAUD Xx 


bq .8556—- DB 


2ND BAUD xXx 


14.9658 DB 


3RD BAUD Xx 


bounces UB 


GTH BAUD Xx 


14.7450 DB 


= UOVERAESCREAL + IMAG) STATISTICS FOR THE STH BAUD Xx 
OTA HUMBER OF, POINTS; 2K = Zo 
BAUD MEAN = 342.34058 
BAUD VARIANCE = 5904.97 
BAUD SHROUT = 726472 = 2297 7 D-eDB 
XXX TOTAL OVER ALL 5 BAUDS WITH EPSILON = 0.25 XX 
TOTAL NUMBER OF BAUDS = 5 
MEAN SHNROUT = 149.4549 DB 
STANDARD DEVIATION OF SHROUT = 0.37350) DB 
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SNROUT VS. DOPPLER TEST == EPSTEON —SeteeD 
INPUT SNRNB = S625 49— 15.000 D5 
BAUD TYPE 3: KX = 1024 SAMPLE POINTS; K = 64 TONES 


*X OVERALL CREAL + IMAG) STATISTICS FOR THE 1ST BAUD xXx 


TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = S51). Ge2l5 
BAUD VARIANCE = DoT o tod 
BAUD SHNROUT = 22.9909 = 136 VSGaDB 


¥¥ OVERALL CREAL + IMAG) STATISTICS FOR THE 2ND BAUD xx 


TOTAL WUMBER OF GEOL Sy 2 Ke 128 
BAUD WEAN = 350.0766 
BAUD VARIANCE = DOGS eo) 
BAUD SHROUT = 22,0825) = 13.4904 DB 


¥*¥ OVERALL CREAL + IMAG) STATISTICS FOR THE 3RD BAUD XxX 


TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = 346 .26807 
BAUD VARIANCE = 4283.94 
BAUD SNROUT = 27.9886 = 14.4698 DB 


*X OVERALL CREAL + IMAG) STATISTICS FOR THE STH BAUD Xx 


TOTAL NUMSER UF POINTS, 2k = 128 
BAUD MEAN = Soca 69? 
BAUD VARIANCE = 5469.51 
BAUD SHNROUT = 22.6220R- 13.5953 )08 


xX OVERALL CREAL + IMAG) STATISTICS FOR JHE SS UR ere 


TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = Jicpe iba 
BAUD VARIANCE = 7 e0GeT 
BAUD SNROUT = P57 6o 0 11.9482 DB 
x¥X TOTAL OVER ALL 5 BAUDS WITH EPSILON = 0.50 xXx 
TOTAL NUMBER OF BAUDS = > 
MEAH SHROUT = 13,9039en8 


STANDARD DEVIATION OF SNROUT 0.9112 D5 
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Siw vo, DOPPEER VEST —-- EPSILON = 0.75 
INPUT SNRNB = Sl 6250 — 15.000 305 


BAUD TYPE 3: KX = 1024 SAMPLE POINTS; K = 64 TONES 
XX OVERALL CREAL + IMAG) STATISTICS FOR THE 1ST BAUD xx 
TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = 343.56445 
BAUD VARIANCE = aon al) 
BAUD SNROUT = 16.4005 = 12.1486 DB 
RX OVERALL CREAT + IMAG) STATISTICS FOR THE 2ND BAUD xx 
TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = seo cUlG6 
BAUD VARIANCE = 7479.94 
BAUD SNROUT = PS .95re2. = 12.0229" D5 
ax OVERACE CREAL + IMAG) STATISTICS FOR THE SRD BAUD xx 
TORAE NUMBER OF POINTS, 2k = 128 
BAUD MEAN = 336 .72485 
BAUD VARIANCE = Dig 2c0'4 
BAUD SHROUT = 21.0260 = PS er oe lo 
*¥X OVERALL CREAL + IMAG) STATISTICS FOR THE GTH BAUD Xx 
POGAL NUMBER OF ROMHTS, 2K = ine 
BAUD MEAN = 342 .88135 
BAUD VARIANCE = 7409.56 
BAUD SHROUT = Po.oord = 12.0049 DB 
XX OVERALL CREAL + IMAG) STATISTICS FOR THE 5TH BAUD xXx 
NOVAL NUMBER OF POINTS, 2k = ize 
BAUD MEAN = S125 20076 
BAUD VARIANCE = INU 3) aes 
BAUD SHROUT = 11.6614 = 10.6675) DB 
Peal OVER ALL 5 BAUDS WITH EPSILOH = 0.75 Xxx 
TOTAL NUMBER OF BAUDS = 5 
MEAN SHROUT = 12.0142 DB 
STANDARD DEVIATION OF SHROUT = 0.9088 DB 
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ONROUT VS. DOPPLER TESTI -= EFSREON - Sigaim@ 
INPUT SNRNB = $136 25% 15,000 DB 
BAUD TYPE 3: KX = 10246 SAMPLE POINTS; K = 646 TONES 


XX OVERALL CREAL + IMAG) STATISTICS FOR THE 1ST Rau ee 


TOTAL NUMBER OF POINTS, 2K = 1238 
BAUD MEAN = 532.496948 
BAUD VARIANCE = 9581.14 
BAUD SNROUT = Il. 536g = 1046209 D5 


XX OVERALL CREAL + IMAG) STATISTICS FOR THE 2ND BAUD xXx 


TOTAL NUMBER OF POINTS... 2k. = 128 
BAUD MEAN = 555 256 56:2 
BAUD VARIANCE = 10027 .16 
BAUD SHROUT = Tl. TS¢37 = 10.4752 DB 


xX OVERALL CREAL + IMAG) STATISTICS GOR [HE SRDS EAer eS 


TOTAL NUMBER OF POT See 2k = 128 
BAUD MEAN = 325777 U6 
BAUD VARIANCE = BOOT co 
BAUD SNROUT = kG 9690S L177 Ser eve 


xX OVERALL CREAL + IMAG) STATISTICS FOR ThE 4 TiebAvre 


TOTAL NURBER OF SPOTS. ck = 128 
BAUD MEAN = SW ACERS eee) 
BAUD VARIANCE = POOLS 46 
BAUD SHROUT = 10.8725 = 10mS6555 05 


*¥*¥ OVERALL CREAL + IMAG) STATISTICS FOR TEs otiebAGie 


TOTAL HUMBER OF POINTS, 2K = WR 
BAUD MEAN = 310.86450 
BAUD VARIANCE = 11450.02 
BAUD SHROUT = CPG ows = 9 2es4 te 
*¥¥*% TOTAL OVER ALL 5 BAUDS WITH EPSILON = 1.00 Xxx 
TOTAL NUMBER OF BAUDS = 5 
MEAH SNHROUT = 10.4949 DB 


STAHDARD DEVIATION OF SHROUT 0.88497 DB 
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simeuievo. DOPPLER TEST == EPSTBON = 1.25 
PRPUT sonRINB = SoZ a= 15,000 DB 
BAUD TYPE 3: KX = 1024 SAMPLE POINTS; K = 64 TONES 


¥X OVERALL CREAL + IMAG) STATISTICS FOR THE 1ST BAUD xx 


TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = 318.38794 
BAUD VARIANCE = 12548 .15 
BAUD SHNROUT = 8.0786 = 90755 DB 


*¥X OVERALL CREAL + IMAG) STATISTICS FOR THE 2ND BAUD xXx 


TOTAL NUMBER OF POINTS, 2K = 126 
BAUD MEAN = SEIU 5 Castel al Ores 
BAUD VARIANCE = oo 0 57.20 
BAUD SNROUT = i707 = 8.8867 DB 


** OVERALL CREAL + IMAG) STATISTICS FOR THE 3RD BAUD xXx 


TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = SUgcool/ 
BAUD VARIANCE = econ.) 9 
BAUD SHROUT = 10.3947 = 10.1472 DB 


*X OVERALL CREAL + IMAG) STATISTICS FOR THE GTH BAUD Xx 


TOTAL NUMBER OF POINTS, 2K = 26 
BAUD MEAN = Sao 20005 / 
BAUD VARIANCE = oto Lod. 3S 
BAUD SNROUT = 7.38749 = 8.6849 DB 


Po eOVER ALE CREAL + IMAG) STATISIIECS FOR THE STH BAUD Xx 


TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = Zoe 616.2 
BAUD VARIANCE = Logoo. 31 
BAUD SNROUT = 229671 = 1 le Ge DD 
pe TOTAL OVER ALL 5 BAUDSSHITHeEPSIUGNW = 1.25 xx 
TOTAL NUMBER OF BAUDS = 5 
MEAN SHROUT = Cg 099 . DB 
STANDARD DEVIATION OF SNROUT = 0.8568 DB 


ae 


ONROUT VS. DOPPLER TEST -="EPSTEUN) = ieee 
INPUT SNRNB = SI Io) 2 1.090705 


BAUD TYPE 3: KX = 1024 SAMPLE POINTS; K = 64 TONES 
XX OVERALL CREAL + IMAG) STATISTICS FOR THE 1ST BAUD xx 
TOTAL NUMBER OF SPOIiSye ck = 128 
BAUD MEAN = 301.69111 
BAUD VARIANCE = LGUs da. 5 
BAUD SHROUT = 5.6727 = Te oot oebe 
xX OVERALL CREAL + IMAG) STATISTICS FOR THE 2ND BAUD xx 
TOTAL NUMBER OF POLDRTS, 2k = 128 
BAUD MEAN = 30S 277565 
BAUD VARIANCE = 17292.49 
BAUD SRROUT = 5.20) ae= 7.3136 DB 
XX OVERALL CREAL + IMAG) STATISTICS FOR THE 3RD BAUD xx 
TOTAL HUMBER OF POUNDS e cih= i26 
BAUD MEAH = 28/:.729i2 
BAUD VARIANCE = D7 See 
BAUD SHROUT = 7.0999 = 8.9788 DB 
XX OVERALL CREAL + IMAG) STATISTICS FOR THE 4TH BAUD xx 
TOTAL NUMBER OF POINTS, 2K = b2c 
BAUD MEAN = 29G2 5547 6 
BAUD VARIANCE = Li 2G 309 
BAUD SHROUT = a Cio = 7, 00S2 ere 
XX OVERALL CREAL + INAG) STATISTICS FOR THE Sin Eau = 
TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = 271, 16650 
BAUD VARIANCE = Li? Gwe 
BAUD SNROUT = 4.1504 = 6.1209 DE 
XXX TOTAL OVER ALL 5 BAUDS WITH EPSITMON = espe 
TOTAL NUMBER OF BAUDS = 5 
MEAN SHROUT = 125059 ie 
STANDARD DEVIATION OF SNROUT = 0.85492 DB 
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SNROUT VS. 


DOPPLER TEST -- EPSILON = 1.75 
INPUT SNRNB = 31.623 = 15. 
BAUD TYPE 3: 


000 DB 


KX = 1024 SAMPLE POINTS; K = 


—* OVEROUEMGREAL + IMAG) STATISTICS FOR THE 


KX 


KX 


TK 


KR 


TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 

BAUD VARIANCE 

BAUD SNROUT 


CVERAUE CREAL + IMAG STATISTICS FOR THE 


TOTAL NUMBER OF POINTS, 2K 
BAUD EAH 

BAUD VARIANCE 

BAUD SHROUT 


128 

282.71875 
19969. 26 
4.0026 = 


Hou ou 


OQVERALEESCREAL + IMAG) STATISTICS FOR THE 


TOTAL NUMBER OF POINTS, 2K 
BAUD MEA 

BAUD VARIANCE 

BAUD SNROUT 


UE eee CREA + IMaG) STATISTICS FOR THE 


WOTAE NUGBER OF POINTS, 2h 
BAUD MEAN 

BAUD VARIANCE 

BAUD SNROUT 


OVERAibu (REAL + IMAG) STATISTICS FOR THE 


TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 

BAUD VARIANCE 

BAUD SHROUT 


*xX TOTAL OVER ALL 5 BAUDS WITH EPSILON = 


TOTAL NUMBER OF BAUDS 
MEAN SNROUT 
STANDARD DEVIATION OF SNROUT 


64 TONES 


1ST BAUD xx 


Gri oo DB 


2HD BAUD xXx 


= 23 
= 286 .64027 
= eS 5 a7 0 
= Se, oul = 537750 DB 
SRUeBAUD. xx 
= 2s 
= 26576060535 
= 14806.48 
= G.7673 = G.7026 DB 
GTH BAUD xXx 
= Ze 
= CLEEROIO 5S 
= 21720.69 
= 3.49236 = 5.349499 DB 
Sh paw xx 
= Zs 
= 247 .10403 
= ZG 5 lee5)5 
= 2.8964 = G- 5450. D5 
eek KH 
= 5 
= Sos . UB 
= (aco, DS 


1 


SNROUT VS. DOPPLER TEST == EPSTEO}H =e 2am 
TP UT SURE = Slee coy 15.000 DB 


BAUD TYPE 3: KX = 102% SAMPLE POINTS; i=" S645 Ges 
XX OVERALL CREAL + IMAG) STATISTICS FUR THE 15a oe 
TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = 20 Wee ao7 35 
BAUD VARIANCE = 249362.90 
BAUD SNROUT = 2.8076 = 4.4834 DB 
®X OVERALL CREAL + IMAG) STATISTICS FOR THE 2ND BAUD Xx 
TOTAL NUMBER OF POINTS, 2K = i238 
BAUD MEAN = 266 .46338 
BAUD VARIANCE = COO14 ay ® 
BAUD SNROUT = 2.6479 = 4.2290 DB 
xX OVERALL (REAL + IMAG) STATISTICS FOR THE S3RDSSAUDEe 
TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = 242295509 
BAUD VARIANCE = 18343712 
BAUD SHROUT = 5.2188 = 5.0769 D6 
*¥ OVERALL CREAL + IMAG) STATISTICS FOR THE GTH BAUD XxX 
TOTAL NUMBER OF 1POMIGS 2k 128 
BAUD MEAN = 244 59 
BAUD VARIANCE = 20)) Dowa 
BAUD SHROUT = 2,5 cge = 3.64906 DB 
xX OVERALL CREAL + IMAG) STATISTICS FORMTHE SinmeaAU Dee 
TOTAL NUMBER OF (POLS) ear. 7 128 
BAUD MEAN = fare (MEIN =. 
BAUD VARIANCE = C5565 40 2 
BAUD SHROUT = 22 oS Coe 3.4884 DB 
**¥X TOTAL OVER ALL 5 BAUDS WITH EPSILOH = 2.00 Axx 
TOTAL NUMBER OF BAUDS = 5 
MEAN SHROUT = 6.1637 2D, 
STANDARD DEVIATION OF SHROUT = 0.6458 DB 


PAy 


Ttewivo. OUPPEER@ Es) --"“EPSITEON = <¢.25 
INPUT SNRNB = 31.623 = i 0 O05 
Dawe i rmemot KX = LOZ SAMPCE POINTS; K = 


KX 


XX 


XX 


KX 


KX 


OVERALL CREAL + IMAG) STATISTICS FOR THE 


TORAL NUMBER OF POINTS, 2K = 12é 
BAUD MEAN = 243.64471 
BAUD VARIANCE = Boers L 
BAUD SNROUT = 2.0374 = 


OVERALL (REAL + IMAG) STATISTICS FOR THE 


TOPAL NUMBER OF PODETSs 2K = 128 
BAUD MEAN = 246299190 
BAUD VARIANCE = Bye Bye! eae 
BAUD SHROUT = Peolrss = 


OVERALL CREAL + IMAG) STATISTICS FOR THE 


TOTAL NUMBER OF POINTS, 2K = 126 
BAUD MEAN = 2097 cc 
BAUD VARIANCE = 2ecoo. 06 
BAUD SHROUT = 2.6064 = 


OVERAEL CREAGE +  IMAGIeSTATISTICS FOR TRE 


TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = ZOU 2366 
BAUD VARIANCE = I2coo «1 6 
BAUD SHNROUT = F265)2 = 


OVenoee CREALY + JHAG) "STATISTICS FOR THE 


64 TONES 
1ST BAUD xXx 
SOI Uew B 


2ND BAUD Xx 


Zuo. DB 


3RD BAUD Xx 


4.1603 DB 


GTH BAUD xXx 


22.76). 5DB 


BTA BAUR xx 


TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = 2494 .09976 
BAUD VARIANCE = SUG OL 
BAUD SNROUT = 1.9834 = Cava DB 
7> TOTAC OVER ALL 5 BAUDS WITH EPSILON = 2.25 xx 
TOTAL NUMBER OF BAUDS = > 
MEAN SHROUT = 5. UG502DB 
STANDARD DEVIATION OF SHROUT = UZ oo= LD 
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SNROUT VS. DOPPLER TEST -=rPSSRCN = seegeee 


INPUT SNRNB = ose 2) 7= L5¢ 


000 DB 


BAUD TYPE 3: KX = 1024 SAMPLE POINTS; K = S645 0Gites 


XX OVERALL CREAL G 


XX 


KX 


XX 


KX 


TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 

BAUD VARIANCE 

BAUD SHROUT 


OVERALL CREAL + IMAG) STATISMICSSEOR ST HE 


TOTAL HUMBER Ch SPURh is cack 
BAUD MEAN 

BAUD VARIANCE 

BAUD SNROUT 


OVERALL CREAL +c1MAG) STATESTECS Seta ane 


TOTAL NUMBER OF “PO Mii soa ek 
BAUD MEAN 

BAUD VARIANCE 

BAUD SHROUT 


OVERALL CREAL + IMAG) STATISTICS GOR Si HE 


TOTAL NUMBER OF POINTS, 2K 
BAUD MEAN 

BAUD VARIANCE 

BAUD SHROUT 


OVERALL CREAL + MAG) STATISTI I CSSR Ghar BE 


TOTAL HUMBER OFS ROT isaac. 
BAUD MEAN 

BAUD VARIANCE 

BAUD SHROUT 


IMAG) STATISTICS FOR THE 1ST 


128 

240.09703 
39182.69 
Le Goa e- 


Wats 

234.67413 
56299265 
1.49398 = 


128 

2549 .86904 
26564261 
2.9952 = 


oul of ou 


128 

ooo C6 
380499.78 
eo d= 


hort out oa 


Lae 

25 5p 9229 
34706 .496 
Te oe 


*xXX TOTAL OVER ALL 5 BAUDS WHITH EPSPECH - Se2ee0 


TOTAL NUMBER OF BAUDS 
MEAN SHROUT 
STANDARD DEVIATION OF SNROUT 


> 
2. ofa D 
1 002 eDB 


LZ 


BAUD xXx 


22 cot 9 e0n 


2ND BAUD Xx 


1. SoS 08 DE 


SRD BAUD xXx 


Vater 10/3 


GTH BAUD Tx 


] 38:2 7D 


5TH BAUD Xx 


2.7400 DB 


KKK 


Samwon Vo. @DOPRPLER TEST --VERSZUEON = 2.75 
INPUT SNRNB = Sle cs = 152 010,0» DB 
BAUD TY¥FE,S: KX = 1024 °SAMPLE POINTS; K = 


XX 


XX 


KX 


KX 


XX 


OVERALL CREAL + IMAG) Sp esyees el sis 
TOTAL NUMBER OF POINTS, 2K = ZS 

BAUD MEAN = 250,07 709 

BAUD VARIANCE = 9 So 09 

BAUD SHROUT = 1.5817 = 

OVERAECREALS + IMAG) STATISTICS FOR THE 
TOTAL MUMBER OF POINTS, 2K = 128 

BAUD MEAN = 225mg 51.08 

BAUD VARIANCE = 44515.95 

BAUD SHROUT = 1.1469 = 

OVERALL CREAL 4 KiAGy otal loriCs FUR THE 
TOTAL NUMBER OF POINTS, 2K = PZ 

BAUD MEAN = 265.65649 

BAUD VARIANCE = Pea 1 Oo 2 

BAUD SHROUT = 2.2600 = 

CVERAPiGREAL telMAG) “STATISTICS FOR [FHE 
TOTAL NUMBER OF POINTS, 2K = b23 

BAUD MEAN = Zoc.09 147 

BAUD VARIANCE = 44222 .02 

BAUD SHROUT = 1.2669 = 

Verse EeGreAl += IMAG) STATISTICS FOR THE 
TOTAL NUMBER OF POINTS, 2K = 128 

BAUD MEAN = Coe. 1 7295 

BAUD VARIANCE = Se) eg ae ie) 

BAUD SHROUT = race = 


XXX TOTAL OVER ALL 5 BAUDS WITH EPSILON = 


TOTAL NUMBER OF BAUDS 
MEAN SNROUT 
STANDARD DEVIATION OF SHROUT 


thou 
J fet 
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64 TONES 
DSi ebA UD ex x 
1.9913 DB 
22ND BAUD xXx 
0.5951) DB 
3RD BAUD xXx 
S54 DB 
Ga rien WD 
lS tiz7 3. DB 
5TH BAUD xX* 
2.4119 DB 
Zoe *X 
5905 
& DB 


INPUT SNR = 5 dB 


SNROUT VS. DOPPLER TES) --— EPSISGH = Scr 
INPUT SNRNB = Sel O2ae= 5 20100 DB 
BAUD TYPE 3: KX = 102% SAMPLE POINTS; K =" o Ge Gites 


XX 


XX 


KX 


KX 


OVERALL CREAL + IMAG) STATISTICS FOR THE 1ST BAUD xx 


TOTAL NUMBER GF PORiISs 2K == 128 
BAUD MEAN = 348 .85937 
BAUD VARIANCE = 36188 .43 
BAUD SNROUT = Sooo Um— 52.26 Poa 


OVERALL CREAL + IMAG) STATISTICS FOR THE 2ND BAUD Xx 


TOTAL HUMBER OF sPODTITS 2k e= ae 
BAUD MEAN = 352.664979 
BAUD VARIANCE = G3007 219 
BAUD SHROUT = 2.67 16e— 4. 61ers 


OVERALL CREAL + IMAG) STATISTICS FOR, JHE SRB SpAUD 


TOTAL WUMBER OF POINTS, 2K = 126 
BAUD MEAH = eRe, 1 Urey lh 
BAUD VARIANCE = SO a lew! 
BAUD SHROUT = 5.0961 = G.9082 DB 


OVERALL CREAL + IMAGO STATISTICS "FUR IRE 4thebAUpe 


TOTAL HUMBER OF PORMTS 2 ch = is 
BAUD MEAN = Bye ye ao 2) S 
BAUD VARIANCE = 56077 30 
BAUD ShROUT = See = 5 LZ Sepp 


X% OVERALL CREAL + IMAG) STATISTICS FOR TRE SS ihe saws 
TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = 50522077 6 
BAUD VARIANCE = S20 soe 
BAUD SNROUT = 1 ee 2.5256 DB 
XXX TOTAL OVER ALL 5 BAUDS WITH EPSIEGHS-= 0900. 
TOTAL NUMBER OF BAUDS = 5 
MEAN SHROUT = 4.4873 DB 
STANDARD DEVIATION OF SNROUT = lL. t22q Er 
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Simuwlevso. DOPPLER @ TEST -- EPSILON = 0.25 
INPUT SNRNB = 3.162 = 5.000 DB 
DAUD TYEtsS: KX = 1029 SAMPEE POINTS; K = 


KX 


KX 


KX 


KX 


XX 


OVERAEL TC REALE + °IMAG) STATISTICS FOR THE 


TOTAL NUMBER OF POINTS, 2K = 12Zé 
BAUD MEAN = S47 = 2 5566 
BAUD VARIANCE = yf eee ig 
BAUD SNROUT = 3.2294 = 


OVER ES CREM +°IMAG) STANISTICS FOR THE 
TOTAL NUMBER OF POINTS, 2K = 128 

BAUD MEAN = Soe 0.67 

BAUD VARIANCE = 43143.09 

BAUD SNROUT = 2,007 4 = 


OVERAEE CREAM + IMee) STATISTICS FOR THE 


TOTAL NUMBER OF POINTS, 2K = Wes 

BAUD MEAN = 37741895 

BAUD VARIANCE = $5979.64 

BAUD SNROUT = 2.0896 = 

Vern ieCREALS + =iMaG) STAPISTICS FOR THE 
DORAL HUMBER OF POINTS, 2K = 128 

BAUD MEAN = 352.7494938 

BAUD VARIANCE = 39246 .14 

BAUD SHROUT = On b7 USe = 

UVEReet (REAL + IMAG) STATISTICS FOR THE 
TOTAL NUMBER OF POINTS, 2K = ras 

BAUD MEAN = 303.096494 

BAUD VARIANCE = DOD 0412.99 

BAUD SHROUT = 1.7154 = 


=~ TOTAL OVER ALL 5 BAUDS WITH EPSILON = 


TOTAL NUMBER OF BAUDS 
MEAN SNROUT 
STANDARD DEVIATION OF SNROUT 


How 
e UO 
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6% TONES 


1ST BAUD Xx 


sy te 30), 


2ND BAUD xXx 


4.5597 DB 


SRD BAUD XxX 


4.8992 DB 


GTH BAUD xXx 


Se Ne BD, 


5TH BAUD XxX 


2.3436 DB 


0.25 Xxx 


SNROUT VS. DOPPLER TEST == ERSWeON: = aires 
INPUT SNRNB = Selo = >. VOOM) EB 
BAUD TYPE 3: KX = 10294 SAMPLE POINTS3 K = SiaeguOiies 


X*% OVERALL CREAL + IMAG) SST ues FOR THE ISI SBAUDEA 


TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = 3492.34521 
BAUD VARIANCE = S92 On ou 
BAUD SNROUT = 2.9890 = G.7/553 DB 


X¥ OVERALL CREAL + IMAG) STATISTICS FOR THE 2ZHBREAUD 


TOTAL NUMBER OF POINTS, 2K = 126 
BAUD MEAN = 346 .349448 
BAUD VARIANCE = 44183.36 
BAUD SHROUT = 2./149 = 4 3276 DE 


xX OVERALL CREAL + IMAG) STATISTICS FOR THE SRD BAUD 


TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = S27 oO Sone 
BAUD VARIANCE = 56 lSGnec 
BAUD SNROUT = C. 70008 — 4.7260 DB 


xX OVERALL CREAL + IMAG) STATISTICS FOR THE SGT eAUDe 


TOTAL WUMBER OF URDU iiiSee 2 Ks= 128 
BAUD MEAN = 346 .97949 
BAUD VARIANCE = G0537.41 
BAUD SNROUT = 2.90 - G.7275 05 


x% OVERALL CREAL + IMAG) STATISTICS FOR THE STNREAUDe. 


TOTAL NUMBER OF PODMI Secs = 128 
BAUD MEAN = 20 GeO SOc 
BAUD VARIANCE = DS Deon 
BAUD SNROUT = 1.5849 = 2.0000 DE 
XXX TOTAL OWER ALL 5 BAUDS WITH EPSIDGie= (05505 == 
TOTAL NUMBER OF BAUDS = : 
MEAN SNROUT = S209 5.005 
STANDARD DEVIATION OF SNROUT = lL oayeebDb 
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SNnGwamevo, DOPPLER TEST -- “EPSILON = 0.75 
INPUT SNRNB = S, 162° 5.000 DB 


BAUD TYPE 3: KX = 10246 SAMPLE POINTS; K = 64 TONES 
=x OVERAEESCREAL + IMAG) STATISTICS FOR THE 157 BAUD Xx 
TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = So Gege ls 77 
BAUD VARIANCE = Seo 50 
BAUD SNROUT = E07 51. = 4.2735 DB 
A SCVeraU Meer EAL + IMAG) STATPSTICS FOR THE 2ND° BAUD Xx 
TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = 338 .46997 
BAUD VARIANCE = G6133.04 
BAUD SNROUT = 2.49833 = oo 5 0S ae0B 
XX OVERALL CREAL + IMAG) STATISTICS FOR THE 3RD BAUD xXx 
TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = Sb6709514 
BAUD VARIANCE = $0676.21 
BAUD SHROUT = 2.7438 = G7 oco>° DB 
—* CUvVERALL CREA + IiAG) STATISTICS FOR THE GTH BAUD xx 
TOWAL HUMBER OF POINTS, 2k = 128 
BAUD MEAN = Sof no 5.55 
BAUD VARIANCE = 42550.491 
BAUD SNROUT = 2.0107) = 4.2776 DB 
XX OVERALL CREAL + IMAG) STATISTICS FOR THE 5TH BAUD xx 
TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = 266 .04077 
BAUD VARIANCE = Sy art 
BAUD SHROUT = 1.4129 = eel DB 
XXX TOTAL OVER ALL 5 BAUDS WITH EPSILON = 0.75 XxX 
TOTAL NUMBER OF BAUDS = > 
MEAN SNROUT = oor 7) 2. 0B 
STANDARD DEVIATION OF SNROUT = Perc cno ew. 
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SNROUT VS. DOPPLER TEST == 9EF SeRgii serene 
INPUT SNRNB = 3162 = 5.000 DB 


BAUD TYPE 3: KX = 1024 SAMPLE POINIsS>9 kK = Seoce @mes 


XX OVERALL CREAL + IMAG) STATISTICS FOR (Re Ist fauie 
TOTAL NUMBER OF POINTS, 2K = 126 
BAUD MEAN = SCI ee) 
BAUD VARIANCE = 44902 .68 
BAUD SNROUT = aR ila SS 3.6647 DB 
XX OVERALL CREAL + IMAG) STATISTICS FOR THE 2RD BAUD xXx 
TOTAL NUMBER OF ROTI ES 2 128 
BAUD MEAN = S270 S 
BAUD VARIANCE = 48909.04 
BAUD SNROUT = 2.1964 = >-4170° DB 
XX OVERALL CREAL + IMAG) STATISTICS FUR THE Srv sAUle 
TOTAL HUMBER OF POLI Ses 2k = ie 
BAUD MEAN = 20 Sel 09S 
BAUD VARIANCE = Sod Ss has Ja 
BAUD SHROUT = 2.44910 = Bete) 5.7 SNE 
XX OVERALL CREAL + IMAG) STATISTICS FOR THE GTH BAUD xXx 
TOTAL NUMBER OF SEO t tic — 128 
BAUD MEAH = 329.55516 
BAUD VARIANCE = Q517 iae2 
BAUD SHROUT = 22 ooo = 3.6760 DB 
XX OVERALL CREAL + IMAG) STATISTICS FOR THE STH BAUD Xx 
TOTAL HUMBER OF VROLTHT Sec e— 12é 
BAUD MEAN = 271.624976 
BAUD VARIANCE = GUS oT aa7 
BAUD SHROUT = 1.21/77 = 0.8555 205 
XXX TOTAL OVER ALL 5 BAUDS WITH EPSIUGN = Si gue 
TOTAL HUMBER OF BAUDS = 3 
MEAN SNROUT = Si 0.9 Pee 
STANDARD DEVIATION OF SHROUT = Yy2o4l7 v5 
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See o. DOPPLER TEST =="EPSILON = 1.25 
INPUT SNRNB = Slo ce= 5.000 DB 
PAwOTpVPews: KA = 1024 SAMPLE®@POINTS; K = 64 


XX 


KX 


KX 


KK 


OVemaseeeCReAL + IMAG) STATISTICS FOR THE IST 


TOPALY NUMBER OF POLHTS. 2ke= Zs 
BAUD MEAN = S05 105 
BAUD VARIANCE = G8664.34 
BAUD SNROUT = 1.9627 = 


OVeMreeencnl + IMAG) STATISTICS FOR THE 2ND 


TOTAL NUMBER OF POINTS, 2K = ize 
BAUD MEAN = SG en oas.S 
BAUD VARIANCE = peoOCumog 
BAUD SNROUT = OF 9 0 


OVER Ee CREAR IMAG) STATISTICS FOR THE 3RD 


POVAL HUMBER OF PORNTS.. 2k = 128 
BAUD MEAN = Ze omOos 7c 
BAUD VARIANCE = De 2 eo 
BAUD SNROUT = 2.0861 = 


OVERALL CREAL + IMAG) STATISTICS FOR THE 4TH 


VORAE HUMBER OF POENTS, 2k = bao 
BAUD MEAN = Boi \Uesrene k 
BAUD VARIANCE = G8482.35 
BAUD SNROUT = I Sy os 


TONES 


BAUD xXx 


Cu92ooe0B 


BAUD XxX 


Cro Cab 


BAUD xXx 


Se 9o2DE 


BAUD XxX 


ees (2! Oe) b a 


0.04609 DB 


KR 


XX OVERALL CREAL + IMAG) STATISTICS FOR THE STH BAUD xx 
TOTAL NUMBER OF POINTS, 2K = 128 

BAUD MEAN = 253.492087 
BAUD VARIANCE = GooZ0ruc 
BAUD SHROUT = ie oogS. = 

ae rOVAEROVER ALL 5 BAUDS WITH EPSILON = 1.25 

TOTAL NUMBER OF BAUDS = 5 
MEAN SNROUT = 2.3646 DB 
STANDARD DEVIATION OF SNROUT = ee S10 CL Os 
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SNROUT VS. DOPPLER TEST -= ERS tig) — See 
INPUT SNRNB = 3.162 = 5.000 DB 
BAUD TYPE 3: KX = 1024 SAMPLE RODENTS; K = oGuioiies 


XX OVERALL CREAL + IMAG) Si oaics FOR THE 1ST BAUD xx 


TOTAL NUMBER OF POTTS 2 — as) 
BAUD MEAN = 29 yao 1 
BAUD VARIANCE = a Yay Boyd a 
BAUD SNROUT = 6 135 2.077 Gwe 


XX OVERALL CREAL + IMAG) STATISTICS FOR THE 2ND BAUD xx 


TOTAL NUMBER OF SROLIS se2h a= 128 
BAUD MEAN = 298 .09663 
BAUD VARIANCE = DOO 7G aac 
BAUD SHROUT = Toke = L..935S5720E 


XX OVERALL CREAL + IMAG) STATISTICS FOR THE 3RD BAUD xXx 


TOTAL NUMBER OF POINTS, eK = 128 
BAUD MEAN = Zo) obey 
BAUD VARIANCE = G2304.40 
BAUD SNROUT = Poe 2.3340 DB 


XX OVERALL CREAL + IMAG) STATISTICS FOR THE STH BAUD Xx 


TOTAL NUMBER OF POINTS, 2K = lie 
BAUD MEAN = 20 en gu2l 
BAUD VARIANCE = eyalsn isle, 
BAUD SNROUT = iS 2oa- 2.0216 RUE 


XX OVERALL CREAL + IMAG) STATISTICS FOR THE SihebAU Dex 


TOTAL HUMBER OF oF Oi S Zhe les 
BAUD IiEAN = Cou ao. 
BAUD VARIANCE = Syleyrei yo) FINS: 
BAUD SNROUT = 0.8452 = —0<7304. DE 
XXX TOTAL OVER ALL 5 BAUDS WITH EPSILON = 1.50 Xxx 
TOTAL NUMBER OF BAUDS = 2 
MEAN SNROUT = 125727 Sais 
STANDARD DEVIATION OF SNROUT = lL 20 Ome 


130 


clinouipevo, UOPPLER TEST -= EPSILON@e = 1.75 
INPUT SNRNB = S70 Ca= a U00) DB 


BAUD TYPE 3: KX = 1024 SAMPLE POINTS; K = 64 TONES 


XX 


KX 


XX 


XX 


| xX 


OVERALL CREAL + IMAG) STATISTICS FOR THE 1ST BAUD xXx 


TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = 216, 1 5000 
BAUD VARIANCE = D(O92 ner 
BAUD SHNROUT = 1.3420 = ect? i DB 


OVERALL CREAL + IMAG) STATISTICS FOR THE 2ND BAUD xx 


TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = 266 .96216 
BAUD VARIANCE = G9 09 75 
BAUD SNROUT = 1.34867 = Peeooe DE 


OVERALL CREAL + IMAG) STATISTICS FOR THE 3RD BAUD xXx 


TOTAL HUMBER OF POINTS, 2K = 128 
BAUD MEAN = 247 .049657 
BAUD VARIANCE = G17 2a 
BAUD SNROUT = We Ssaii = 1.3068 DB 


OVERALL CREAL + IMAG) STATISTICS FOR THE 4TH BAUD Xx 


COPA UMDEReOrerROLNT So, 2k = bac 
‘ BAUD MEAN = 261.21 61-2 
BAUD VARIANCE = 56809. 34 
BAUD SHNROUT = ie 25) 5 a= 0579950 D5 


OVERALL CREAL + IMAG) STATISTICS FOR THE 5TH BAUD Xx 


TOTAL NUMBER OF POINTS, 2K = 128 
BAUD MEAN = 241.30634 
BAUD VARIANCE = 105¢5.50 
BAUD SHROUT = 0.8278 = miro 210 DB 


pe TOTAL OVERTAEE 5 BAUDS WITH EPSILOH = 1.75 Xxx 


TOTAL NUMBER OF BAUDS = 5 
MEAN SNROUT = Oe ei > DB 
STANDARD DEVIATION OF SNROUT = G7 921l6 DB 


heh 


APPENDIX E. STATISTICS OF THE DOPPLER ESTIMATION 
ALGORITHM 


A. INPUT SNR = 10 dB 


ALPHA-HAT VS. ALPHA -> INPUT SNK = lORDE 
BAUD TYPE 3: KX = 1024 SAMPLE POINTS; K = 64 TONES 
Ax EPS PEOIee=3 0.0 DELTA~ALPHA = 0.000376506 
Ee ALPHA ALPHA-HAT EPSILONSRE 
] 0. 000 0.00008162 0.21678060 
Fa 0.000 =0 7000055 ie -0.14657104 
3 0.000 0.00000018 0.00046881 
u 0.000 =( 00029125 -0.77353498 
5 0.000 -0.00000678 =0)..01.8, 056.26 
NUMBER OF BAUDS 5 


MEAN OF EPSILON-HAT 
VARIAHCE OF EPSTEOH SHAT 


-0.196174457 
0.140806437 


AX EF SLEOH = 20.425 DELTA-ALPHA = 0.000376506 
iE ALPHA ALPHA-HAT EP ST EUHSHAT 
] 0.000094128 O-0 0017936 0.497636899 
Va 0. 0000941 2s 0. 00003786 OFLU 0S6 Ua 
3 0.000094123 OF GCOOLRU2S 0.29278034 
4 0.00009412é6 -0.0002018% =0 5360657 2 
> 0.000094122 0.00009492 0.25209652 
NUMBER OF BAUDS = 5 
MEAH OF EPSTPONSHAT = On yaiasi ss 
VARDTANCE OF EPSTUOH=HAT = “Ol ISi2y ceo: 
xX EPSILON = 0.50 DELTA-ALPHA = 0.000376506 
ie ALPHA ALPHA-HAT Eno LUO HAT 
] O.,00018S255 0200 0Ti 252 =05 29886555 
Z 0.000188255 0.000132492 0. 35169344 
3 O2000T Se 255 0. 0002361 0,58 S272e7 
4 O-000Tso25s5 0.2.0.0 Olele7l 0. 311 0aSe6 
5 0 20@0186255 0.00019830 0.52667250 
NUMBER OF BAUDS > 


MEAN OF EPSILON-HAT 
VARITAHCE OF “ERSTUOHSEAd 
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0. 299766512 
0). 1 23232ne2 


ACenemnnah V>. ALPHA -= INPUTSSNR = 910 DB 
BAVUmIre 5: KX =9F024 SAMPLE POINTS; K = 64 FONES 


ex TEP ot LON =20 475 DELTA-ALPHA = 0.000376506 
ae ALPHA ALPHA~HAT EPSILON-HAT 


mUe 02825535 -0.00024254 ~0.69918107 
1000282383 0.00003668 0.097491002 
.000282385 0% OO es 307 0.35349494987 
.000282383 C0021 555 C257 2510 HOC 
“000262563 02 0C0s0crT I Ceo s7900L 


in DNWNFe 
ooceo 


NUMBER OF BAUDS 
MEAN OF EPSTLON-HAT 
VARTANCE OF EPSILON-HAT 


5 
0.236632884 
0.310890456 


XX EPSILON = 1.00 DELTA-ALPHA = 0.000376506 
Lie ALPHA ALPHA~HAT EPSILON-HAT 
th 0.000376506 ~Q0.00034564 =O; 9 303193 
Z 0.000376506 OS 0G0TS37¢ O2959578 559 
S 0.000376506 0.00023758 Oso 200257 
Gq 0.000576506 O2C0GSt6E79 0.84140372 
5 0.000376506 O2C00rg TS> Wi 5037 5167 


NUMBER QF BAUDS 
MEA® OF EPSILOH-HAT 
VARIANCE OF EPSILON-HAT 


5 
Ogos vocir al 
0.48291354994 


poe eo ON = 11.25 DELTA-ALPHA = 0.000376506 
iE ALPHA ALPHA-HAT EPSP BONEHAT 
] 0.0004970638 -0.000881495 Se vod 5505 
Z 0.000470638 O07. U0G23501 O76 2618530 
3 0.000470638 0.000549109 Oe7U a7 ooo! 
q 0.000470638 One OMiZe > | 0.61489588 
S 0.0004970638 aU .0UCH So07 Uso Soe ICO" 


NUMBER OF BAUDS 
MEAN OF EPSILON-HAT 
VAeanee UE EPSP UON=HAT 


5 
SU TU Sea07 0 
ie oo 46 
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ALPHA-HAT VS. ALPHA ==> INPUT “Sith e=s0mpsS 
BAUD TYPE 3: KX = 1029 SAMPLE POINTS; KK = sGeeOncs 


XX -EPST EON =e oo DELTA-ALPHA = 0. 900576506 
ee ALPHA ALPHA-HAT EPSILON-HAT 
] 0.000564765 =O. 0010067. -2.674496518 
2 0.000564765 -0 (00081077 -1] 55668050 
3 0.000564765 -0.00061805 “1 64) 55102 
% 0.000564765 =O. OID 06107 =O rh6.4 59 Feo 
5 0.000564765 -0.00002642 -0 070MioG 11 
NUMBER OF BAUDS = 5 
MEAN OF EPSILOWSHAT = =letelavos? 
VARIANCE OF EPSILON-HAT = Poleo lus 
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B. INPUT SNR = 15 dB 


ALPHA-HAT VS. 


KX 


XX 


XX 


BAUD TYPE 3: 
ERSLEON = 0.0 
aE ALPHA 
z Oe alee 
Zz 0. OGD 
3 0.000 
4 0.000 
5 0.000 


NUMBER OF BAUDS 
MEAN OF EPSILON-HAT 
VARIANCE OF EPSILON-HAT 


ERom EON. = 0525 
ae ALPHA 

1 0.000094128 
2 0.000094128 
5 0.000094125 
4 0.000094128 
5 0.000094128 


ERometrH 


ee AL 


tn IN WN AD 
Srioierc 


NUMBER OF BAUDS 
NEAL OF EPSELON=HAT 
VARIANCE UR EPS EEOH HAT 


Oy 510) 
PHA 


POUOMNSS 25D 
~00 6168255 
- VIOUS 255 
-UCO185255 
-UUOTSE255 


HUMBER OF BAUDS 
MEAN OF ERS LEON =RAT 
VARIANCE OF EPSILON-HAT 


DELTA-ALPHA = 
ALPHA-HAT 


0. 
= 08 
=O: 
=0F 

OF 


DELTA-ALPHA = 


ALPHA -- INPUT SNR = 


00004005 
00002999 
00000758 
000034988 
00000226 


ALPHA-HAT 


SQ) oe) t=, 


DELTA-ALPHA = 


OO Os72 15 
.00006437 
»000094953 
- 00005696 
SOUS 0 


ou a 


ALPHA-HAT 


ooo°oo 


.00023408 
.00016005 
O0GTS 7.05 
WIC Sas} 
.00020043 


Oise2 
O17 
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Le B 


KX = 1024 SAMPLE POINTS; 


Ox 
Er 
0: 
= 0% 
=O 


=i0r 
ON 


0. 
EF 


0 
0 
0. 
0 
0 


OF 


K = 


000376506 
SILON-HAT 


P0637 712 
07964146 
O2Z0929:77 
09264153 
00600086 


5 


-0.160069466E-01 
Ores 50USa7e5E-02 


000376506 
SILON-HAT 


. 364934567 
slo 901 


Z5u06921 


ole? 569 
~coro 1760 


5 


0.241130829 
U.7 29961 256E-U02 


000376506 


ERSTLON=har 
O6-2e7 1990 


Soee, 


007701 
Joloud 


HzZo09770 
Bb) DSS ey age 
.40144742 
Joe Jo2c00 


5 
tao 
(ale = 02 


64 TONES 


ALPHA-HAT VS. ALPHA -=SENPUT Sonn ieee Gb 
BAUD TYPE 3: KX = 1024 SAMPLE POINTS; K = 64 TONES 
XX EPSILON «sie DELTA-ALPHA = 0.000376506 
LL ALPHA ALPHA-HAT EPSILONSEAi 
l 0 S00 C2Ze2ae.5 0. 00035001 0.87650g55 
2 0. 0002623525 0.00025643 0.6810/72e3 
5 0.. 000 ZS 256.5 0 Oi0z 39 02 0.79417539 
4) 0. OU UZS25e.5 0.000249693 0.65530 5579 
5 0 0UUZS2565 OO (OU abe es) 0.79678440 
NUMBER OF BAUDS 5 


XX 


XX 


MEAN OF EPSILON-HAT 
VARIANCE OF EPSILON-HAT 


EPSILON 
ae AL 
if 0 
Z 0 
2 0 
7 0 
2 0 


EPSILEOR 


ee AL 


in DN IAQ F 
5 fi cam JE cams fil ca | ca 


= 0 
PHA 


00 


.000301205 


.00UsS0T 205 


NUMBER OF BAUDS 
MEAN OF EPSILON-HAT 
VARIANCE OF EPSILON-HAT 


= 0465 
PHA 


- 000320030 
000 SZ00S0 
000 S200 0 
000320050 
.000320030 


NUMBER OF BAUDS 
MEAN OF EPSTEGHN=nak 
VARIANCE OF EPSILON-HAT 


DELTA-ALPHA 


ALPHA-HAT 


0.00034829 
0.00027 502 
OF 
0 
0 


00031849 


.00026549 
-UU0 STZ 


ORs 
v.26 


] 
> 


DELTA=ALRTA 


ALPHA-HAT 


oo 2 2 & 


-U00S6E756 
.000294966 
.U00S3S904 
-UUU0Ze515 
COCs 50a 


wou ot 
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0. 


0.760878146 
0 6c59 esr e-Uze 


000376506 


EPSI LORS Hea 


oo 2 © © 


Aas Veo) Nese) 
.730494205 
odoo gon 
TUS ie2o5© 
rOd7 20097 


5 


0824394 
1642002E-02 


Oy 


000376506 


EPSILON-HAT 


(SSS ey (=) 


“976 25ton 
of ©OCOMG 
.90049165 
-1 57359798 
9 0107 Tee 


5 


0 S635 55iea2 
0. 8325092435 sua 


Ce dhe VS. ALPHA == INPUISHR =e15. DB 
BAUD TYPE 3: KX = 1024 SAMPLE POINTS; K = 64 TONES 


KX (ERS ILONg@=00. 9.0 DELTA-ALPHA = 0.000376506 


LL ALPHA ALPHA-HAT ERSTE ON SHAT 

I 0000555855 0.00038678 Viera 02 

2 0.000338855 000051935 0.83485985 

S O. 005358855 O7CU0Ss9 51 0.95485568 

uy 0. 000335655 0.00030481 0280957507 

2 OF CUUSS5655 UC 00S 59:59 O2959> 286 
NUMBER OF BAUDS 5 


MEAN OF EPSILON-HAT 
VARIANCE OF EPSILON-HAT 


0.9162249658 
ORCS 226 PISS e- UZ 


wow ou 


xx EPSILOMe= 0.95 DELTA-ALPHA = 0.000376506 


Cie ALPHA ALPHA-HAT EPSILON-HAT 

] UnmU00S5 7/651 0.00020686 0.54943258 

Fi 0.000357681 0 200035396 0.88705266 

3 07000557661 0, 00057992 1.00906944 

4 O 2000557651 0.000352449 0 660eSo19 

5 UUs 57 661 0.00037949 ] -Q0792Zes 5 
NUMBER OF BAUDS > 


MEAN OF EPSILON-HAT 
VARIANCE OF EPSILON-HAT 


0.8630653574 
(ool 207 2ee0) 


wom oat 


xX EPSTLOMe = 1.00 DELTA~ALPHA = 0.000376506 


Ee ALPHA ALPHA-HAT EPSILON-HAT 
l 0.000376506 000005566 0.09470624 
Z Or 00S 76506 0.00016103 OG 2768975 
5 C7 000376506 0.00020522 OF Sso07 101 
G OSU S76506 0.00034329 0.91177434 
5) 0000576506 0. 0004756 1 ISS 21677 
NUMBER OF BAUDS = 3) 
MEAN OF EPSILON-HAT = 0.607591391 
VARIANCE OF EPSILON-HAT = 0.148841441 
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ALPHA-HAT VS. ALPHA -- INPUT (She Soa 
BAUD TYPE 3: KX = 1024 SAMPLE POINTS; Ko= scGuRgiics 


x% EPSILON = 1.25 DELTA-ALPHA = 0.000376506 
LL ALPHA ALPHA-HAT EPSILON-HAT 
1 0.000470638  -0.00069591 -1.31713200 
2 0.000470638 0.00025919 068840146 
3 0000470638 0.00030580 0°81219983 
4 00004706 38 0.00005479 0. 14553676 
5 0 | 000470638 0.00028259 0:75056285 


NUMBER OF BAUDS 
MEAN OF EPSITUOQNShAT 
VARIANCE OF EPSILON-HAT 


5 
0.2159 E2996 
0.804957628 


xX EPSILON@=si50 DELTA-ALPHA = 0.000376506 
LL ALPHA ALPHA-HAT EPSILONSHAT 
1 0.0005649765 -0.00082413 -2.18889141 
2 0.0005649765 -0.000496605 “lL. 237 65607 
3 0.0005649765 =0. 00020933 -0,. 557557 
x 0.0005649765 -0.000249733 -0 65690029 
5 0.0005649765 -0.000049868 ~QSr29s02eu 


NUMBER OF BAUDS 
1EAN OF EPSILON-HAT 
VARIANCE OF EPSILON-HAT 


5 
-0.954047918 
0. 632gmG7 32 
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C. INPUT SNR = 20 dB 


ALPHA-HAT VS. ALPHA ~-- INPUT SNR = 20 DB 
BAUDEOYFE 3: KX = 1024 SAMPLE POINTS; K = 64 TONES 


ela ON o—i0). 0 DELTA~ALPHA = 0.000376506 


PL ALPHA ALPHA-HAT EPSLCON-HAT 
1 0.000 O.00002055 0: 0:5457-59'6 
Z 0.000 —07, 0000i6 21 ~0.04305860 
3 0.000 —O 000006 31 =0 Oley 5629 
4 0.000 =O Socdig rs 5 =0, 050s576>5 
S, 0.000 O,OoUo0Si) 0.00826014 


NUMBER OF BAUDS 
MEAN OF EPSILOM=HAT 
VARIANCE OF EPSILON-HAT 


5 
=) Sooo VE >0Z 
OProZ GUGr E. U2 


a TER SLEOM =-0.2¢5 DELTA~ALPHA = 0.000376506 


LL ALPHA ALPHA-HAT EPSTUCHSHAL 

HI 0.000094128 0.00011643 0.303246392 

2 0.0000949128 O. C0C07 S15 0.207 56477 

S 0.000094128 0.00009204 0.2449446658 

4 0.000094128 0.00007348 0. er? 

5) 0.00009491285 Ommmnog959 0.26451526 
NUMBER OF BAUDS 2 


MEAN OF EPSILON-HAT 
VARIAHCE OF EPSILON-HAT 


0.2449188309 
0.209734496E-02 


xX EPSILON = 0.50 DELTA-ALPHA = 0.000376506 


EL ALPHA ALPHA-HAT EPSTLON= ter 

iL 0.000188255 O20002h26 2 0.564971968 

2 0. 000TeS2Z55 000017358 0.46102297 

2S 0.000188255 0.00019098 0.50724143 

4 COU S255 0.00016762 0.44520897 

S, 0.000188255 O20U0TI67 5 0. Sez 5269S 
NUMBER QF BAUDS > 


MEAN OF ERS EBGN= HAT 
VARIANCE OF EPSILON-HAT 


0.500145766 
Oras Sco0UE-02 
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ALPHA-HAT VS. ALPHA == SDNBVUiI==S ecu. 
BAUD TYPE 3: KX = 1029 SAMPLE POINTS? 9K = soi 


AX EP STLUN = Ur > DELTA-ALPHA = 0.000376506 


Lae ALPHA ALPHA-HAT EPSILOR=aw 

1 0 000232565 0.00030833 0.81892306 

2 0.000282383 0.00026935 0.71539066 

3 0.000282383 0.00028965 O:7693 0a 

4 0.00UZSZ5c5 0. 000267252 0.697 Zola, 

5 0. U0UZS2363 0.00029374 0.78016865 
NUMBER OF BAUDS 5 


MEAN OF EPSIEG HS AT 
VARIANCE OF EPSILON-HAT 


0.7 50cUGe5 
0.249549934964E-02 


xx EPSILONP=S un DELTA-ALPHA = 0.000376506 


Et ALPHA ALPHA-HAT EPSTUOn Sn AT 
] 0.000376506 0.00040266 L206 eS7e17 
2 0.000376506 0.00036461 0.968490972 
5 0000376506 0.00038703 1.027967 0 
4 0.0003/76506 0. 000357 26 0.9488/573 
5 0.000376506 0.00038970 1.03504181 
HUMBER OF BAUDS = a) 
MEAN OF EPSTEOHSRAtT = PREOos as 7 5 


VARIANCE OF EPSILON-HAT 0.248605061E-02 


XA EPSILON, = 1205 DELTA-ALPHA = 0.000376506 


ae ALPHA ALPHA-HAT EPSILEON SHAT 
] 0 -000395336 0.00042184 1.12041664 
ze 0: 000395256 0.000384908 1. 0201UG32 
3 0 00NS9 S556 0.00040684 1.08056164 
4 0. 000395336 0.00037666 1.00041580 
is) 0000595336 0.000490925 1. 0S67 7a 
NUMBER OF BAUDS = 2 
MEAN OF EPSILON-HAT = 1. Cote 7a ie 


VARIANCE OF EPSILON-HAT 0.248205289E-02 
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Metitemenio vo. ALPHA -= INPUT SNR = 20.DB 
BAUD TYPE 3: KX = 1024 SAMPLE POINTS; K = 64 TONES 


AA (ERS GONe= 91.10 DELTASALERA = 0..01003/76506 


LL ALPHA ALPHA-HAT EPSILON-HAT 
1 feonalelel 0.00046105 1.17141533 
2 0.000414161 0.00021188 0. 56274819 
3 0.000414161 0.00042667 1 Yesese 79 
G 0.00061416] 0.00039613 1 OSen IG 3 
5 0.000414161 0.00042883 1.13897419 
NUMBER OF BAUDS = 5 
MEAN OF EPSILON-HAT = 1.01169682 


VARIAHCE OF EPSILON-HAT 0.6499175048E-01 


ee ES MO i=) 1.5 DELTA-ALPHA = 0.000376506 


ae ALPHA ALPHA-HAT EP oO i= iA T 
l O20 00452987 0.00026110 O69 5 ro De 
2 0.000432987 0.00023134 0.61445159 
3 0.000432987 OPU0GZ 5252 U. Gru s 254 
< 0, 0004632957 0.00021900 Usoter ore 
» 0.000432987 0.000449832 1.19074154% 


NUMBER OF BAUDS 
MEAN OF EPSILON-HAT 


5 
Us 50099657 
VARIANCE OF EPSILON-HAT Urrer2 


6925225E-01 


ae EF STLOWN = 1.20 DELTA-ALPHA = 0.000376506 


EE ALPHA ALPHA-HAT Es PPC HAT 

] 0.000451812 0.00008400 Upe2esi ces 

Z O,000G 51812 0.00024993 0.66560346 

3 0.000451812 U-O002Z7115 0.72011048 

4 0.000451812 0.00004699 0.124804902 

5 ORR alot 2 0.00046693 1.249015808 
NUMBER OF  BAUDS 5 


Near ERSLEONSHAT 
VEREAHCE OF EPSILOW=HAT 


0.594397545 
Or ior 79D. 4 


nou 
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ALPHA-HAT VS. ALPHA == INPUT S@Siike= 


KX 


XX 


BAUD TYPE 3:7 UKXe= 


EPSILON == ie2s 
Ce ALPHA 


.000470638 
.000470638 
.000470638 
.000470658 
.000470638 


in DNA Fe 
ooooao 


NUMBER OF BAUDS 
MEAN (OFS Er oll ne} 
VARIANCE OF EPSILON-HAT 


EPSLLEOGH]- 1750 
ib ALPHA 

il 0.0005649765 
(a 0.000564765 
5 0.000564765 
4G 0.0005649765 
5 0.000564765 


NUMBER OF BAUDS 
MEAH OF EPs rlonrenAT 
VARIANCE OF EPSILON-HAT 


20 DB 


1024 SAMPLE POINTS; K = 64 TONES 


DELTA-ALPHA 


ALPHA-HAT 


Oaoaq Oo 


DELTA-ALPHA 


.00010312 
.00026997 
.00010148 
.00006651 
.000586496 


ALPHA-HAT 


Se 
=U 
=e. 
=(0. 
=; 


00062383 
00024093 
COCZ2IU7 
00023840 
00006303 
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= 0.0003/76506 
EPSILON-HAT 


0.27389467 
0.7 P57 eae 
0.26905 2iyea 
0.17663908 
1.29204941 


5 


0.595565367 
O.2160753951 


= 0.000376506 
EPSILON=BAT 


-]1 65653676 
-U 26550 Ue 
-0.60310274 
me ees Su lhres(0) Sos 
—O lov Sea 


5 


-O7 790 CS 3 
0, 30217 e2o7 


D. INPUT SNR = 40 dB 


AlSineenle vs. ALPHA ~= TMRUIeSHR = $0 DB 
BAUD TYPE 3: KX = 1024 SAMPLE POINTS; K = 64 TONES 


<< EPoLeEON = "0.0 DELTA-ALPHA = 0.000376506 


ee ALPHA ALPHA-HAT EPS TECHS AT 
] 0.000 0.00000226 0 00G0tSs5 
2 0.000 -0. 00000109 -0.00289558 
3 0.000 -0.0g00003¢ =0 . 0000656 
4 0.000 -0.00000142 -0 .0037626¢ 
5) Our 0.00000092 0.:0029557 1! 


NUMBER OF BAUDS 
MEANT OE lERSTEGN=HAT. 
VARIANCE OF EPSILON-HAT 


3) 
908704E-03 
832447E-04 


07156 
VG IS 


“A EPSsleON = 0.25 DELTA-ALPHA = 0.000376506 


LL ALPHA ALPHA-HAT EPSILON-HAT 

il 0.000094128 0.000096 23 0. 2bSS5c 5 

Z 0.000094128 0.00009274% 0.296314953 

5 0. vgongcl Ze 0.00009578 0.249083829 

q 0.000094122 0.00009219 0.244985534 

2 0.000094128 0.0000948/7 0.25196874 
NUMBER OF BAUDS 5 


MEANwOr CEPSILOH-HAT 
VARIANCE OF EPSILON-HAT 


0 cS 6G booy 
0.167399326E-05 


a ER SEO =="0'750 DELTA-ALPHA = 0.000376506 
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ee ALPHA ALPHA-HAT ERS DUGM=trAr 
i 0.000188255 0.00019089 0. 507200274 
2 0. 0U0GTeéc255 OOS 55 0.4997557946 
5 O- UUGPS3255 O.- 00013357 0- 50085011 
4 0. 00G1SS255 0.00018658 0.499555252 
5 0. U00TSe255 0.00018942 0.503094753 
NUMBER OF BAUDS = 5) 

Mean Or EPSEEOHN=HAT = 0.500809073 

VARIANCE OF EPSILON-HAT = 0.204471289E-04 


ALPHA-HAT VS. ALPHA -= INPUT SiiR=eetRDs 


BAUD TYPE 3: KX = 1024 SAMPLE POINTS, N= Coates 

xX EPSIUG@H - 0275 DELTA-ALPHA = 0.000376506 

tele ALPHA ALPHA-HAT EPSILON RAT 

] 0. 0002823555 0.00028547 0.756220 

2 0.000252563 0.00028194 0.74884433 

S 0 C0025 2e55 0. 00Ucoal2 0.75216144% 

q 0.000282383 0. 00026101 0.79635658 

5 0 C0UZeZsss 0, UUUZS 535 7 23 o65e> 

NUMBER OF BAUDS 2) 


KX 


MEAN OF EPSTEOis hay 
VARIANCE OF EPSILON-HAT 


EPSILON 
ae 


tin DOWN ee 
o0.0 oo 


= ]1.00 


ALPHA 


000376506 
000376506 
0005765106 
.V00S7 6506 
.000376506 


NUMBER OF BAUDS 
MEAN OF ER oO eh 1 
VARIANCE UF EPSTEOR-=AAT 


EPSIL Ghee 1.06 
ae ALPHA 

1 OO 00IS 95556 
2 0.0003955 56 
5 0 O0USe Ss sc 
“ 0.000595336 
5 OU OUS9Sss6 


HUMBER OF BAUDS 
MEAN OF (EPSTEGI SoA 
VARIANCE OF “EPSTEORShA tT 


DELTA-ALPHA 


ALPHA-HAT EPSILON-HAT 
0.00037910 1.00688934% 
0 000s 7 2679 0.997 65507 
0.00037675 1 C00Gse 
0.000374961 0.99496496 
0.00037728 UU ZU Serre 

= 2) 

= 1.000974993 


DELTA-ALPHA 


ALPHA-HAT EPSILON-HAT 
0. 000Sg618 lL 4057 Sec 
0.00039480 1. 04360105 
0.00037560 1.051249378 
OOO S S570 1.04568005 
0. U0USIG54 1.05267429 

= S) 

= 120s 4 4 
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= Q. 


= Q. 


O° 7 Sloe sai 
0. 2099997267 EauG 


000376506 


0. 2UZS2s5 Rs EsuG 


000376506 


0.2000204490E-0% 


ALPHA-HAT VS. ALPHA -- INPUT SNR = GO DB 
BAUD TYPE 3: KX = 1024 SAMPLE POINTS; K = 64 TONES 


ex EROLLON = 1.10 DELTA-ALPHA = 0.000376506 


LL ALPHA ALPHA-HAT EPSILLONSRey 
i 0.000414161 0.00041730 1.10836124 
a 0.000614161 0.00041396 P0296 71009 
S 0.000614161 0.00041488 10192490 
4 0.000614161 0.000412846 1.096490386 
5 0.000414161 0.00041541 Tes 53556 
NUMBER OF BAUDS = 5 
MEAN OF EPSILON-HAT = PevOolo res. 


VATOIAICE “OF EPSILOU-=HAT Oro Gl Sore = 04 


re ERS] LON =aiedS DELTA-~ALPHA = 0.000376506 


ee ALPHA ALPHA-HAT EPSIEOH-=HAT, 
I! 0.000432987 0.00043634 T 15691647 
2 0.000432987 0.00043304 1) NS 0MGs 6 
5 0.000432987 0.00043389 lS 2au zed 
4 0.0004329387 0. UGU 25556 UO. 62o124o! 
5 0.000432987 0.00043442 Ie esol 5 
NUMBER OF BAUDS = 5 
MEAN OF EPSILON-HAT = 1.04808331 


Vaneeee OF EPSILON=RAT OF 57149 icor=0 1 


oe ee EO = 2,20 DELTA-ALPHA = 0.000376506 


ee ALPHA ALPHA-HAT EPSILON-HAT 

1 OPO SS13812 0.00045451 EUG eRe, 

2 OPO O0SG51812 0.00045124 erlI6 G97 7 7 

5 0.000451812 OU, COUZEZh< Uo 750925 

4 U20004s1eiz Ue U0U0G299 Op lorsu2/ 3 

5 0.000451812 0.00045256 2019901 3 
NUMBER OF BAUDS io) 


Moen cOfeer sLLORSHAT 
VARIANCE OF EPSILON-HAT 


0.894494414G , 
Oc lsUs7 ols 
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ALPHA-HAT VS. ALPHA -= INPUT SSR — aime 
BAUD TYPE 3: KX = 1026 SAMPLE POINTS) =—soae ie 


KX ERSTE ON = lec DELTA-ALPHA = 0.000376506 
ae ALPHA ALPHA-HAT EPSILON-HAT 
l 0.000470638 0.00047359 1.257 35255 
2 0.000470638 0.00008576 0.227 tego 
> 0.000470638 0.00028165 0.743055 
y 0.000470638 -0. 00011694 -0 ste 
5 0.0004970638 0.00047162 AR 25), | C2 
NUMBER OF BAUDS = 5 
MEAN OF EPSILONSHAT = O.65519 77 a0 
VARIANCE OF EPSILON-HAT = 0.460538387 
xX -EPSICON “=r DELTA-ALPHA = 0.000376506 
Ge ALPHA ALPHA-HAT EPSILON-HAT 
l 0.000564765 -0.000449042 -1l 1697 S402 
Zz 0.000564765 -O-0p Cuts -0.04919485 
e) 0.000564765 -0.000244938 -0.649908123 
4 0”. 0U00>es765 = 0) CUtZ2515 - 0). 397 Fear 
5, OS 00USo 47-65 -0.00006345 =0 les o2cle~ 
NUMBER OF BAUDS 5 


MEAN OF EPSTLOMNSHAT 
VARIGHCE OF ErsicuiienAd 


—O 268537 42 
0.197466671 
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C 
c 
c 


QAM 


OM 


oe eu os ee i a 


APPENDIX F. THE SIMULATION CODE: RXSIM 


(Hi SeseiviEATION PRODUCES A PREDICTION OF THE RECEIVED SIGNAL THROUGH 


AN ACOUSTIC CHANNEL. THE TRANSMITTED SIGNAL IS MFQPSK SIGHAL. 


Re mC LO0; 1550), PHICIO0, 1 S50) jul, DELU) 
REAL FREQ, FF, ABSORP, AAUNDB, TL, AAC100,1350) 
REAL YRX, YY, TIMECG10000), YYRX(C410000) 


REAL RXKMNC0:4096), RXKP(0:4096), RXKPD(0:4096), KFREQ(0:4096) 
ReaeevELTC 100), DELFCI00) 


REN x0100), XC100) 

REAL VYC100), YC100) 

Reameve (100), 20100) 

REAL CC100), RO100), THETADC1I00) 

REAL TAULCIO0), ALPHAC100), ALPHANC100), DELALFC100) 
REAL UHATIC100), LBAUDC100), PTOTCIQ0), NOSVAR(100) 

REAL TAUKCI350), NOSAVG, HOISE 

DUTECEReeDITFECIO0), BDTOTL, NBAUDS(50), IDFT 

INTEGER KMINC100), KHAXC100), KxX(100), Mc100), KPTSC100) 
INTEGER ITIPHIC100,1350), IPHI 


COMPLEX RXItH(0:4096), RXOUTC0:4096) 
DOUBT EweeeCISTON DSEED 
CHARACTER*90 PLABEL 


USER INPUTS FOR THE SIMULATION 


PRET ECG 000) 
READC5,*) X0, YO, 20 


WRITECS0,1000) 
Rie XO = *,xX0,' YO = *,¥0,' 20 = ',20 
REE SUe * Ct 


HRITECG, 1001 ) 
READC5,*) VXAVG, VXVAR, VYAVG, VYVAR, VZAVG, VZVAR 


RITE Co0n! 001) 
HRP ECoCr DV XAVG 
HRITEC30,%) ' VYAVG 
HRITECSOe es? * VZAVG 
HRITECS0R A) * * 


",VXAVG," VXVAR 
',VYAVG,' VYVAR 
',VZAVG,' VZVAR 


', VXVAR 
',VYVAR 
',VZVAR 


HR TECos 1010 2) 
READ(C5,%*) CO, CVAR 


VRIVEGS0,1007). . 
rele Godp-o me “CO = *)C0,* CVAR = *,CVAR 
Lik EGS Ores? 4 


LIR TE Core 0 0S) 
READ(C5,%*) THETAO 


HRITECS 0, 1003) 

WRITECS0,X) * THETAO = ',THETAO 
HRITEC SO; 497% * 

MWRITEC6,1004) 

READC S32) TAP 


HRITEC30,1004) 
HWRITEC30,%) ' IAMP = ',IAMP 
HRITEC30,x%) ' ! 

INITIALIZE VARIABLES 


DOTOMS = THE TOTAL HUMBER OF BAUDS IHN THE SIGNAL 
RUS= "CUO0SEST POINT OF APPROACH 
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BDTOTL = 0 
RO = 1OXX15.0 


X(1) = XO 
YC1) = YO 
Z(1) = 20 
IF (X0 .GT. 0) VXAVG = -ABSCVAAVG) 
IF (XO .LE. 0) VXAVG = ABSCVXAVG) 
IF (YO .GT. 0) VYAVG ="-ABS Gi 
If (YO .LE. 0)  VYAVG = ABSCy Avo 
IF (20 .GT. 0) VZAVG = -ABSC(VZAVG) 
If (ZO .LE. 0) VZAVG = ABS(VZAVG) 


WRITEC6,1010) 
READ(5,%*) NPAKS 


LRITECS0, Pon) 
WRITEC30,%) © PAKS = has 
HRITEC( 30; 40 


10-FIRDTE 
NPAKS 


HRT TECSO, 
HRITECSO, 
HIRT PEGS O51 
HRI TECS Oy x9 
HMRITECS 0, *) 


— ps 


HPAKS 
METHOD ==>" METHOR 


UG ee eli gO 8) 

IF ChHEPHOD 

DO 100 
50 WRITE 


E. 0)°GO 7010 
EE. 3 GO) Or 0 
NPAKS 


READ( 


MRE ECS0; 
WRITEC30, 
LR LPECS 0; 
HRUEEC SO, 
HRITEC30,1]1 
MRITEC So.Tg 
HIRITECS0,10 
WIRITEC3SO,¥) 
HRIVECS Oy) 


0 
] 
@ 
5 
q 
5 
6 
D 
102 
102 
102 
102 
2 
z 
Va 


= =-<AWlnD GN © 
Nee ee es oe 


BDIVPEC', 17°). = 60) ie 
1 


EP eCEDI VEG) ele. 0) GO Oe >0 

Li eGo Cl) 6G). 5) 760s L0e Su 
60 HRIMECG, 1027) 1 

READ(C5,*) NBAUDSCI) 


EGU 027) et 
HRITEC30,%) * NBAUDSC(C',1,°) = * BAUDSUT) 
WRITEG@S0,4) * * 


IF CilepauDpSCI) LE. 0) GORlG {66 
100 CON PIO: 
HE TEPECGs 1.05.0) 
READCS, 43-10 I 


MRITEC30,10350) 


HWRETRECSOn 4) SDE = Geet 
RoE Gs Ort? 
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Pie BD, 


ALLOW 


Peeeoor? .EQ. 1) THEN 
WRITEC6,1090) 
READC5,%) IWNDO 


WRITEC30,1040) 


WRIFEC30,%) *§ INNDOW = ', IWNDOW 


NrECSO,x)~* * 
ENDIF 
WRITEC6,1050) 
READ(5,*) SNRDB 


WROMEGs0, 1050) 


MRETECSO, *CAL0,FI0.49)") * SHRDB = 


VR TE o0 , X) tt 
SWelive= 10. 0X*CSHRDB/1I0.0) 


LL = 0 
DO 300 I = I, HPAKS 
DOwe00 J! = 1, WRAUDSCI) 
ne) i ae 
Peer eC eee g. 1) THET 
DEETGELO =91.0 7 240.0 
DEE CUED =—290..0 
REINCLED = 68 
KMAXCLL) = 83 
fe GEL r= 256 
ENDIF 
Peeeopr yi Rec ly eQ. 2c) THER 
DER GEE = 1.057 7120.0 
DEEP GEL) = 120.0 
RiItGhE) = 135 
RMAACELE) = 166 
RACEL) = 512 
EVOL 
Demos OnrreC]) 260. 3) THER 
Peet) = 2105760 .0 
BEM Cli) = 60.0 
KMTWCh) = 269 
KMAACLL) = 332 
KXCLL) = 1024 
ENOL F 
Tech Dir PeCl) EO. 46) Thehl 
Dee PCEly = 170-7" 30.0 
DERE OME): = 50.0 
GED. = 537 
KMAXCLL) = 664 
KACLED = 20986 
ENDIF 
tmmeobine ECT EC. 5) THEM 
DEC r = 0) 7 ks. 0 
DEL ECLED = 25.0 
Kittie) = 1073 
Riera el == S25 
KXC(LL) = 4096 
ENDIF 
COURTEUE 


300 CONTINUE 
BUDTTAEEZE VARTABLES 


PI = G.OXATAH(I.0) 
BDIOL. = UL 
Eee DOT Gh. 1 00) nie 


PRE CGs coe eee KX CERROR: 
PIRITECS0,%)") XXX ERROR: TOTAL 
CO Migs 9 99 
ENDIF 
Pee Uist TO-ENCODE THE PHASES 
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ron DB 


TOTAL NUMBER OF BAUDS EXCEEDS 100 


HUMBER OF BAUDS EXCEEDS 100 


> 


Le CMe mnOD ee Ot 2) rain 


ee 0 
DO 490 I = 1, NPAKS 
HR ITE C6 Ee Geel 
WRITECG, iow 
410 READ(5,%) MTHDPK 
G 
HRITE CSC; Lite Tt 
HWRITECSU iio. 1 
HRITECS0,*%)" @ ShtHDPK =e hii DEk 
HRI VEG30 ee 
C 
IF (NTHDPK .LEs- 0) “CO shea. 
IF VCHIHDPRK .GE= 50) | GOR i Vaan 
IF OCHIHDPK 7 EQ.) Stier 
DO 430 J = 1, NBAUDS(I) 
LiGe= srt 
DO 420 K = KMINCLL), KMAXCLL) 
DSEED = CK *33257 29) 4 seer 
CALL PHASECDSEED, RUDPHI, Tene 
PHICLL,K) = RMDPHIT 
PUPHECEL,K) =. Cem 
420 CONTINUE 
G30 COHTINVE 
ENDIF 
LE SCHWiHDPRK. ~E9 4, 20> aitsete 
DO 480 J = ], NBAUDS(CI) 
i= bt 
G40 READC5,%®) NITHDBD 
C 
HKD ECS 05 tev) le oJ 
DIR TLE< 30 Petomels oy Loe 
HRIFECS0;%) ° TARHDBD. =] oT Hbow 
NRE A x: * 
C 
LF CiERDED «LE. °O)> -GUs 107440 
DF Ci aDeD..CE...5) ~COmg saat 
TF CMTIHDBD .EQ:. 1) THE 
DO) G50 K = KATUCLL),. KHAX Ci 
DSEED = CKy e527 tn Viste 
CALL PHASEGDSEEDRIIDPHL I it) 
PHlObe, Kh) = RUDE NA 
Pic bOEL eK = ert 
G50 COMM PIVE 
ENDIF 
Lea Crne eb EQ. 2) THEN 
Hi vec6, | YOO) 
MRAaTEC6, 1101) 
C 
WET TEGO 100) 
HRI PECSUS Tio) 
e 
DO) 470" Ko = KMTHCLED, KMAX CHS) 
READCS,*) IPHI 
Peet Ko) s=) CPi 
C 
NREL EC SO, tL SO) ly a 
HRIMECS0+ *)—* TERi= ie 
HR ULECSO; =) : 
® 
IF CIPH]. .LEw 0). GOST Ge4e0 
IF CEPHI..GE. SI55CU FUR 
IF CIPHI .EQ. 1) PHICULZRD = C45 0. eae oo 
IF CIPHI] .EQ. 2) PHICLE, KR) = (S500 eee 
IF CIPHI .EQ@. 3) PHICLL,K) = (-o5500 er eee 
IF CIPHI .EQ. 4G) PHICLL,K) = (€-4500)0% pee 
470 COHTINUE 
ENDIF 
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LD BS Ge Ya pep) 


iG Ce Ge) Gm tj ee Si @ i ee LS @ 


Cre? OO © 


Grits <> 


Se Pie i 1 


480 CONTINUE 
ENDIF 

490 CONTINUE 

ENDIF 


Domsov, Ll = 1, BDIOTL 
INITIALIZE VARIABLES 


Dotmere= <UL * 15.5987 )>+ DSEED 


LBAUDCLL) = LL 
COMPUTE: 


robe? =sDISTAHCE THE TRANSHITTER TRAVELED Il THE 
A=DIRECTIONSDURING THEME BAUD 


G 
) = ZRiiD 
L 


COMPUTE: 


AUSSCDSEED, VXAVG, VXVAR, ZRND) 
Pol we Owe Cea CV XCEL) x DELT CLL) 


TOLLS = DISTANCE THE TRANSHIITPER TRAVELED IN THE 
Y=DIKECHION DURING THE EETH BAUD 


CALL GAUSSCDSEED, VYAVG,VYVAR, ZRHD) 


VCE = ZRIMD 


eee Oo 1) eC oe Ca) eee eC) 3 (DEL T CLL) 


COMPUTE: 


ZOE =S DISTANCE THE IRANSIITIER TRAVELED It! THE 
C-DIRECIIGCHU DURTIG THE LLIn BAUD 


CA 
Venue ye= ZRD 
lls ae 


COMPUTE: 


LL GAUSSCDSEED, VZAVG,VZVAR, ZRUD) 
( 
Riel ee EL) = ZUR at V2 lL) ex DERTCLL 2 ) 


Coli? == INE SPEED UF SOUNDS DURING THE LIH BAUD 
CALL GAUSS(CDSEED,C0,CVAR, ZRUD) 


CCLL) = ZRIUD 


COMPUTE: 
XXX = ABSCXCLL) 
T1Y = ABSCTCULD 
ZL6) = ABSC LCL 
UG Oe Ee 
Peace CELE 


COMPUTE: 
PHeELAD Crt 


Viet) 1 eet 220% 20 X05 
ek Ua CED 


ANGLE BETWEEN THE SLANT RANGE AND ZO 


ARG =eABSCZCLL)) 7 RCLL) 
Pe welt = ACUSCARG) ¥ (180.07P1) 
VOU PH TAQ). HEN 


Pee ine TADCL Ly 


vi DECG., x ) 

%, " XX WARNING: 
WRITEC6,%) 

yA HR KK XX X 


kD Tete, =) * 
VIRITEC30, X) 


ws ' ¥* WARNING: 
WRITECS0,%) 
6 1 HK KKE NRK KH 
WRIVECSO, %)" 
ENDIF 


Pie veelosORenmEReTHAN *,THETAQ,* DEGREES 
THieereeeiveRwcAl HOT RECEIVE THE SIGHAL 


Poel omorRen tee THAI ',THETA0,* DEGREES 
THe REECE Vermeer HOT VREGEIVE THE SIGHAL 
' 


% 
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C 
c 


Mt) 


Cy Cnty Gy) Gis ee Gey Se te 1 ep ne 


el Bk 


THITIALIZE CINITS FORVeRAGhS 


TEC -EQ, 1) 
XXMIN DT 
XXMAX 
YYNITH 
YYHAX 
2201 
Z2ZNAX 
CChinn 
CCNAX 
RRMIN 
RRMAX 
THT i= 
THTMAX = 

EMD F 


SET LCIMDIStrOReCkARiis 


XXMIN) 
XXIAX ) 
Toh) 
YYNAX) 
22MIN) 
Z2ZHAX) 
CCE) 
CChiax } 
RRHNIN) 
RRIMAX) 
lle 
Ole 


THEN 
arr it 


ogee le 


be tb 08 Od ae de 

eal 

+ 
A ADAOONN NK “4K 
—] OR RN EN RIN TNS TS 
SO eet at eet bee bt et bd 
MMewewwe wwe ewe 
ee ee tae oe es es a 


XXMNIN 
XXHMAX 
YO 
YYHAX 
2ecnin 
2Z2NAX 
CChini 
CCHiAX 
RRNIN 
RRMAX = 
TA) THTMIN 
THTMAX) THIMAX 


ALL THE PHASES RAHDOMLY FOR EVERY 
500 K = 


o~ LX TR LRN LR TN FRTN TN T™ 
DADOAOANN << KOK 
LR LOLOL LO LO FRO TN TN 
Ce seme | ee pet rere ee me Fe ee 
[ae oe pee ee ee ee ee 
Nee Nee Ne ee ee ed ee ee 
= 
ee Ee be ES 
220 OM a 
row rn rT TX TX TN rN er 
ee eel see Ge Gee ee LS 
A a a ll 


L 
L 
L 
L 
L 
L 
L 
L 
L 
CL 


(THETADCLL) 
CTHETA DICED 


Ph at hd be bd Rd bed Pd 4 
WMMNADNADNADA ADNAN 


ASSIGH 
DO KNINCLLD, KRAXCLL) 

LFSC it HOD SE Cee et Hie 
DSEED = ON ee) 21676 550-+ USEED 
CALLS PHASECDSEED RNUCHI Lehi) 
PHECEE ho) =—RHDEMT 
PIP COR he) = rae 

ENDIF 


500 CONTINUE 
550 °COMLTIVE 

DO.-660° “b= BOLO mt 
COMPUTE: 


ALPHACLL) = 


AVX = VXAVG * 
BMY “= VyAN Gis 
AVZ = VZAVG * 
ALPHACUL) y=—4 AO Voost eh 2) ee, 
COLPUTLE* 

DELALECLU) = 

MCLL)D. = THE DOPPLERP CHAE es 

ALPRHAHICLLD = THE SDOPPEER CHAiiItEs 
NOTE: 


DOPPLER -COMPRESS 107 
TO THE MOVING TRANSIiiiiES 


THETADCL 
THEVADCL 


L) 
L) 


PACKET 


FACTOR 


CROLL IMEC LED 


THE MAXINUM CHANGE IN ALPHA FOR THE ETH paw 
NEEDED FOR  1HESGi hepa 


FACTOR FOR THE Lik BAUD 


ALPHAMCLL) IS USED 10 ESTIMATE THE SAP Rie 
FREQUENCY OF THE RECEIVEDRSICHAE 


DELALFCLL)D. =] -7 ©3607 AMA Cie 
ANTH.= CALPHACLL) 7 -DELALE CEL Di Urs 
AMAX = CALPHACUL) 47 DEVALECE Dias 
MOLL) = INTCCANIN + ANAX) 7 2.0) 
ALPHAMNCLL) = MOLLY * DEDALEURE 
COMPUTE: 
TAULCLU) = THE TINE DELAY DUE 
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> 


THE SIGHAL TRAVELING 


Oo ee 


BODO 2 Owe tk. 


Cy Oi eae OO) Oo) Oa OO a 


OOO 


C 
C 


(LHROvVG SI Ree MEDIUI OR CHAniWEL AND THE RECEIVER 
FOR THE LTH BAUD 


PAE Ge ae= ROLL 7 COLL ) 


Peel eG. 1) THEN 

ALFMIN = INTCALPHA(C1)) 

ALFNAX = JHTCALPHAC1)) + 0.1 
ENDIF 
Decent nAClL io LT. ALFUIND GALEMIN: = ALPHACLL) 
Peace aaACkl) .G), ALEWMAXDSAUERIAX = ALFHACLL) 
PLOTS = 0.0 
HOSVARCLL) = 0.0 

M 


DO 640 K = Kt 


COMPUTE: 
Vine RECEIVESSIGHAL SAMPLE START TIME 
DEEUIN= THE STHCHROHIZATION ERROR 
WRITEC6, P18 0} Dts 
1140 FORMAT C7 Sen PeRe Het omhiwieD RECEIVE SIGHAL START TINE *,7,2X%, 
i FOR BAUD NUMBER ',149,7,4%, 
Ao ee AND THE, TIE DELAY THROUGH THE Dea kane *,7,49X%, 
oa ELECTRUNIGS FOR MONE HMUNBER *,;14,° ... 
READ(S;4) UHATICLL), TAUKCK) 
TAUKCK) = 0.005 
Ua UL COLL) 4° CTAURCK). x% Cl + ALPHACLE))) 
UHATICLL) = Ul 
BEET == Ul > URATICEL) 
COMPUTE: 
NUCLL,K) = A RANDOM SEQUENCE OF TINING JITTERS 
Geis —  DELULZDELTCLL ID x CKACLEIACIVALPHANKLL) )) 
COMP UNE 
Pies > She ANPOITUDE OF THE RECEIVED SIGHAL, WHICH I5 
PePOR ALS Civ DS LEStlORTIALIZED AMPLITUDES ARE 
VDEstReED TC Cyl Atih- lO RAT TENUATION FACTOR 
DUES TOSTHE TRANSNISssTONSIUSS TF HORMALIZED 
Hime wes ake WON  DESTRED (I.E... TANP=0) 
PeGCLANE S00) THE 
COMPUTE: 
PSoORr = THE AMOULT OF ABSORPTION Ji DB7FT AT G DEGREES C 
AeA DEPTH OF APPROXIMATELY 3000 FT 
BeeGe= Chex MDELFC LID) 7 1000.0 
fe = oPReQx<2 0 
PESO Rpe ee CU.00S +e CNCU rE) ¢ Cl+FF) 2 
i ome ttmlrr ors (oS  100trr) ) 
a PoC MOmOU02T Ss x tr )) 7 3000 
(CVERe= te eDEPTH THE TRANSMITTER 1S BELOW THE SURFACE IN FEET 
IADER = 1000.0 
CLOUT At TE XDER + ABSCZCLL)) 
AAU DB = ~10 *¥* CABSORP 7 20.0) 
te oo ce cTOIAL = 3000.0) 7 1000.0) 
iieeGtpaceetl. 0) FACTOR = 1.02 
Dee lince ol. 0) FACTOR = 0.98 
IIFAC = ABSCIFAC) 
DeeCierrac NE. 0) THE 
DO G20" 1 = I, JIFAC 
AAUNDB = AAUNDB * FACTOR 
620 CONTINUE 
ENDIF 
ABSORP = 20.0 * ALOGIOCAAUNDB) , 
CORE UTE: 


IHCLLD, KNAXCLL) 
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oO  CHO@g 


OOO Vag 


OCGA 


c 


TL = THE TRANSMISSION LOSS IN DB DUE T0 SPHERICAL-SPREADING 


AND ABSORPTION OF THE ACOUSTIC CHANNEL 
FTL = (20.0 * ALOGIOCZTOTAL)) + CABSORPRe a7 iat 


AACLL, HG S210 0xk C201 U7 2029 
SESE 

AACLL,K) = 1.0 
ENDIF 


COMPUTE THE NOISE VARIANCE FOR THES DESIRED 


640 CONTIIVE 
KPTSCLL) = CKMAXCLLD S RMT CEL Sie 
NOSVARCLL) = CPTOTCLLD *® RXCLLI) 7 (2 08 x SSHR Te Gree 
660 CONTINUE 
Il) =20 
RXMAX = 0.0 
RAHI = 0.0 
HNOSAVG = 0.0 
HOTS ER= 2050 
DO> 700 LL = looper 
DO 690 H = I, °hxCLED 
DSEED. = (CLL 4e522 767 Soe DS EED 
CALL GAUSSCDSEED,NOSAVG, NOSVARC(LL),HOISE) 
YRX = 0.0 
NH =H = 71 
et a a rk 
TIMEC ITD v= sti 
DO 680 K = KNINCLLD, KMAXCLL) 
COMPUT es 
YRX = THE RECEIVED SIGHAL Hite AGieG: 1 He AbuVe 
PARAMETERS COMBINED FOR 1HESUine BAUD 
AT TINE = HN OVER ALL KMIN TO KHAX FREQUENCIES 
YY = AACLL, RK) me COSC CUC2ZAPI SR) EXCEEDS 
CCIFALPHAMCLL)) 7 AI+AEPRACEL) De 
CHE NWUCL Ish)" 4. PHT CMI a) 
YRX = YRX + YY 
680 CONTINUE 
YYRXCII) = RY + HOSE 
IF CYYRXCII) .GT. RXAHAX) RXAMAX = TY RAXGITD 
TF -CYYRXCIID .LT. RXAMDIDCRANIM =. Terre 
690 COUTINVE 
700 COMTINVGE 
BPiSus= 91) 
KLL = Q 
DO 710 Ut = 1, S01OTL 
KEG nD Lot hE TSE? 
710 CONTINUE 
HRITECS0> 2000) REL 
WHRITEC SU, e00 ly LL K PREC. PPh 
DO 740 tl = 1, BDITOTL 
DO 720 K = KMINCLL)O, KMAXCLL) 
WRITEC30,2010) LL, Ks PRICLU, Ko eee 
720 COMPLUHVE : 


HIDE 


Pp 
N 


BAND SIGNAL-TO-NOISE RATIO 


TOTCLLS = THE TOTAL POHER OF THE LER Bae 
OSVARCLL) = THE NOISE VARIANCE OF THE LLTH BAUD 


PTOTCLL3 = PTOTCLLY 4 © CAACEL KO xx2)7 ae 


740 CONTINUE 


C COMPUTE THE DFT OF THE OUTPUT CRECEI VER SSI Ciiew 
C 


IP CRDET - EQ. 1) SIGE 
RXKMIH = 0.0 
RANIBA = 020 
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DO 800 LL = 1, BDTOTL 


If = 0 

Melt .EO. ioe athe 
JJ =I 

lola 
ee) ot) 

ENDIF 


mel = KACLL) Sl 
Je J OK XA 
WOe7e0 “J°= JJ, J)0 
ReeriCil) = CMPUXCVYRA( J) 2020) 
Pie Ii st 1 
760 CONTINUE 
CALL DFT CKAMI, RAIN, RXOUT) 


CS CONVERT OUIFUT = DATA (10 EXPONENTIAL FORM, I1.E., RXKMCDXEXPCIXRXKPC)) 
& 


0, KXMI 

SQRTC CREALOCRXOUTC IJ) )**¥2 + CAIMAGCRXOUTCI)))xX2 ) 
RUCK AOUT CIID) CE. Ives 5) En 

= ATAH2CAIMAGCRXOUTCI)), REALCRXOUTCI))) 


PieCAbo CA TIAGUR OUI) 2) SLE. 12E-15),  RXKPCI) = 0.0 

PeeGAlLhAGCRAGURCIOD@ Gi. T,E-I15) RXKPRIDe= Pi72.0 

eee oe CR SOU RCI ei a =1.E-15) RXKPCI) = -PI72.0 
ENDIF 


Reel) = “CRXRPCL) * 180.0) 7 PI 
770 CONTINUE 


RE 
) 


e 
e TeCES OE LHRUT AND OUTPUT TO A FILE 


WRITEC30,2200) LL 
DO 780 I = 0, KXMl 
eet ou; 27 lO) el REAL CRAIHCID), ALMAGCRXINCI)) 
780 CONTINUE 
HivtveGa0, 2ez0) 7 LL 
PeeOCrHumon .EQ> 0) ° THEN 
INIH = 0 
IMAX KXM1 
MIMI 0 
MINA X 
ELSE 
INH 
IMAX 
NMIN 
WMAX 
ENDIF 
Ae GAA = — THIN) + 1 
DO 790 I = If4Itt, IMAX 
Wit oU5 2 2500) 1h REAMCRXOUTICID), AINAGORXOUTCI)), 
a Cia) eee DOCI), RARPCI) 
Dee GrnCl) Pol. RAKMIH) RARIMIH = RXKMCI) 
CE CR ele G) = RAKHAXSRARHIAX = RXAKMCI) 
KFREQCI) = I 
790 CONTINUE 


KXM1 


KMIHCLL) 
KMAXCLLD 
COL 
Cro 


tou tr al 
s 
—~——~ 
— 
ee 
— 


C 
C PROT SIRE DET GUILPUT 
C 
Beier LABEL, 2300) LL 
ieee 60 eel) CALL COMPRS 
CARES AGES 1, 625) 
CALL HOBRDR 
CALL AREA2D(8,6) 
ChE eit REQ CK) $*,100) 
CALL YHANEC' MAGHITUDE $',100) 
Sve eeiecOie ET OUTPUT OF THE RECEIVED SIGHAL ¢$',100,49,2) 
CALE HEADLHCPEABEL, 100, 3,2) 
CALL GRAFCWIMIN, 'SCALE',WMAX, RXKMIN, 'SCALE',RXKMAX) 
CARL GRIDGIF 
CARES SETCERC MAGENTA) 
Dee ChAx Tee. 100) THE 
Cr PePeCURVECKP PREG (TIN) »RAKRMOCININI, KAX, -1L) 


> 
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c 
C 
C 


CALL VLINECIMIN, IMAX, KFREQ,RXKM) 
ESE 
CALL CURVECKFREQCIMIN),RXKMCIMIND, KXX,0) 
ENDIF 
CALL ENDGR(O) 
CALL ENDPLCO)D 


CAPES PAGE IT S25? 
CALL NOBRDR 
CALL AREAZDCS, 6) 
CALL XNAHIEC™ ~FREQ] CK) 3337 
CALL YHAHEC™ SS Puaob CiEGoms 
CALL HEADTINC. DET VON TEU lr 
CALL HEADINCPLABEU, 1007s; 2) 
CALL GRAFCHIMIN, *SCALE',HMAX, - 
CALL GRIDCI], 7) 
CALE SETCURC HAGE IAS > 
IF PCKXA. CE SD ae el 
CALL CURVECKFREQCIMIND,RXKPDCIMIND,KXX,-1) 
CALL VLINECIMIN, IMAX,KFREQ,RXKPD) 
aye 
CALL CURVECKFREQCIMIH), RXKPDCIMIND,KXX, 0) 
ENDIF 
CALL ENDGRCOD 
CAL VENDPIECOD 
800 CONTIIIVE 
ENDIF 


) 
RECEIVED SIGiAiS U0 42> 
0.0,95.0, le0.0> 


FORMATS 


1000 FORMAT(7, 2X, 

% ‘PLEASE ENTER THE INITIAL POSITION 20770) 2ORGr ee 

~ ' (OF MRE RECELVER REUAIT VEs Ose TRANSMITTER ee 
1001). .FORMALCZ, 2x, 

% “ENTER THE TRAUSHITTER™S VELOCITY AVERAGE VANDA AR AiG ae, 

“% * JW THE X, Y, Z-DIRECTIONS CFIZSEC) ANDNCRT7 Sbarro 

ae (2) ik ote VXAVG, VXVAR, VYAVG, VYVAR, VZAVG, VZVAR) ') 
LOO2 VEQRiiViG7 wc x 

Zoe Te Wilok THE AVERAGE SPEED OF SOUND Ii Fl7SEGse 72 

“~-*. AHDUIRE VARIANCE IN GCEi7~SeGuxxc 4 
1003 FORMAT(CZ, 2x, 

7 ENCE THE TRAUSHITTIER™S. DON: UItik TRANSMITTED YF ne 

% * HALE BEAM PIDT He ANiGUESCbeEGs Seren 
1009 FORMATC/, 2x," DO YOU WANT NORMALIZED AMP Tat Wiig: 5 

Fon eR THE RECEIVED STOMA S47, 2X, 

sg ahs PLEASE ENVER 2s ,ES OR 0110 eee 
1010 FORMATC2X%s (ENTER THE # UR SPACKETS 1h THE TRANSMITTED SIGNAL ...1) 
LUZ0 FORMAT CA 2. 

OUVE Nes Te BAUI) TYPE # CORRESPONDING WOR ine FOLLOWING 
1021 FORHATC2ZxX,* BAUD LENGIE: FOR VPAGRE SUNDER: seo eae 
1022 FORHATC/, 10X,* 1: BAUD RENG TS Cle = 720 SECONDS : 
1023 FORHATCIOX,.' 2 > BAUD LENGIY CDELID = 4) 20s SECU 
1024 FORMATCIOX,® 3 : BAUD VENGin- CoERT 1760 SECONDS 
10Z5 FORRATCIEOx,' 4 3 BAUDSUENGTH (CDE) 1/30 SECONDS 
1026 PORMATCIUX, * Sa: JPAUOSCENGTRSCIELIO 1715 SECCHD Sas 
LOZ? FORHAICZ, 2X, 

% ‘ENTER THE NUMBER OF BAUDS I! PACKET SNUB ERS Phas a. 
1030 FORMATC/Z, 2x, “SHOULD YOU LIKE THESDTES GEE cae TER TRANSFORM We 

4 2X%,' OF THE OUTPUT SIGNAL 2 CENTER es sO) OCD eee 
1040 FORHATCZ, 2X, 

% ‘WOULD YOU LIKE THE DFT OUTPUT WINDOMED Se es 

i AEE Rey E Se OK 0 One ee 
1050  FORIVAT GAA Zx. 

% "PLEASE ENTER THE DESIRED TNPUI HILDE VAD 

ras SIGHAL>=TO-NOISE. RATIO DOE ee 
1100 FORMAT C7, 2x, 

% *THIS PROGRAM ENCODES A QPSK MULTIFREQUENCY SSG hs 

4 * JHE PHASES ARE SHONIN BELOW FOR CHES FREQUEI(G ea 
L101 FORMAICZZ, 

1A Gee eres es ae 

4 6%," ¥ ire ee ¥ ae ger 


) 


Vw SY 


oe 


7» 


7,2X, 
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LOK 2 ra ae 1 te 
4 ions ' ay as eo ies 
YG mera ga es 
a 5x a: tee a, 
LOK a ee ore LL 
ZZ 6X58 ay ted Sens 35 
me GX By eke Uae en ae 
ie) © Ua Te 
Y Goce x Uy oat ar Ve, 
DG oa 5 ee 4 ',7, 
“6 Mae ee to ',///) 


1102 FORMAICZ7, 2x, 
ioe eee OnE er the FORLOUTHG METHODS FOR ENCODING *,7,4%, 
Ze hee enoes FURTEACH FREQUENCY FOR EVERY BAUD °,7,4%, 
ZA Mm eACH PACKET ....7) 

OS: EURMAiGz ox ENTER ...%n7, 7X, 
oe eile ROoRA RAIDOMEY SELECTS ALL THE PHASES °,7,7x, 
ea Poke i, 19," PACKETS °,7,7%, 
coe OUumeyolvI DUALLY SELECT THE PHASES *) 

1110 FORMAT(C//, 2x, 
A Deeeer ONE OF THE FOLLOWING METHODS FOR ENCODING ",7,4X, 

Tose, 


A) ineeenases FOR PACKET NUMBER: ',14,! 
Pri lr ORMAD Czy ox CENTER ...'o 757%, 
% "Ll: THE PROGRAN RANDOMLY SELECTS AUU THE PHASES) *,75 7X, 
aoa Siinwrecke! HUMBER: ',14,7,7X, 
oe eee DUALLY SS EEECT THE PHASES ',7,7X, 
roy DORR GOCKET HUMBER: *, 14) 


120° FORMAT 77; 2X, 


Ze oe eeeeOWE SUES THE FOMBROIING HETHODS FOR ENCODING *,7,49X, 

Ze eee semen OoR PACKET MHUMBER: *,19,* BAUD NUHBER: *,14,' 
Weel ORM ieee tee EHIER ...'57,7X, 

7 emir rem RAl RANDONULY SELECTS ALL THE PHASES *,7,7%X, 

ae FOR PACKET NUMBER: ',149,* BAUD NUMBER: ',149,/7,7X, 

i eum LYI DUALLY SELECT THE PHASES *,7,7X, 

oe FOR PACKET HUMBER: *,19,° BAUD HUMBER: ',14) 


0 (PORMAnGA7 2x, 
Valen eenrerR the DESTRED QUADRANT FOR THE PHASE *,7,49xX, 
Zo UG ere ieniUinbeER: ',14,° BAUD NUMBER: *,19,7,9X, 
Pe eREQUrid@ nm enUlBER €K)J: *)14," s..°) 


2000 FORMAT(C18) 

2001 FORMATCA38) 

2010 FORMAT(C2X,16,2X,16,2X,F10.3,2X,149) 

eve FORMAICZ7, 2x, 138) 

2100 FORMAT(2%,16,49%,F19.9) 

2oUOmrPORMAIC 77 ex, DET LNPUT DATA FOR BAUD" ® *,149;,7, 5x, 
Calo oe een yok,  TNAG PART") 

210m ORMALC2x > 1S. 2xX,El2.6,2X,E12.6) 

ZecverORMelG77 ex, DEY OUTPUT DATA FOR BAUD © ',14,7,5X, 


ewe ReMi AR I 7745 TIAG PART", 7X, *HAGHITUDE',16X,*PHASE',7 


“4 649X%,'CDEG)',11X,'CRAD)') 
ao) SOUR ebee we eal 9. 8,2X,E149.8,2X,E14.8,49X%,E19.8,2xX,E14.8) 
2500 FORMATO’  FORS SAUD NUMBER *,14,° $*) 


IS 

c PLOTTING 

C 
POH eS=IITCXXHIN)D = 0.2 
Pre = LNT CXXMAX? + 0.2 
Penne Ler yt) =- 0.2 
Rogan = SLUG YYTIAX) +.°0.2 
COU ee ela Ce2ZhilH) = 0.2 
COWAX = HEC ZZNAX) + 0.2 
See ENPCCCHIN) — 1.0 
CCiax 8= VHP CCCHAX) + 1.0 
Ret = IT ORRIN) -— 1.0 
RRWAX = INTORRITAX) + 1.0 
AER IH = VHICALFHNITID - 0.2 
CUA X= aT CALFIIAX) + 0.2 
THTIMIN = INTCTHTHIHD - 0.9 
TH Ax = TNC THTHAX) + 0.4 


CED el etary CALL COMPRS 
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CALLYPAGCECI 17 25) 
CALL NOBRDR 
CALL AREA2D(8,6) 
CALL XHANEC. TENE Ch)? aap 
CALL YHAMEC® HAGHITUDE.$ 'y oon 
CALL HEADINC' RECEIVED SIGs. VUtre 
CALL GRAF GO, SCALE Ginihi Ss, RXMIN, *SCAL ES 
CALL GRID(),1) 
CALL SETCLR('GREEN') 
CALL CURVEC TTHESYYRX,HPISs, 0) 
CALL ENDGR GOD 
CAUL ENDPIECG? 
IF (BDIOIE 2Gleo awe 
CALL PAGE CII, a 250 
CALL NOBRDR 
CALL AREAZDCS; 6) 
CALL. XNABIEC*— “BAUD CIL) => "71005 
XXM = XXMAX - XXMIH 
IF OXAM 2. GEO Gre) he , 
XXMIN = XXMIN 7 1 
XXMMAX = XXMAX 7 1 
CALL YHANECGC 9 xXGL 
DO 900 1 = I, BD 
ACID -2ACP) 7 A 
CONTINUE 
Eboe 
CALL YHAMEC*® XCEL SCE $',100) 
ENDIF 
CAUL HEADDUIC Ss X-PUs righ 7 'ohO0, 3,19 
CALL GRAF(1], "SCALE", BDIOTL, XxMI SCA EE ae 
CALL GRID(1,1) 
CALE SETCER Gre ine) 
CALL CURVECL BAUD A BDIOTL 
CALL ENDGR CO) 
CALL ERUe Cd) 
CALL PAGEG@I1 73.050 
CALL NOBRDR 
CALL AREA2D(8,6) 
CALL XPANEC” “BAUD CLUE )mo* 3 10) 
YUN = erhiaxX Se a 
LPO Yi GE a D0 0) be 


gl) 
»>RXHAX) 


YOLEN = Vries 000 

YYMAX = YYMAX 7 1000.0 

CALL YHANEC Byes Chel $',100) 

DO 905 lt -> Il, oBDIGIE 

Y¥CJO = 40) 7 TU00Re 

CONTINUE 
ELSE 

CALL. YHAVIEC* TCI. Cri.) $*',100) 
ENDIF 
CALL HEADUING" YoROstT ron op Urol 


CALL GRAF(1, "SCALE", BDIOTL, YYIITH SCARE hr ee 
CALS GORI DC1s 12 

Call SETCER CC MAGERT SE) 

CALL CURVECLBAUD,Y,BDTOTL,0O) 

CA ER SEHDGR CG) 

CEUR Ero Co) 

CALISEACE CI I.6. 5) 

CALL NOBRDR 

CALL ARENCDCS, 62 

CAUL  XHAMEC’ “BEUD CLE) 37> L0Ge 
EZNige= a Cel lnn ss CO 

Lee CZ Ce GE) OUT Oar aE 


ZEN. = 2c 7 100020 

Z20AX = ZZMAX 7 1000.0 

CAD HAVE See2 Gl eek rie $',100) 

DO.9)0 1. ="1, BOPOTL 

ZOE 2= 201) Zevocoet 

CONTINUE 
ELSE 

CALL YHAMEC’ ZCLTISICE I? > > VOD 
ENDIF 
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915 


920 


TX 100,351) 


CxEL HEADING t Z-POSITION et UCs oD 
Sabb GRAFCIT; "SCALE", BDTOTL, ZZ, "SCALE*, ZZMAX) 
CALL GRID(1,1) 
eye SelerR Cy BLUE * ) 
CALL CURVECLBAUD,Z,BDTOTL,O) 
CALL ENDGR(O) 
CARL ENDPLUCO) 
Cab PAGEC11,8.5) 
CALL NOBRDR 
CALL AREA2D(8,6) 
CAEE AWAMEC® BAUD (LL) $*,100) 
Cen = -Cenax = CCMIT 
Meee CAxe Ce. 100.0) THEN 
CChiItt = CCHIN 7 1000.0 
COMAX = CCHMAX 7 1000.0 
PoeecninAnew: | CCLE) SGmET/ SEC) $*',100) 
DOS tose - 72, BDIOTL 
C@lpe=—C(l> 7 1000.0 
CONTINUE 
ELSE 
CAP eeaueanieG” ~COLE) CPiZASEC) 7'7 100) 
ENDIF 
CAbeteedinC’ SPEED OF SOUND $°,100,3,1) 
Chee CRAECI. *SCALE*, BDIOIL,CCMHIN, "SCALE", CCHAX) 
CALL GRID(1,1) 
GAPE VSETCERC * RED") 
CALL CURVECLBAUD,C,BDTOTL,0) 
CALL ENDGRCO) 
Peer DE iC 0) 
Coble rAGECI 1,8 .5) 
CALL HOBRDR 
CALL AREA2D(8,6) 
CALL AWAHEC*® BAUD (LL) $*,100) 
RRM = RRHAX - RRMIN 
Pee okRi -CE. 100.0) THE 
RRMIN = RRMIN 7 1000.0 
RRMHAX = RRMAX 7 1000.0 
Camb yitanlec": ROLL) CKFT) $*',100) 
DUS ecu sie. ky BDIOIL 
Polvo € 1) 4. 2.000..0 
CONTINUE 
EESE 
Cle tAMEC’ “RELL ) CFI) 7 00) 
ENDIF 
CALL HEADINCG* SLANT RANGE TO RECEIVER $',100,3,1) 
Cem ernnGin, scnLe*, BN DIOTLZRRIIN, *SCALE*, RRMAX) 
CAULEY GRIDC], 1) 
Ceo ClR< GREEN) 
CALL CURVECLBAUD,R,BDTOTL,0O) 
CALL ENDGR(C0) 
Call ENDPLCO) 
CrP PAGEC 11,3 .5) 
CALL NOBRDR 
CALL AREA2D(8,6) 
CALL XNAMEC' BAUD (LL) $*',100) 
Combe iteinie« os, ALPHAC(LLJ 9 3*,100) 
Cen aimtiC COMPRESSION FACTOR DUE 10 THE HOVING 
Come Cinihcl. oCALE*, BDIOTLZALFININ, ‘SCALE*, AL FIAX) 
Ceri Or 1 DC lie) 
Clie cEICEKC RED ) 
CALL CURVECL BAUD, ALPHA, BDTOTL,0) 
CALL EUDGR GO) 
CALL ENDPL(0O) 
Cobb ACEC], 8.5) 
CALL NOBRDR 
CALL AREA2D(8,6) 
Cee AAMEG* BAUD (CLL) $*,100) 
CABLES THANE’ “THETACEL) CDEG) $*,100) 
Cee HEADING ANCHE BETWEEN ROLL) AND 20 $',100,3,1) 
ChEU Ghee ti SCALE -sDDTOTESTAINII, *SCALE’, THTMAX?) 
CALLeGRIDGCI 1) 
Geely seICeRC* BLUES) 
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CALL CURVECLBAUD, THETAD, BDIOTL, 0) 
CALL ENDGR(0) 
CALL ENDPECOD 

ENDIF 

CALL DOHNEPL 


9999S Ot 
END 


SUBROUTINE GAUSS(CDSEED, AVG, VAR, ZRHD) 


THIS SUBROUTINE GENERATES A GAUSSIAN RANDOM NUMBER HITH 
MEAN = AVG AHD VARIANCE = VAR 


INPUTS: DSEED = A DOUBLE PRECISION SEED THAT SHUS Tepe ey nee 
AVG = THE NEAN OF THE GAUSSIAN RANDOM VARIABLE 
VAR = THE VARIANCE OF THE GAUSSIAN RANDOM VARIABLE 


OUTPUTS: DSEED = THE SEED IS CHANGED DURING EXEGUIT aerhie hee 
REQUIRED FOR JHE HEX] SUBROWTIVESe AE 
ZRHD = A GAUSSIAN RANDOM NUMBER WITH MEAH = AVG AND 
VARIANCE = VAR 


REAL AVG, VAR, ZRND 

DOUBLE PRECISION 22." DSEED? URIID 

ZZ 2 0.0 

DO 000 -f =. 1, 22 
DSEED:=-DSEED + 2357954571597 557 
URHD = GGUBFSC(CDSEED) 
LA = CURD =" 0 Spee 

5000 CONTINUE 
ZR HSU LZe * (VAR ASO. 500 AVG 


OOOO OOOOCO Teele ee 


RETURH 
END 
C 
SUBROUTINE PHASECDSEED, RNDEHIL, lei 
C 
C THIS SUBROUTINE SELECTS A PHASE RANDOMNLEY FROH QUADRA Tews) Ome 
C 
C THe US: DSEED = A DOUBLE PRECISION SEED THAT MUST BE A VARTABE 
C 
é OUTPUTS: DSEED = THE SEED IS CHANGED DURINIG EXECUTIONS INT Crimes 
c REQUIRED FOR THE HEXiTSSUBEUUT Tie alan 
G RNDPHI = A RANDOM PHASE FROM ONE OF THE FOUR QUADRANTS 
C IN RADIANS 
C IPHI = THE QUADRANT NUMBER OF THE RANDOM PHASE 
cS 
REAL PHIRIDS RNUDPHI 
DOUDECE PRECIS On. DSEED 
INTEGER STP AI 
C 
Pl =) 4300 ATANC)] 0) 
PHIRHD = 0.0 
{Pal = "0 
PHIRWD =-GCUBFSCDSEED) 
DSEED = SDSEED + bas. cos 
FHIRHD = SGRHIRND. = -4., 00 
TER = Se RD 
IF CiPlil JEG. lL) RUDPHI = (95.07 * Pier = 
IF CIPHI .E@. 2) RHDPHI = (135.0 0.%5) 10 ee 
IF CIPHI .E®. 3) RHDPHI = ( [ee Pie 130. 0 
IP SCIPH] .EQ@. S42) (RMNDPHI = (-45.0" % Pile Tove 
PURE 
END 
C 
SUBROUTINE. DEVGHITT = 270 0U I 
c 
c THIS SUBROUTINE COMPUTES THE DISCRETE FOURTERR IRA EURO 
C A COMPLEX DATA SET OF W C= NH14+1) POINTS STORED Timi teak 
c JHE RESULT IS STORED Il! THE COMPUEX SARA rT : 
C 


COMPLEX XINCOzNMY), XOUTC Os Tia 
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©) C) CY ee 


Hos NM) + 1 

PI = G4.0XATANC] .0) 

EN = N 

Hi2 = WM1l - | 

Penk EQ. 02" THER 
XOUTCO) = X1HCO)d 


BETA = 2,0 5% PI 7 EN 
Ho= CMPLXCCOSCBETA), ~-SINCBETA)) 
DO 3110 K = O, hl 
HI = bIX*K 
XOUTCK) = XINCHMI) 
DOU S100) 6k = M2, 60, 1 
XOUTCK) = XOUTCK)¥IM + XINCL) 
CONTINUE 
CONTINUE 


ENDIF 

RETURN 

END 

SUBROUTINE VLINECJ1, JN, XARRAY, B) 


REAL XARRAY(0:4096), BC0:4096) 


DRAM VERTICAL LINES 


DOV SeUU JUINE = Ji, J 
XFROM = XARRAYCJLINE) 
ALOy = XEROH 
YFROM = 0.0 
TLO, = BCJLINE) 


ChbpreRevECCXEROM, YFROM,XTO,YTO, 0000) 


S2UO CONT EHUE 


RETURN 
END 
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